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INTRODUCTION 


1. CONSIDÉRATIONS PRÉLIMINAIRES 


L'évolution de l’informatique conduit naturellement à l’utilisa- 
tion de l’ordinateur pour la résolution des problèmes de Résistance 
des Matériaux. Le présent a pour objectifs : 


— Rappeler quelques résultats pragmatiques concernant les 
structures, sans revenir sur les théories de base qui sont exposées 
dans de nombreux ouvrages. 


— Présenter les modules de calcul correspondants dans un 
langage de programmation facilement adaptable sur différents 
micro-ordinateurs. 


C’est aussi une initiation au calcul automatique des structures, 
discipline qui a eu un essor particulier dans les dernières décennies, 
et qui est au confluent de la Résistance des Matériaux, de l’Analyse 
Numérique et de l’Informatique. Le développement de la micro- 
informatique a rendu possible la solution de problèmes qui nécessi- 
taient auparavant des moyens puissants et le calcul des structures va 
se vulgariser. 


Si la majorité des projeteurs utilisera des programmes achetés 
à des sociétés de services, dont le coût ne fera que décroître avec le 
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temps, il est néanmoins souhaitable de savoir sur quels principes 
sont bâtis ces programmes, ce qui permet d'éviter bien souvent 
des erreurs lors de leur exploitation. Il est aussi parfois utile de 
savoir écrire un programme spécifique mieux adapté à ses besoins 
ou à ses habitudes de travail. Le présent opuscule pourra alors être 
utile. 


2. CHOIX DES PROBLÈMES TRAITÉS 


Nous avons choisi quelques problèmes classiques fréquemment 
rencontrés en Génie Civil, et dont la solution est donnée par la mé- 
thode des matrices de rigidité. 


1. Le calcul des ossatures planes ou spatiales. Ce sont des struc- 
tures composées de barres travaillant uniquement à l’effort normal. 


2. Le calcul des poutres continues. Le programme est déduit du 
programme d’ossatures avec un minimum de modifications. 


3. La résolution des problèmes d’élasticité plane à l’aide d’élé- 
ments finis triangulaires (les plus simples). 


4. La répartition des contraintes de cisaillement et l’inertie de 
torsion pour une barre soumise à un moment de torsion. Ce pro- 
blème se ramène à l’étude d’une structure grâce à l’analogie de 
membrane. Il est traité à l’aide d’éléments triangulaires. 


Le lecteur ayant bien assimilé l’organisation générale du pro- 
gramme pourra l’adapter sans trop de difficultés aux portiques 
plans, aux portiques spatiaux, aux structures plan-spatial, aux dalles 
et aux coques. 


Il pourra également développer les programmes présentés 
suivant deux voies diamétralement opposées : création d’un super- 
programme unique, pour intégration de tous les modules nécessaires, 
ou au contraire simplification et spécialisation pour des usages bien 
précis. Par exemple, un programme de calcul d’ossatures ne traitant 
qu'un seul cas de.charge constitué par des forces nodales comportera 
beaucoup moins d'instructions et pourra résoudre des structures 
plus importantes. 
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3. PROGRAMMATION 


Les programmes présentés sont écrits en BASIC pour micro- 
ordinateur HP 9845. Nous avons recherché une présentation aussi 
claire que possible, étayée par de nombreux commentaires. Il en 
résulte en contrepartie des programmes volumineux. Pour des raisons 
de lisibilité, nous avons utilisé des variables en tableaux portant des 
noms mnémoniques. Cette façon de procéder est peu rationnelle 
quant à la bonne utilisation de l’espace mémoire de l’ordinateur : un 
problème de structure dépend de plusieurs paramètres (nombre de 
nœuds, nombre de liaisons, nombre de chargements différents, etc.), 
et on réalisera un gain appréciable de mémoire en rangeant toutes 
ces variables dans un seul tableau géré à l’aide d’une table d’adresses. 
C’est l’allocation dynamique de mémoire, systématiquement em- 
ployée dans tous les grands programmes de calcul disponibles sur le 
marché. 


Il existe de bonnes raisons de penser que le coût de la mémoire 
disponible ira en diminuant et la vitesse d’exécution, en augmentant. 
Aussi n'est-il pas rentable de se livrer à des astuces de programmation 
si, pour une poignée d’octets supplémentaires et quelques secondes 
de plus, on obtient un programme plus lisible et plus facilement 
modifiable ou adaptable, ce qui se traduira par une économie finale. 


La convivialité d’un programme, en particulier la facilité d’intro- 
duction des données et la présentation des résultats, conduit à dou- 
bler ou à tripler le volume des instructions strictement nécessaires à 
la solution du problème technique. Il n’est pratiquement pas possible 
de décrire une structure autrement qu’à l’aide d’un fichier facilement 
accessible et modifiable. Cet aspect n’a pas été abordé dans l’ouvrage, 
et nous nous sommes contenté d’une introduction des données 
toute théorique. 


Nous avons évité l’usage des mémoires auxiliaires (disques ou 
bandes) de façon à ne pas particulariser les instructions utilisées. 
Celles-ci ont été réduites au strict nécessaire, et les correspondances 
avec les autres langages de programmation sont évidentes. 


L’emploi de variables mnémoniques nous a parfois conduit à 
donner deux noms différents à une variable suivant la destination 
du programme : il est naturel de parler du ‘nombre de barres’ pour 
une ossature et d'introduire la variable ’Nbarres', alors que la même 
variable employée dans le programme d’élasticité plane est désignée 
par ‘Nelem’ (nombre d’éléments). 
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Suite à une habitude héritée du FORTRAN, les variables 
relatives à des nombres entiers commencent par I, J, K, L, M, N. 

Les numéros des lignes des sous-programmes dont les listings 
sont reproduits dans le reste n’ont qu’une signification locale et 
évitent l'emploi des étiquettes. 


4. ORGANIGRAMMES 


Les organigrammes présentés reflètent l’enchaînement des 
instructions et permettent de mieux appréhender l’organisation 
logique des sous-programmes. 


Nous avons utilisés les symboles suivants : 


— Bloc d’instructions arithmétiques 
ou d’affectation 
— Appel à un sous-programme. 


— Boucle. La variable I prend successivement 
les valeurs 1 à N. 

I=N,1,— 1 signifie que I prend 
successivement les valeurs N, N — 1,1 


— Test 


— Sortie d’un sous-programme. 





Nous avons également utilisé des instructions relatives à des 
matrices. Celles-ci sont précédées de l’abréviation MAT. 
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1 PRÉSENTATION 


1.1- RAPPEL SOMMAIRE DE LA MÉTHODE DES DÉPLACEMENTS 


Toute structure peut être considérée comme formée des sous- 
structures élémentaires. Pour les ossatures, les portiques, la sous- 
structure élémentaire est un élément linéaire dont la définition 
s'impose à partir de la géométrie. Pour les structures de type plaque 
ou coque, les sous-structures sont des petits éléments découpés 
presque arbitrairement dans le milieu continu. Nous utiliserons 
des éléments de forme triangulaire. Dans tous les cas, on établit une 
relation entre les forces appliquées en certains points (nœuds) de la 
sous-structure élémentaire et les développements qui en résultent. 


Ces déplacements vont constituer les inconnues du problème. 
En exprimant que chaque nœud est en équilibre sous l’action des 
efforts induits par la déformation des éléments incidents et la charge 
directement appliquée, on constitue un système d’équations linéaires 
dont la solution fournit les déplacements des nœuds. On utilise 
alors ces déplacements pour calculer les efforts dans les sous-structu- 
res élémentaires et ensuite les réactions d'appui. 


1.2 - DEGRÉS DE LIBERTÉ 


C’est le nombre de paramètres nécessaires pour définir le dé- 
placement d’un nœud. Ce nombre dépend du type de la structure 
étudiée. Il sera désigné par ‘Ndl’. Donnons quelques exemples : 


1. Ossature plane 











Fig. 1. - 2 degrés de liberté par nœuds : 
déplacement horizontal u 
déplacement vertical v. 


2. Ossature spatiale 





Fig. 2. - 3 degrés de liberté par nœudés : 
déplacement horizontal u 
déplacement vertical v 
déplacement horizontal w. 


3. Poutre continue (charges perpendiculaires à la fibre moyenne) 


P 
Dy 





@ 


Fig. 3. - 2 degrés de liberté par nœud : déplacement vertical Dy, rotation 6. 


4. Elasticité plane 





Fig. 4. - 2 degrés de liberté par nœud : déplacement v, déplacement v. 


5. Analogie de membrane (problème de torsion) 


Z 


X 


Fig. 5. - 1 degré de liberté par nœud : déplacement vertical w de la membrane. 


1.3- DÉFINITIONS RELATIVES AUX NOEUDS 


1. Nombre de nœuds par élément : c’est le nombre de points de 
connexion de l’élément avec les éléments voisins ou avec le milieu 
extérieur. Ce paramètre est désigné par ‘Nnelem’ dans les program- 
mes. Pour les éléments de type linéaires utilisés dans cet ouvrage (os- 
satures, poutres), Nnelem=2. Pour les éléments triangulaires (élasti- 
cité plane, torsion), Nnelem=3. 


Le produit (Nnelem X Ndi) donne la dimension de la matrice 
de rigidité de l’élément. 


2. Nombre de nœuds de la structure : c’est le nombre total de 
points de la structure pour lesquels on va déterminer les déplace- 
ments. Il est désigné par 'Nnœuds’ dans les programmes. 

Le produit (Nœuds X Ndi) donne la dimension du système 
linéaire à résoudre. 


3. Définition des appuis : ce sont les liaisons de la structure 
avec le milieu extérieur, constituées par des nœuds dont on connaît 
certains déplacements qui peuvent être nuls ou non (déplacements 
imposés). Les degrés de liberté de chaque nœud sont spécifiés par un 
code : ‘Isup’. 


Nous avons choisi le code suivant : 


1 si le déplacement est bloqué 
0 si le déplacement est libre. 


Pour un seul degré de liberté, on a Isup=0 ou Isup=1 


Pour deux degrés de liberté par nœud, par exemple u et y, 
on peut avoir les combinaisons suivantes : 






Libre 
Bloqué suivant Y 1 
Bloqué suivant X 10 
Bloqué suivant X et Y 11 


Pour trois degrés de liberté par nœud, par exemple , » et w : 












Libre 0 0 0 0 

Bloqué suivant Z 0 0 1 1 
Bloqué suivant Y 0 1 0 10 
Bloqué suivant X 1 0 0 100 
Bloqué suivant Y et Z 0 1 1 11 
Bloqué suivant X et Z 1 0 1 101 
Bloqué suivant X et Y 1 1 0 110 
Bloqué suivant X Y Z 1 1 1 111 


Les appuis doivent obligatoirement assurer la stabilité de la 
structure vis-à-vis de tout chargement compatible avec la définition 
de celle-ci. Dans le cas contraire, la résolution du système linéaire ne 
pourra être menée à bien. 


1.4- DÉFINITIONS RELATIVES AUX ÉLÉMENTS 


1. Nombre d'éléments : c’est le nombre total de sous-structures 
élémentaires constituant la structure globale. Nous supposerons dans 
chaque programme un seul type d’élément, mais il est évidemment 
possible d’en avoir plusieurs. On pourra par exemple créer des élé- 
ments linéaires dans le programme d’élasticité. 


Les grands programmes généraux possèdent une bibliothèque 
d'éléments permettant de modéliser un grand nombre de structures 
différentes. Le nombre total d'éléments est désigné dans les program- 
mes par ‘’Nbarres’ pour les ossatures et les poutres continues, et par 
‘Nelem' pour l’élasticité plane et la torsion. 


2. Définition des éléments : les éléments sont définis par les 
numéros des nœuds placés dans le tableau ’Liais’. Pour les éléments 
triangulaires, il est particulièrement important de définir les nœuds 
sommets dans le sens trigonométrique positif (sens contraire des 
aiguilles d’une montre). Dans le cas contraire, la surface du triangle 
est négative, et le pire est à craindre ! 


Le tableau ‘Liais’ contient également les références de certains 
attributs de l’élément tels que la section et le module E pour les 
ossatures. 


1.5- LE SYSTÈME LINÉAIRE GLOBAL 


Chaque ligne du système linéaire exprime l’équilibre d’un nœud 
pour un degré de liberté donné. Plutôt que générer le système linéaire 
ligne par ligne (très peu logique, car il faudrait rechercher l’occurren- 
ce d’un nœud parmi tous les éléments), on le forme globalement 
en cumulant les coefficients de chaque matrice élémentaire à la 
bonne place. On sait que la matrice du système linéaire est symétri- 
que. Afin d'économiser la place en mémoire, on ne conserve que les 
coefficients situés au-dessus de la diagonale. De plus, nous verrons, 
lors de la description de l’algorithme de résolution, qu’il n’est pas 
utile de conserver les coefficients nuls situés avant le premier terme 
non nul de chaque colonne. Le gain de place (et de temps de calcul) 
ainsi réalisé permet d’accroître les possibilités du programme, mais 
implique une organisation des calculs plus sophistiquée. 


1.6 - MÉMOIRE AUXILIAIRE 


Nous avons choisi de conserver en mémoire un certain nombre 
de valeurs caractéristiques des éléments afin de ne pas avoir à les 
recalculer. Le rôle de cette mémoire de masse est tenu par le tableau 
’R'. Les micro-ordinateurs disposant d’une large mémoire disque 
pourront Sauvegarder tous les coefficients des matrices de rigidité 
des éléments. Notons que ces matrices sont utilisées deux fois : 
pour former le système linéaire et pour calculer les efforts à partir 
des déplacements. 


1.7- CHARGES EXTÉRIEURES 


Nous avons envisagés plusieurs types de chargement en fonction 
de la destination des programmes. 


1. Ossatures planes et spatiales : 
- charges appliquées aux nœuds 
- poids propre 
- température 
- déplacements d’appuis. 
2. Poutres continues : 
- charges appliquées aux nœuds 
- charges concentrées sur les barres 
- charges réparties sur les barres 
: déplacements d’appuis. 
3 - Elasticité plane : 
- Charges appliquées aux nœuds 
- poids propre 
- température 
- déplacements d’appuis 
4. Torsion (analogie de membrane) : 
- pression uniforme, générée automatiquement. 


Les données relatives aux charges ne seront pas conservées en 
mémoire : elles sont reproduites immédiatement après introduction, 
les effets correspondants sont calculés et cumulés dans les seconds 
membres du système linéaire. 

Cinq cas de chargement distincts ont été prévus : ils peuvent être 


constitués d’une combinaison quelconque de types de charges (char- 
ges aux nœuds, poids propre, température ou déplacement d’appui). 


LES SOUS-PROGRAMMES 
GENERAUX 


2.1 - LE PROGRAMME PRINCIPAL 


A des modifications mineures près, les quatres programmes 
présentés sont construits de manière semblable. 


Le programme principal comporte : 
— Des instructions de déclaration telles que 
OPTION BASE 1 signifiant que 1 sera le premier indice 
d’un tableau. 
PRINTER IS 7,1 spécifiant l’unité d'impression. 
— La description des variables principales et des variables 
auxiliaires, ainsi que la déclaration de ces variables en ‘commun’, afin 
d'êtres accessibles par différents sous-programmes. Les tableaux 


recevant des nombres entiers sont déclarés de type INTEGER afin 
de gagner de la place en mémoire. 


— L'initialisation. des tableaux en ‘commun’ au cas où l’on 
enchaîne un autre problème. 
— Les sous-programmes, qui sont de deux types : 


e les sous-programmes généraux (indépendants du type de 
structure étudiée), dont les arguments sont transmis en paramètres, 


e les sous-programmes spécifiques au type de structure étudiée, 
dont les principaux arguments sont dans les zones en commun’. 


Lecture des données 


Indapp Indices des degrés de liberté bloqués 


Reproduction des données 


marres) Koucle sur les barres (ou éléments) 


Elements de la matrice de rigidité 
Hauteur des colonnes 


Indices des termes diagonaux 


Assenb Assemblage de la matrice globale 


Foucle sur les barres (ou éléments) 


atrig Matrice de rigidité de la barre 


Calcul des adresses des coefficients 


Cuaul dans la matrice globale 


Second Formation des seconds membres 


Generetenp Forces nodales dues à la température 


Forces nodales dues aux déplacements d'appuis 





Modification en fonction des appuis 


Résolution du système linéaire 


ur ) Houcle sur les cas de charges 


Impression des déplacements 


Calcul et impression des efforts et des réactions 


Fig. 6. - Organigramme général 
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L’organigramme du programme principal est donné figure 6, 
Les sous-programmes généraux décrits dans le présent chapitre 
sont encadrés d’un trait fort. Le listing du programme général, par- 
ticularisé pour chaque problème traité, sera reproduit dans les 
chapitres correspondants. 


2.2- RANGEMENT DE LA MATRICE PAR LIGNE-DE-CIEL (SKY-LINE) 


Avant de poursuivre la description des programmes généraux, il 
convient de préciser la manière dont sont rangés les coefficients 
de la matrice du système linéaire. La méthode utilisée est appelée 
‘ligne-de-ciel’ et consiste à ranger dans un tableau à un seul indice 
les coefficients situés au-dessus de la diagonale, en ne tenant pas 
compte des termes nuls situés au-dessus du premier terme non nul 
d’une colonne (ce qui définit une ‘ligne-de-ciel'). Cela est illustré 
figure 7. 








Fig. 7. - Rangement d'une matrice symétrique par ligne-de-ciel. 
1 - Tableau complet à deux indices 
2 - Tableau à un indice. 
3 - Täble des adresses des termes diagonaux. 


Les coefficients seront rangés dans le tableau ’A' et les adresses 
des termes diagonaux dans le tableau "Maxa'. Pour définir ce tableau 
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d’adresses, on effectue un premier balayage des éléments au cours 
duquel on détermine les hauteurs de colonnes (ainsi que les caracté- 
ristiques des éléments qui serviront à former sa matrice de rigidité). 
Par cumul des hauteurs de colonnes, on obtient les adresses des ter- 
mes diagonaux. 


2.3 - CALCUL DES INDICES DES DEGRÉS DE LIBERTÉ BLOQUÉS 


Oui: le noeud n'est pas uñ appui 














Nombre=Isup(1) 


Recherche des degrés de liberté 
bloqués du noeud 


CALL Decode 
LENdLKCI-1) 


Boucle sur les degrés de liberté 






Numéro du K ème degré bloqué 


Ê 
K=K+i 
Lappui(K)=L 


Ndib : Nombre total de degrés 
de liberté bloqués 


Fig. 8. - Organigramme ‘indapp’ 
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Il est pratique, pour la suite des calculs, de savoir si un degré de 
liberté est libre ou bloqué. Cela est effectué par le sous-programme 
’Indapp', qui génère à partir du tableau des codes appuis ‘Isup' le 
tableau auxiliaire ‘Iappui’. Ce travail peut être également exécuté au 
moment de la lecture des données concernant les nœuds. Rappelons 
à ce sujet que, dans le programme bien connu "STRESS, on spécifie 
le nombre de supports en données générales, on désigne les supports 
par la lettre ‘S', et on précise ensuite si il y a des possibilités de 
déplacement dans certaines directions. En informatique, le même 
résultat peut être obtenu par bien des moyens différents. L’organi- 
gramme du sous-programme ‘Indapp' est donné figure 8. 


Listing des sous-programmes ’Indapp' et 'Decode’ : 


2290 1 XHPHNNNNENHEHNNNNHNNNNNHNMENHNMNNNRNNRHNNMNERNNNENNNRHNIANNRNNNNRNX 
2300 SUE Indapp(Nnoeuds,Nd1,Ndlb,INTEGER Isup(k),Lappui(*)) 

23140 1 KA HEMENANHNNNNNMENNMNNNNRNNRNNNNNNNNMNNNNRNERNNNNNRNNNRX 
2320 ! 

2330 ! Appelé par Prog. Princ. 

2340 ! Appelle : Decode 

2350 ! 

2360 OPTION RASE i 

2370 ! 

2380 ! Sous-progranme créant le vecteur Iappui des degrés de liberté connus 
2390 ! 

2400 DIM Coef(3) 


2410 K=0 

2420 FOR 1=1 TO Nnoeud5 

2430 IF Isup(1)=0 THEN 2530 
2446 Nonbre=Isup{(1) 

2450 CALL Decode(Nombre,Nd1 ,Coef(*)) 
2460 LENd1K(I-1) 

2470 FOR J=1 TO Ndl 

2480 LEL+i 

2490 IF Coef(J)=0 THEN 2520 
2500 K=K+1 

2510 Jappui(K)=L 

2520 NEXT J 

2530 NEXT I 


8540 Ndlb-K 
2550 SUBEXIT 
2560 SUHEND 
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2570 
2580 
2590 
ce60 
2610 
2be0 
cb30 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
e7e0 
2730 


LR OHHENNRHNHNNHNRNRHHHRAHNRNNNNRNNRANONNRRRINNNONENNRNNRPHNNNRNINRHIONNNERONX 
SUE Decode(Nombre ,Nchifr ,Coef (#)) 

LR OBHNNNACNNNANENHNNNNNNNNMIONRENNONNRRREONARPONNNRRRINNRNRNONX 
} 


| Appelé par Indapp 
l 


OPTION BASE 1 

! 

1 Ex : Nchifr=3; Nonbre=110 -----) Coef(1)=1, Coef(2)=4, Coefi3)=Ù 

1 

FOR I={ TO Nchifr 
N=INT(Nombre/10) 
Coef(Nchifr+{-I)=Nombre-{0xn 
Honmbre=N 

NEXT I 

SUREXIT 

SUBEND 


2.4- CALCUL DES HAUTEURS DE COLONNES 
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On détermine les hauteurs de colonnes à partir de la définition 
des éléments. L’algorithme peut être résumé de la manière suivante : 


Pour chaque élément : 
— Déterminer le plus petit parmi les nœuds définissant l'élément. 
— Pour chaque nœud de l’élément : 


— Déterminer la hauteur de la colonne du premier degré 
de liberté. 

e Si cette hauteur est inférieure à la hauteur déjà déter- 
minée pour ce degré de liberté, passer au nœud suivant. 
e Sinon, pour chaque degré de liberté du nœud, réactua- 
liser les hauteurs de colonne correspondantes. 


— Passer au nœud suivant. 


L’organigramme correspondant est représenté figure 9. 


Mini=9999 











Nhaut=(Liais(l,J)-Mint)KNdl+i 


Li=(Liais(l,7)-1)XNd1#1 


Maxa(LitL-1)=Nhaut+L-{ 







Recherche du noeud ayant le plus 
petit cardinal 


Boucle sur les noeuds de l'élément 


Hauteur de la colonne pour 
le noeud considéré 
Degré de liberté correspondant 


Comparaison avec la hauteur actuelle 
de la colonne 


Hauteur supérieure à la hauteur actuelle 


Actualisation des hauteurs 


Fig. 9. - Organigramme "Hcol” 


Listing du sous-programme ’Hcol’. 


3890 1 XP OPONOENNNNNMNHNNNNNNNMNNNNENNNNNNNNINONNENNNRNNNNNX 
3900 SUB Hcol(I,Nd1,Nnelen,INTEGER Liais(#),Maxa(t)) 
3910 1 XL EPONNENNNNNHNHNNNNRPNRONNNNENNRHHNPEEHNRNRRRNRPONNENNX 


3926 ! 


3940 ! 
3950 OPTION RASE 1 
3960 ! 


3930 ! Appelé par Frog. Princ. 
i 


3970 ! Calcule les hauteurs des colonnes 


3960 ! de la matrice A. 
3990 ! 
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4000 Mini=9997 


4010 ! 

4026 FOR J=1 TO Nnelen 

4030 IF Liais(i,J(Mini THEN Mini=Liais(1,J) 
4040 NEXT J 

A0S0 ! 

4060 FOR J=1 TO Nnelen 

A070 Nhaut=(Liais(l,J)-Mini)kNdl+i 
4080 Li=(Liais(1,J)-1)KNd1+i 

4090 IF Maxa(Li))=Nhaut THEN 4130 
4100 FOR L=1 TO Ndl 

Aâ10 Maxa(Li+tl-{)=Nhaut+L-1 
4120 NEXT L 

A130 NEXT J 


4140 SUBEXIT 
AÂ50  SUBEND 


2.5 - CALCUL DES INDICES DES TERMES DIAGONAUX 


Une fois toutes les hauteurs de colonne déterminées, le sous- 
programme ‘Indiag’ calcule les adresses des termes de la matrice ‘A’ 
situés sur la diagonale. Il suffit d'effectuer le cumul des hauteurs de 
colonnes. 


On détermine ainsi le nombre total de coefficients utiles de la 
matrice du système linéaire ‘Ntot’. Comme on connaît le nombre 
de cas de charges, on peut vérifier si la taille du tableau ’A' est 
suffisante pour traiter le problème. On doit avoir : 


Dimension de ‘A’ }= Maxa(Nnt 1)—1 + Nchar X Nn 


Si cette condition ne se trouve pas vérifiée, deux essais peuvent 
être effectués : 


— Essayer de renuméroter les nœuds de la structure, afin de 
diminuer les hauteurs de colonnes. 


— Redimensionner le tableau ‘A’ au détriment des cas de char- 
ges, ou du tableau ‘R'’. Sinon, la structure est trop importante. 


L’organigramme du sous-programme ‘indiag' est donné fi- 
gure 10. 
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Maxe(l+i)=Mara(T)+M 


Ntot=Maya(Nn+i)-f 





M=Maxa() 







Hodification de Maxa(X) par cumul 
des hauteurs de colonne 





N=Maxa(T+1) 






MEN 






Ntet : nombre total de coefficients utiles 
de la matrice À 





Fig. 10. - Organigramme lindiag’ 


Listing de ‘Indiag' : 


4160 
4170 
4180 
4190 
4200 
4240 
4220 
4230 
4240 
4250 
4280 
4270 
4280 
4290 
4300 
4310 
1320 
4330 
4340 
4350 
1360 
4370 
1380 


RETITICES CT ET OCT TT ET TT TER CE TIR TE LT TS TT PT CT TTL TR LT TT LT ES 
SUB Indiag(Nn,Ntot,INTEGER Maxa(X)) 

ME LLEUILE CSS TU TER PTT ET TS TE TELE LT TES TT LOTS TT ET ET TT ETES 
| 


fppelé par Prog. Princ. 
! 

OPTION BASE i 

1 


! Calcul des adresses des termes diagonaux de la matrice À 
! 
M=Maxa(1) 
FOR I=i TO Nn 
N=Maxa(1l+1) 
Maxa(l+i)=Maxa(1)+M 
MEN 
NEXT 1 
l 


1! Ntot=Nonbre de coefficients de la demi matrice du système 
l 


Ntot=Maxa(Nn+i)-i 
; 


SUBEXIT 
SUBEND 
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2.6 - ASSEMBLAGE DE LA MATRICE DU SYSTÈME LINÉAIRE 


Un deuxième parcours des éléments permet d’effectuer les opé- 
rations suivantes : 

— Former (ou lire sur mémoire auxiliaire) la matrice de rigidité 
de l’élément. 

— Calculer les adresses des coefficients de cette matrice dans le 
tableau A’. 

— Cumuler ces coefficients dans le tableau 'A'. 


Deux tableaux auxiliaires sont dimensionnés dans ‘Assemb' : 


— Lm(*), tableau des adresses des degrés de liberté de l’élément 
dont la dimension est : Ndi X Nnelem 


— Rig(*), qui contient les coefficients de la demi-matrice de 
rigidité de l’élément et dont la dimension doit être : 


(Ndi X Nnelem)X (Ndi X Nnelem+ 1)/2 


Listing du sous-programme ’Assemb’ : 


4390  ! XX AONNNNNNNNNNNHNNNNNNNNNNNNNNNNNNNNNONEETOPOONNNN 
A400 SUK Assemb 
4440 KANNNNNMHONNNNNNERNNMNNNNNNNNNNNNNONNNNNNNNORIONNINNNE 
A420 ! 
4430 ! Appelé par Prog. Princ. 
4440 ! Appelle : Matrig, Adress, Inject 

| 


4450 

4460 OPTION KASE À 

4470 ! 

A480 ! Assemblage de la matrice globale 
4490 ! 


AS00 COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelen,Nnodul ,Nsect? 
4510 COM Titre$[80),Charget#(5)[80] 

AS20 COM X(200) ,Y(200) ,Z(200) ,E(10) ,Sect(30) 
4530 COM INTEGER Lia1s(400,4) ,Isup (200) 

4540 ! 

4550 CO Ndib,Nn,Ntot 

4560 COhi INTEGER Iappui(50),Maxa(60i) 

4570 COM A(6000) ,R(400 ,15) 

AS80 !- 

4590 DIM Ln(b),Rig(21) 

4600 ! 

4640 FOR I=1 TO Nbarres 

4620 | 
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4630 CALL Matrig(R(H),Rig(#) dl 1) 
[l 


1640 | 

4550 CALL Adress(Nnelen,Nd1,1,Lm(X),Lia1s(#)) 

A660 ! 

4676 CALL Inject{I,La(#),Rig(x) ,Hnelem,Nd1,ACK) ,Maxa(x)) 
1680 Û 

4690 NEXT I 


A700 SUBEXIT 
4710  SUKEND 


2.1 - CALCUL DES ADRESSES ET INJECTION 
DANS LA MATRICE GLOBALE 


L'emplacement des coefficients de la matrice de rigidité de la 
barre dans le tableau ‘A’ est défini à l’aide d’une table d’adresses Lm. 
Cela est illustré par l’exemple ci-dessous (fig. 11), relatif à une ossa- 
ture plane. Comme une barre (ou un élément) n’est pas forcément 
décrite par des nœuds numérotés suivant un ordre croissant, et que 























im=[sTalr7Tse 


Barre définie par 2-4 Barre définie par 4-2 


Fig. 11. - Emplacement des coefficients de la barre 2-4. 
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l’on ne considère que la demi-matrice supérieure, il convient d’en 
tenir compte en effectuant une éventuelle transposition. Dans 
l'exemple (fig. 11), on remarque qu'il faut distinguer deux cas, 
suivant que le numéro du nœud origine est plus grand ou plus petit 
que le numéro du nœud extrémité. 


Il est à noter que, dans le cas particulier des ossatures, cette 
distinction est inutile, compte tenu des symétries de la matrice de 
rigidité (une simplification pour un programme spécifique aux 
ossatures!). Il n’en serait pas de même dans le cas de la matrice de 
rigidité d’un élément de poutre travaillant en flexion. 


Les organigrammes des sous-programmes ‘Inject’ et ’Adress’ 
sont donnés figure 12. 






Nd12=NnelemkNal 


Transposer ci J1)]] 


li=La(l) 
Jj=Latk) 












Ta=Maxa(Jj)+J3-Ii 


B(Ja)=A(Ta)+Rig(L) Cumul dans la matrice 


glohale 


Fig. 12. - Organigramme ‘Inject’ 
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Indice=(Liais(l,l)-4)KNd] 


K=K+1 
LatK)=lndicet]-i 





La) contient les numéros des 
degrés de liberté de l'élément 





Fig. 12(bis). - Programme dress’ 


Listing des sous-programmes ’Adress’ et ‘Inject' : 


5310 
5320 
5330 
5340 
5350 
S360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 


14290900 000NHnpHNNOHNNHNNRNHIONNHINNMNNNONNNNIONEONENNENEENEONENNONONE 
SUE Adress(Nnelem,Nd1,1,Ln(%) ,INTEGER Liais(*)) 

La opHNNNHNNHAONHIHNNNHHNNNHNRNNHHNNNNONONNNNRONONEONEONONNNTENENEX 
I 


| Appelé par Assemb 
j 


OPTION BASE 4 
1 


! Calcule les correspondance des indices entre Rig et À 
l 
K=0 
FOR L={ TO Nnelen 
Indice=(Liais(I,L)-1)KNd1l 
FOR J=4 TO Ndl 
K=K+1 
Ln(K)=Indicet] 
NEXT J 
NEXT L 
SUBEXIT 
SUHEND 
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S510 
5520 
5530 
5540 
5550 
5560 
5570 
5560 
5590 
5600 
Séit 
5620 
5630 
5640 
5650 
5660 
Sé7t 
5680 
5690 
5700 
5710 
5780 
5730 
5740 
S750 
5760 
5770 
5780 
5790 
5800 


1 KPNNNNHNNNHIHNNNHNNNNNNENNNNNENNNNONENNNENNONNNNONEOTEEN 
SUK Inject(I,LmCXk) ,R1g(X) ,Nnelem,Nd1,ACK),INTEGER Maxa(X)) 

1 KOOMNNNNNNNNNNNNNNPNNNONNNNAONIOONNNNNNNNCIOOENENNNNNNONNONEOEINX 
l 


| Appelé par Assemb 
1 


OPTION BASE à 
| 


| Injection de Rig dans À 

1 

Ndl2=NnelenkNdl 

L=0 

FOR J=1 TO Nd12 

FOR K=J TO { STEP -{ 

LEL+{ 
Ti=La(K) 
Jj=Ln(J) 
IF li(=Jj THEN 5750 
l 


| Si Ii ) Jj 
! Transposer 
[ 
Hi=La(]) 
Jj=Lmik) 
la=Maxa(Jj)+]j-1i 
Atla)=A(Ia)+Rig(l) 
NEXT K 
NEXT J 
SUREXIT 
SUEEND 


2.8 - CHARGES NODALES DUES AUX DÉPLACEMENTS D'APPUIS 


Si le déplacement d’un des nœuds est connu, il en résulte 
des efforts dans la structure que nous nous proposons de déterminer. 


Considérons d’abord à titre d'exemple un système linéaire 3X3 : 


2 —1 -05|1X l 
— 1 2 —05||Y|[=10 
—0.5—0.5 1 Z 0 


et supposons Y = 0.1 connu. 
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On voit facilement que le système peut encore s’écrire de la 
manière suivante : 


2 0 —-05S|IX 1 — (—1.0)X 0.1 
0 1 0 Y|=10.1 
— 0.5 0 1 Z 0 — (—0.5)X 0.1 
Pour arriver à ce résultat, il a fallu effectuer les opérations 
suivantes : 


1. multiplier la deuxième colonne par 0.1 et soustraire du se- 
cond membre, 


2. placer la valeur connue 0.1 de Y dans le second membre, 


3. annuler la deuxième ligne et la deuxième colonne de la 
matrice, 


4. donner la valeur 1 au terme diagonal de la deuxième ligne. 
Le sous-programme 'Genererdepl’ effectue les opérations 1 et 
2, pour chaque cas de charge. Les coefficients des seconds membres 


sont placés également dans le tableau A’, à la suite des coefficients 
de la matrice. 


Ar 





L 
Le 
a 





Fig. 13. - Charges nodales provenant d'un déplacement vertical de l'appui 3. 
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Qui: pas de forces nodales 


Eoucle sur les degrés de 
liberté bloqués 


ee 
(, 


K=Iappui(l) 

















Numéro du degré de liberté bloqué 





Min=Maxa(K)+4 
Max=Maxa(K+i)-1 


Ki=Indice+K-1 Kl=Adresse second membre 





Ki=K1-i 
A(K1)=ACK1)-A(L)KD(K) 


Sue 


Nhaut=Haxa(N+1)-Maxa(N) 


Gui: eviter la multiplication 
car coefficient nul 


LEMaxa (N)+N-K 
Ki1=Indice+N-1 
A(K1)=A(K1)-A{LIKD(K) 





Déplacements connus dans 
le second membre 
Numéro du degré de liberté bloqué 


K=ïappui(l) 
A(Indice+K-1)=D(K) 


Fig. 14. - Organigramme ‘Genererdepl' 


Les déplacements des appuis sont rangés dans un tableau 
auxiliaire 'D', initialisé à O0 au début de la description d’un nouveau 
cas de chargement. Un indicateur ‘lc’ permet d’éviter le produit 
des colonnes correspondant à un degré de liberté bloqué si le charge- 
ment ne prévoit pas de déplacement d’appuis. Il est à noter que cette 
opération nécessite la recherche des coefficients de la matrice ’A' 
sous la diagonale, donc à droite de celle-ci (fig. 13). 


L’organigramme de "Genererdepl' est donné figure 14. 


Listing du sous-programme 'Genererdepl' : 


7510 
75e0 
753 
7541 
7550 
7560 
7570 
7580 
7590 
7600 
7610 
7620 
7630 
7640 
7650 
7660 
7670 
7680 
7690 
7700 
7710 
7780 
7730 
7740 
7750 
7760 
7770 
7780 
7790 
7800 
7410 
7820 
7830 
7840 


1 XONNNNNONNANNANAHNNNNNNNNMANNNNNNONNNNNNNNNNNNNNOTACONNNNNNNACNONNX 
SUR Genererdepl(Indice,Ic,Nn,Nd1b,DCk) ,ACK),INTEGER Iappui(#) ,Maxa(4)) 
Lt PONEHANNRTONNRNNNNANTEPOPEONANNNENNIENNNCPENRONNNNORPONNNNX 
Û 
| Appelé par Second 
! 
OPTION BASE 1 
1 
l Inscrit les déplacements d’appuis connus (ou nuls) dans 
| le second membre. Genère les forces nodales dues aux dépla 
! cements d’appuis. 
! 
IF Ic=0 THEN Second 
FOR 1=1 T0 Ndib 
K=Iappui(l) 
IF K=4 THEN 7810 
Min=Maxa(K)+1 
Hax=Maxa(K+1)-1 
Ki=IndicetkK-i 
[ 


! Colonne 

[ 

FOR L=Min TO Max 
K1=K1-1 
ACK1)=ACK1)-A(L)IXD(K ) 

NEXT E 

IF K=Nn THEN 7890 

l 

! Ligne 

l 

Min=K+i 

FOR N=Min TO Nn 
Nhaut=Maxa(N+1)-Maxa(N) 
IF Nhaut(N-K+i THEN 7880 
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7850 LEMaxa (N)+N-K 


7860 K1=Indice+tN-i 

7870 A(K1)=A(K1)-A(L)KD(K) 
7880 NEXT N 

7890 NEXT I 


7900 Second: ! 
7940 FOR 1={ TO Ndlb 


79eû K=Iappui(1) 
7930 A(Indice+K-1)=D(K) 
7940 NEXT I 


7950  SUBEXIT 
7960 SUBEND 


2.9 - PRISE EN COMPTE DES APPUIS DE LA STRUCTURE 


Ce sont les opérations 3 et 4 décrites dans le paragraphe précé- 
dent et réalisées par le sous-programme ‘Modif’, appelé une seule 
fois après que tous les cas de charge ont été décrits et tous les se- 
conds membres générés. 


L'exemple (fig. 15) montre l’état du tableau ‘A’ après cette 
phase dans le cas d’une petite ossature plane soumise à deux cas de 
charge. 








2 4 
A PEN 
B 
o 
o 
Lo! 
o 
o 
2 
o 
lo] 


Fig. 15. - Second membre pour déplacement d'appui et charge nodale. 
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L’organigramme du sous-programme 'Modif' est donné figure 16. 









A(L)=1) 





Dernière colonne 
{aucun terme à droite) 


Nhaut=Maya(l+i)-Maxa(Ll) 
Colonne tron basse : pas de 
pas de multiplication 

K1=Maxa(L)+L-T 


A(K1)=0 


K1=haxatT) 
A(K1)=4 








Fig. 16. - Organigramme ‘Modif’ 


Listing du sous-programme Modif : 


8250 
8260 
8270 
8280 
8290 
8300 


RTTITIIITIT III TITI TITI TITI III TI TITI IIIDIT III CIE 
SUB Modif (Nn,Ndlb ,ACX) ,INTEGER Maxa(4) , lappui (x) ) 
RTTTITIT IEC TA TITI TITI III TITI LILI LIT IIITIT III TEL DII II CETTE 


Appelé par Prog. Princ. 
| 
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8310 OPTION BASE 1 

8320 ! 

8330 ! Annuler les coefficients de la ligne et de la colonne pour le 
8340 ! degré de liberté bloqué et donner la valeur { au terne diagonal. 
8350 ! 

8360 FOR I=i TO Ndlb 


8370 J=lappui(1) 

8380 ! 

8390 ! Colonne 

8400 k 

8410 Min=Maxa(J) 

8420 Max=Maxa(J+1)-1 

8430 FOR L=Min TO Max 

8440 A(L)=0 

8450 NEXT L 

8460 IF J=Nn THEN 8560 

8470 ! 

8460 |! Ligne 

8490 ! 

8500 FOR L=J+i TO Nn 

8510 Nhaut=Maxa(L+i)-Maxa(L) 
8526 IF Nhaut{L-J+1 THEN 6550 
8530 K1=Maxa(L)+#L-J 

8540 A(K1)=0 


8550 NEXT L 
8560 K1=Maxa(J) 
8570 AtKL)=i 


8580 NEXT I 
8590 SUBEXIT 


8600 SUREND 


2.10 - RÉSOLUTION DU SYSTÈME LINÉAIRE 


Nous utiliserons la méthode de factorisation de GAUSS, modi- 
fiée en vue de son utilisation en calcul automatique. Nous suppose- 
rons que la matrice du système linéaire est définie positive : s’il n’en 
est pas ainsi, soit la structure est instable, soit elle ne correspond 
pas à une structure réelle, ce qui provient généralement d’une erreur 
de données. Rappelons que cette méthode comporte deux phases : 

— La première phase transforme la matrice du système linéaire 
de façon à obtenir une matrice triangulaire supérieure. C’est la 
’factorisation’ ou ‘triangularisation’, ou encore ‘élimination’. 


— La seconde phase est la phase de substitution, permettant de 
calculer les inconnues les unes après les autres. 
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Illustrons la phase triangularisation de l’algorithme de GAUSS 
sous sa forme habituelle dans le cas d’un système 5 X 5 où les coeffi- 
cients non écrits sont nuls. 

















Fig. 17. - Système initial. 


On divise la première ligne (et le premier coefficient du second 
membre) par le premier coefficient de la première ligne que l’on 
appelle ‘pivot’. Ce pivot ne peut être nul pour des raisons physiques. 
On multiplie par les coefficients non nuls de la première colonne 
et on soustrait des lignes suivantes, sans oublier d’effectuer cette 
opération sur le (ou les) second membre. Le système est devenu : 











Fig. 18. - Après élimination de la première inconnue. 


On repète ces opérations sur la sous-matrice entourée d’un 
trait pointillé, et ainsi de suite jusqu’à l’obtention d’une matrice 
triangulaire supérieure avec des 1 sur la diagonale. Les étapes succes- 
sives sont : 














Fig. 19. 
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Fig. 19. - Étapes de la triangularisation. 


Le calcul des inconnues se fait alors par substitution. 


Nous constatons que, lors du processus de triangularisation, 
la sous-matrice entourée reste toujours symétrique : les coefficients 
situés sur une colonne sous la diagonale se retrouvent sur'la ligne 
de même indice à droite de la diagonale. On peut donc effectuer 
toutes les opérations sur la demi-matrice supérieure. De plus, on 
remarque que les zéros situés sur une colonne au-dessus du dernier 
terme non nul n’interviennent pas dans le processus d’élimination 
et par conséquent il n’est pas utile de les considérer. Cela justifie 
la méthode de rangement adoptée pour la matrice du système li- 
néaire. 


Malgré le fait que le processus de triangularisation est effectué 
ligne par ligne, il est préférable de mener les calculs colonne par 
colonne. On utilise pleinement les avantages du rangement par 
ligne-de-ciel décrit au paragraphe 2.2. En outre, cette méthode 
présente l’avantage de pouvoir être facilement adaptable à la réso- 
lution de grands systèmes dont la matrice ‘A’ ne peut tenir en mé- 
moire centrale. On effectue alors une segmentation sur disque. 


Exposons d’abord l’algorithme employé dans le cas d’une 


matrice rangée sous forme habituelle, c’est-à-dire où les coefficients 
sont repérés par un indice ligne et un indice colonne (fig. 20). 


I est l’indice de la ligne du premier terme non nul de la colonne J, 
M est l'indice de la ligne du premier terme non nul de la colonne K. 
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Fig. 20. - Indices utilisés lors de la factorisation. 


Schématiquement, l’algorithme peut se résumer de la manière 
suivante : 


J = 1... Nn 
K =H1::7 
fr: MAX(LM) …. K-1 


A(K,J)=A(K,J)—A(L,J)*A(L,K) 


L J-1 
A(J,J)=A(,J)—A(L,D#A(L,D/A(L;L) 
A(L,J)=A(L,J)/A(L,L) 
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La seule différence concerne la diagonale où l’on conserve les 
pivots, toujours nécessaires pour effectuer une nouvelle élimination. 
Cette phase est réalisée indépendamment du second membre et ne 
dépend donc pas a priori des charges. Nous aurons l’occasion de 
revenir sur ce point. 


La figure 21 illustre cet algorithme appliqué au système linéaire 
(fig. 17). 





2-4 


7 






RE 





Fig. 21. - Élimination colonne par colonne. 


L’organigramme détaillé de cette phase est donné figure 22. 


La seconde phase commence par la factorisation du second 
membre, qui est traité comme une colonne de la matrice ’A', et se 
poursuit par le processus de sustitution. A la différence de la métho- 
de de calcul traditionnelle, cette substitution n’est pas réalisée dans 
l’ordre décroissant, mais en parcourant les colonnes. En fait, on ef- 
fectue les mêmes opérations que dans la méthode classique, mais 
dans un ordre différent, moins ‘évident’. 


L’algorithme de la phase substitution est donné figure 23. 
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Jj=Maxa(J) 


TjrMara(J#4)-1 
Nhaut=1l)-Jj-i 














Kÿ=13 
If=J-Nhaut 











Te=icti 
KarKj-1 
Kk=Haxa(K) 
Nhaut=Mayra(k+1)-Kk-1 






Luax=MIN(ic,Nhaut) 





L1=haxa(k) 
CAC )/A( T1) 
ACJj)=A(Jj)-CKA(L) 
A(L)=C 


Foucie sur les rolonnes 


Aucun coefficient au-deseus 
du pivot 


Un seul coefficient 
au-dessus du pivot 


Pas de produit nossihle nour les 
Coefficients des colonnes K et ]J 


Division de la colonne par les 
pivots correspondants et calcul 
du J ème pivot 


ED x structure instahle 
ST OP 


Fig. 22. - Organigramme ‘Soive' (Factorisation). 
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Roucie sûr les seconds membres 


Ta=la+in lai : Adresse du 1 er coetficient 


du second membre 


Ji=Maxa(]i 
Tj=Maxa(J+1)-1 


Division du secand membre 
par les pivots 


Ja=J+1a 


L)=#axa(j) 
A(Ja)=A(Ja)/A(I.) 


Cas trivial 
Substitution 


La substitution ne s'effectue 
pas ligne après ligne, mais 
colonne anrèe colonne 


Ja=T+]a 
Ka=Ja 


Ka=Ka-1 
A(Ka)=A(Ka)-AC(LI*A(Ta) 


Fig. 23. - Organigramme ‘Soive' (Substitution). 


Listing du sous-programme 'Solve’ : 


8610 
8620 
8630 
8640 
8650 
6660 
8670 
8680 
8690 
8700 
8710 
8720 
8730 
8740 
8750 
8760 
8770 
8780 
8790 
8e00 
8810 
e8ot 
8430 
884û 
8850 
886l 
8870 
888ù 
8890 
8500 
8910 
get 
8930 
8940 
8950 
8960 
8970 
8960 
8990 
5000 
9010 
9020 
9030 
9040 
9050 
9060 


LBEPPNNNANANNHENNNHNHNNNHNERHNANNNNNNNRNNRANRNNRNINARNNNNNNNNNRNIENNNNNX 
SUE Solve(Nn,Htot,Nchar ,ACK) ,INTEGER Maxa(k)) 
D HOOONONNOPNNNNNNNNNNNNNNONNNNNENNNNRNNNNNNNNNRNNNNNINNNNRINNNE 


Appelé par Prog. Princ. 


OPTION BASE 1 


1 
1 
| 
! 
1 
! 
1 
l 
l 
1 
[l 
| 


| 
l 
1 
l 
l 
1 
l 
1 
l 


Résolution d’un système linéaire à matrice définie positive 


rangée par ligne de ciel, avec 


Nn = Dimension du système 


Not = Nombre de coefficients de la matrice dans À 


plusieur 


5 seconds 
Méthode de G À U 5 S adaptée au calcul automatique. 


Nchar = Nombre de vecteurs seconds membres 


Exemple : Nn = 5 Ntot = fi 


Matrice 


Ai4) A(3) A(b) 
A(a) A(S) A(9) 


A(4) AÀA(8) 
A(7) (ii) 
Aï40) 


faxa = [1, 8, 4, 7, 40, 1€) 


Nchar 


= 9 
+ È 


Seconds membres 


AIÂE) 
A(13) 
A(14) 
A(15) 
(16) 


Ati?) 
A(13) 
A(19) 
A(&d) 
&(2i) 


| Facterisation de la matrice du systeme linéaire 


FOR J=4 TO Nn 


Jj=Maxa(J) 
1j=Maxa(J#1)-1 
Nhaut=13-J3-1 
IF Nhaut(O THEN 7250 
IF Nhaut=( THEN 9170 
Ic=( 
Kj=1j 
11=J-Nhaut 
FOR K=I1 TO J-1 
Iczlc+i 
K3=Kj-1 
Kk=zMaxa(K) 
Nhaut=Maxa(K+i)-Kk-i 
IF Nhaut(=0 THEN 7140 


nenbres. 
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9070 
9086 
9090 
9100 
9410 
9120 
9130 
9140 
9450 
9160 
9170 
9180 
9490 
9200 
9210 
9ee0 
ges 

9240 
9250 
9260 
9270 
9280 
9290 
9300 
9340 
9320 
734 
9340 
9350 
9360 
9370 
9380 
9390 
9400 
9410 
9420 
9430 
9440 
9450 
9460 
9470 
9480 
9490 
9500 
950 
9set 
953 
9540 
9550 
$560 
7570 
9580 


Lnax=MIN(Ic,Nhaut) 
FOR L={ TO Lmax 
ACKj)=ACKj)-ACKk+LIKAUK j+L) 
NEXT L 
NEXT K 
3 
\ Division de la colonne par les 
l pivots correspondants et calcul 
| du J ème pivot 
1 
K=J 
FOR L=Jj+1 TO 1) 
K=K-1 
Li=Maxa(K) 
C=A(L)/A(L1) 
A(Jj)=A(Jj)-CHA(L) 
A(L)=C 
NEXT L 
IF A(Jj))0 THEN 9290 
PRINT “ STRUCURE INSTABLE :" 
PRINT " ERREUR DE DONNEES." 
PAUSE 


NEXT J 
1 Koucle sur les seconds membres 


La=Ntot-Nn 
FOR 1=1 TO Nchar 


Ja=latNn 
; 
| Factorisation du second membre 
! 
FOR J=1 TO Nn 
‘Jj=Maxa(J) 
1j=Maxa(J+i)-i 
1F 1j-Jj-1(0 THEN 9490 
Ja=J+1a 
Ka=Ja 
FOR K=Jj+i TO 1j 
Ka=ka-i 
A(Ja)=A(Ja)-A(K)#A(Ka) 
NEXT K 
NEXT J 
( 
! Division du second membre 
} par les pivets. 
1 
FOR J=1.T0 Nn 
Ja=J+]a 
L1=Maxa(J) 
A(Ja)=A(Ja)/A(L1) 
NEXT J 


9590 
9600 
9610 
9620 
9630 
9b4b 
9650 
7660 
9670 
7680 
9690 
9700 
9740 
9720 
9730 
9740 
9754 
7760 


IF Nn=4 THEN 9740 
1 


! Substitution 
1 
FOR J=Nn TO 2 STEP -{ 
Jj=Maxa(J) 
1j=Maxa(J+1)-1 
IF 1j-Jj-1(0 THEN 9730 
Ja=J+1a 
Ka=Ja 
FOR L=Jj+1 TO I j 
Ka=kKa-i 
A(Ka)=A(Ka)-A(L)#A(Ta) 
NEXT L 
NEXT J 
NEXT I 
SUEXIT 
SUREND 
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LES OSSATURES PLANES 
ET SPATIALES 


3.1 - GÉNÉRALITÉS SUR LES OSSATURES 


Ce sont des structures formées d’éléments linéaires ne travaillant 
que sous l'effet de l’effort normal : traction ou compression. Ce 
mode de fonctionnement est souvent justifié si la section de la barre 
est faible devant sa longueur, ce qui lui confère une grande souplesse 
vis-à-vis des déformations de flexion. C’est le cas en particulier des 
charpentes métalliques. Pour une structure plane, chaque nœud 
(assemblage de plusieurs barres) possède deux degrés de liberté : 
les deux composantes w et » de son déplacement. Pour une structure 
spatiale, le nombre de degrés de liberté est de trois : u, v et w. (fig. 1 
et 2). 


3.2 - LE PROGRAMME PRINCIPAL 


Le listing du programme principal, correspondant à l’organi- 
gramme donné au paragraphe 2.1, est reproduit ci-dessous. Les 
tableaux sont dimensionnés pour des ossatures de 200 nœuds et 
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400 barres avec 50 degrés de liberté bloqués. Le nombre de modules 
différents est 5, et le nombre de sections différentes, 30. Cinq cas 
de charges sont possibles. 


Listing du programme principal : 


40 ! PRRERRRENEEAEREEUEENEUEEEEENEEEEEEENEEEEEAEEEENEEENE 
30 !* CALCUL DES OSSATURES PLANES OU SPATIALES * 
40  !à+ * 
59  !* PAR LA METHODE DES DEPLACEMENTS * 
CE * 
70 XX TT TTL XXE HKÉEHNE À 
80 OX OX *X  * X OX OX XX à ï * 
90 IX * CE” * * * * * 
100 OX * X  XKYHRX OO HELERR * * 
140 OX * * * YO * x * 
{20 x  k OX À CC CRE k * * 
130 OX  XKk MT TTC ONE: * * * 
140 115 * 
150 L HHHHOPNNEHHPHEHNANMHHNNNNNNNRNNNRTMNNNNNNRHNNNNNNREX 
460 ©! 

170 1 ELOC COMMUN 


{80 OPTION EASE 1 
{90  PRINTER IS 7,1 
200 


240 ! Définition des variables 

220 ! 

230 ! 4°/ VARIAELES DEFINISSANT LE PRORLEME 

240 ©! 

250 1 Nnoeuds........... Nombre de noeuds 

260 ! Nbarres........... Nombre de liaisons 

870  ! Nchar............. Nombre de cas de charges 

280 !'NdE. us. Nombre de degrés de liberté par noeud 
290 ! Ndi = 2 Ossature plane 

300 ! Ndi = 3 Ossature spatiale 

310 ! Nnelem............ Nombre de noeuds d’un élément 

320 ! Nnelen = 2 pour une barre (2 noeuds) 
430 1 Nmodul............ Nonbre de modules d'élasticité 

340  ! Ncect............. Nombre de sections différentes 

450 ! Titre$............ Titre du probleme 

360 ! Charges$.......... Titres des cas de charge 

370 ! X(I), Y(I), Z(1)..Coordonnées du noeud I 

380 ! Isup(Il) .......... Code appui : Nombre de 2 chiffres = ossature plane 
390 ! Nombre de 3 chiffres ossature spatiale 
400 ! = 600 Noeud libre 

410 ! = 100 Noeud bloqué suivant X 

4ëG ! = {0 Noeud bloqué suivant Ÿ 

430 ! = { Noeud bloqué suivant Z 

440 ! = Autres combinaisons : 144 110 104 41 
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450 
46 
470 
480 
490 
500 
510 
52 

531 
540 
55G 
560 
570 
580 
590 
600 
bit 
6eû 
630 
640 
650 
bë&l 
670 
680 
690 
700 
748 
7eû 
730 
740 


760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
LU 
920 
930 
JAÛ 
950 
960 


| Liais(I,4)........ Numéro du noeud origine de la barre 


1 Liais(I,8)........ Nunéro du noeud extrémité la barre 
LLiais(l,3):%55 Numéro de la section Sect 
{ Liais(T,4)........ Numéro du module d’élasticité E 


COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelem,Nmodul ,Nsect 
COM Titre$[801,Charge$(5) [601 

COM xX(200),Y(200),Z(200) ,E(40) ,Sect(30) 

COM INTEGER Liais(400,4) ,Isup{200) 

l 


2°/ VARIABLES AUXILIAIRES EN COMMUN 


1 Ndib.............. Nombre de degrés de liberté fixés (appuis) 
ENT retenir NoeudskNdl : Dimension du système linéaire 
DONTORE AE se Re Nombre de coefficients de la matrice À 
LIappui(I) ........ Numéros des degrés de liberté fixés 

1 Maxa(T)........... Adresses des termes diagonaux de la matrice 
PACE), Le net Tableau des coefficients de la matrice et des 
! seconds membres 

LERCE T4) sta Mémoire auxillaire concernant les barres 


| 4 à 4: EXS/L ,Cosinus directeurs, Cosx,Cosy,Cosz 
| Sa 9: Efforts équivalents à la température 
| 40 à 44 : Efforts équivalents au poids propre 
COM Ndib,Nn,Ntot 
COM INTEGER Iappui(50) ,Maxa(60i) 

CON 4(6000) ,R(400,14) 

i 

? Initisiisation 

[ 

MAT Liais=(0) ! 

MAT A=(0) ! 
MAT R=(0) ! VARIABLES EN COMMUN 

MAT Maxa=(0)  ! 
MAT Isup=(0) ! 

Nsect=i 
Nnodul=i 

Nnelem=2 

l 


3 Lecture des données générales et des coordonnées des noeuds 
1 


CALL Lect 
il 


' Calcul des indices des degrés de liberté bloqués 
l 


CALL Indapp(Nnoeuds ,Nd1 ,Ndlb,Isup(X),Iappui(X)) 
l 


| Reproduction des données relatives à la structure 
1 


CALL Inprin 
l 
FOR I=i TO Nbarres 


! 
1 Formation de la matrice de rigidite de la liaison I 
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97 

9a0 

994 

1000 
1010 
1020 
1036 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1320 
1130 
1140 
1150 
1160 
4170 
1180 
1190 
1200 
1210 
igat 
1230 
1240 
1250 
1260 
1276 
1280 
4290 
{300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 


l 

CALL Format(l) 

! 

! Calcul de la hauteur des colonnes correspondantes 
[ 


CALL Hcol(I,Ndl,Nnelen,Liais(*),Maxa(x)) 
l 


NEXT I 

! 

| Calcul des indices des termes diagonaux 
[l 


CALL Indiag(Nn,Ntot,Maxa(X)) 
1 


| Assemblage de la matrice globale 
l 


CALL Assenb 
il 


! Lecture des charges et création du ou des seconds membres 
1 


CALL Second 
1 


! Modification de la matrice en fonction des appuis 
1 


CALL Modif (Nn,Ndlb AG) ,Haxa(X) , appui (1) 
! 


| Résolution du système linéaire 
l 


CALL Solve(Nn,Ntot,Nchar ,ACK) ,Maxa(X)) 
| 
l Pour chaque cas de charge, impression des déplacements des noeuds, 
l calcul et impression des réactions d’appuis, calcul et impression 
| des efforts dans les barres. 
: 
Indice=Maxa(Nn+{)-Nn ! permet de récuperer les déplacements de chaque 
! cas de charge dans le tableau À. 
PRINT LIN(2) 
PRINT "RESULTATS" 
PRINT 
FOR 1=1 TO Nchar 
! 
PRINT 
PRINT "CAS DE CHARGE N° ";I 


PRINT Charge#(1) 
l 5 


Indice=Indice+Nn 
fl ; 


| Impression des déplacements 
! 


CALL Depl(Indice,Nd1 ,Nnoeuds,A(K)) 
| 


1480 | Calcul et impression des efforts et 
1490 | des réactions d'appuis. 
1500 | 

1540 CALL Effn(Indice,1) 

1520 ! 

1530 PRINT LIN(2) 

4540 NEXT I 

1550 ! 

4560 DISP "FIN PROBLEME" 

4570 PAUSE 

1580 GOTO 730 

4590 END 


3.3 - LECTURE DES DONNÉES 


Ainsi qu'il a été spécifié dans le paragraphe 0.3, la lecture des 
données est prévue de manière toute formelle, en mode conversa- 
tionnel et sans aucune possibilité de correction. Il n’est pas conce- 
vable de générer ainsi les données de structures importantes, les 
erreurs étant inévitables. On pourra par exemple créer un programme 
où les données sont introduites en DATA. Ce préambule est valable 
pour les autres programmes présentés dans l’ouvrage. 


La structure est définie géométriquement par un certain nombre 
de paramètres généraux. Il faut fournir : 


— Le titre du problème, 


— Le nombre de dimensions (structure plane ou spatiale), ce qui 
définit le nombre de degrés de liberté par nœuds ’Nd/', 


— Le nombre de nœuds, y compris les appuis, ’Nnœuds’, 
— Le nombre de barres ‘’Nbarres’, 
— Le nombre de cas de charges, ‘Nchar', 


— Le nombre de sections différentes utilisées dans la structure, 
‘’Nsect', 


— Le nombre de modules différents, 'Nmodul'. 


Ces données générales permettent la lecture des coordonnées 
des nœuds, X(I), Y(D), (Z(L)), ainsi que les codes appuis, Isup(1). 


Les connexions sont définies par la donnée pour chaque barre, 
du nœud origine, du nœud extrémité, du numéro de la section 
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utilisée et du numéro du module d’élasticité. Toutes ces données 
sont regroupées dans le tableau 'Liais’. 


Ces définitions sont illustrées sur l’exemple ci-dessous (fig. 24). 





Nnœuds = 9 Nbarres = 15 Nsect = 1 Nmodul = 1 


(es) 


SD OPROROR 


2 
3 
4 
5 
6 
7 
8 
9 





1] Set | 
E() 


il 


1 
3 
5 
7 
2 
4 
6 
1 
3 
3 
5 
5 
7 
7 
9 


Co 00 A OO BB PB D D 00 A BB © I L 
= et et et et dt et ed et et et dt et jet pt 
= né ed ed et et et et dt et et et et pt jet 





Fig. 24. - Exemple de définition d'une ossature plane. 
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Listing du sous-programme ’Lect' 


1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
4700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 


1840 
1850 


1860 


1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
AU 
2010 
2020 
2030 
2040 
2050 
2066 
2070 
2080 


D NNNNNNNANNNNNNNNNNANTNNNNNNNNNNNNNNNNNNENNNEANNNONANIONNONNX 
SUE Lect 

1 XeonnnnNnNNbHNHHANTHNNNONHNHNHNHNHHNNNNNNNNNNNHINNONEHINNNNERHNNNIONX 
1 


! Appelé par Prog. Princ. 
1 


OPTION BASE { 


| Lecture des données relatives à l’ossature 

\ 

COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelen,Nmodul ,Nsect 
CON Titre$[801 ,Charge$(5)1801 

COM X(200),Y(200),2Z(200) ,E(10) ,Sect(30) 

COM INTEGER Liais(400,4) ,Isup(200) 

[ 

COM Nd1b ,Nn,Ntot 

COM INTEGER Iappui(50) ,Maxa(b0{) 

CON A(6000) ,R(400,44) 

1 

INPUT "Titre" ,Titre$ 

INPUT "Nombre de dimensions (Plan=2, Spatial=3)",Ndl 
INPUT "Nombre de Noeuds" ,Nnoeuds 

INPUT "Nombre de barres" ,Nbarres 

INPUT "Nombre de cas de charge" ,Nchar 

INPUT "Nombre de Sections (1)",Nsect 

INPUT "Nombre de Modules E (1)",Nmodul 
Nn=Nnoeuds*Nd1l 

} 


! Lecture des coordonnées des noeuds 
[ 
FOR 1=1 TO Nnoeuds 
IF Nd1=3 THEN 1940 
INPUT MX, Y, ",XCI),Y(I) 
GOTO 1950 
INPUT "X, YŸ, 2, ",XCI),Y(I),Z(1) 
INPUT "Code Appui ? (4)=Noeud libre",Isup(I) 
NEXT I 
l 


| Définition des liaisons 
1 
FOR I=1 TO Nbarres 
DISP "Noeud origine, Noeud extrémité" 
INPUT Liais(1,4) ,Liais(1,2) 
Liais(1,3)=1 
Liais(1,4)=1 
DISP "N° Section (1) N° Module (4) FAIRE CONTINUE Si { et 4" 
INPUT Liais(l,3),Liais(l,4) 


NEXT I 
Ê 
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2090 ! Lecture des sections 


2100 ! 

2140 FOR 1=4 TO Nsect 

2180 DISP "Section N° ";l 
2130 INPUT Sect(l) 

2140 NEXT I 

2150 : 

2460 ! Lecture des modules 
2470 ! 

2480 FOR I=i TO Nnodul 

2190 DISP "Module N° ";I 
2200 INPUT E(I) 

2240 NEXT I 

2220 ! 


2230 SUBEXIT 
2240  SUBEND 


3.4 - REPRODUCTION DES DONNÉES 


Quelle que soit la manière dont les données ont été introduites, 
il est nécessaire de les reproduire afin de bien savoir quelle structure 
a été calculée. Ces données doivent permettre de reconstituer inté- 
gralement le problème. Dans le cas présent des ossatures, la seule 
instruction ‘logique’ porte sur la distinction entre structure plane et 
structure spatiale. 


Listing du sous-programme ‘Imprim' 


2700 1! KR NNNNMHENNHNNMNINMNPNPNNNNNHNNNMNMNNNNNEMNNNENNIHNNRENNNNNENNNRRINNRNNNNX 
2710 SUB Imprin 
2780 1 XkHNNNNNNNNNNNNNEMMNNNNNNNNNNNTONNNNNNNNNNNNONNNNNEIEONNX 
2730 ! 
2740 ! Appelé par Prog. Princ. 

l 


2750 ! 

2760 OPTION BASE 1 

8770 ! 

2780 ! Reproduction des données générales de la structure 
2790 ! 


2800 COM Nnoeuds ,Nbarres ,Nchar ,Ndi ,Nnelem,Nmodul ,Nsect 
2810 COM Titre$1801,Charge$(5)[801 

2820 COM X(200),Y(200),2(200) ,E(10) ,Sect(30) 

2830 COM INTEGER Liais(400,4) ,Isup(200) 

2840 ! 
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2350 
2860 
2670 
2880 
2a90 
e900 
2910 
2920 
2730 
2940 
2950 
2760 
2970 
2980 
2990 
3000 
3940 
3020 
3030 
3040 
4050 
3060 
3070 
3080 
3090 
3100 
3410 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
2 
3230 
3240 
3250 
3260 
4270 
3280 
3890 
3300 
3310 
3320 
3330 


COM Ndlb,Nn,Ntot 

COM INTEGER Lappui(S0) ,Maxa(604) 

COM A(60@0),R(400,14) 

} 

PRINT "TITRE : "ATitreg 

PRINT 

IF Nd1l=2 THEN PRINT " Üssature plane axe Ÿ vertical" 
IF Nd1=3 THEN PRINT “ Üssature spatiale axe Ÿ vertical" 
PRINT 

PRINT "Nombre de noeuds = ";Nnoeuds 

PRINT "Nombre de barres = ",Nbarres 

PRINT "Nombre de cas de charges = ";Nchar 

PRINT "Nombre de sections différentes = ",Nsect 
PRINT "Nombre de modules d'élasticité = ";Nmodul 


PRINT 

PRINT "COORDONNEES DES NOEUDS" 

PRINT 

IF Ndi=2 THEN PRINT " N° * Ÿ Code appui" 

IF Ndl=3 THEN PRINT " N° X ' 2 Code appui 
PRINT 

FOR I=1 TO Nnoeuds 


IF Nd1i=2 THEN PRINT USING 4070,1,X(1),Y{I),Isup(l) 
IMAGE DDD,3X,MDDD.DDD,5X,MDDD.DDD, 41% ,DDD 
IF Ndi=3 THEN PRINT USING 3090;1,X(1),Y(1),Z(1),Iaup(l) 
IMAGE DDD,3X,MDDD .DDD,5X,MDDD. DDD,SX,MDDD. DDD ,9X,DDD 
NEXT I 
PRINT 
PRINT "DEFINITION DES LIAISONS" 
PRINT 
PRINT " N° Origine Extrémité N° Sect. N° Module" 
PRINT 
FOR 1=£ TO Nbarres 


PRINT USING 3180,1,Liais(1,11,Liais(1,8),,Liais(1,3),Lia1is(1,4) 
IMAGE DDD,SX,DDD,6X,DDD,8Xx,DD,8X,DD 
NEXT I 
PRINT 
PRINT "SECTIONS DES RARRES" 
PRINT 
FOR 1=i TO Nsect 
PRINT "Section N° ",1," = ";Sectil) 
NEXT I 
PRINT 
PRINT "MODULES D'ELASTICITE" 
PRINT 
FOR 1={ TO Nnodul 
PRINT "Module N° ",1," = ";E(I) 
NEXT I 
SUBEXIT 
SUKEND 
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3.5- CALCUL DES PRINCIPALES CARACTÉRISTIQUES DES BARRES 


Dans les grands programmes généraux, c’est souvent au cours de 
la première boucle sur les éléments que l’on calcule les matrices de 
rigidité. Afin de ne pas grossir démesurément la mémoire ’R', nous 
avons choisi de ne calculer que certaines valeurs : la rigidité de la 
barre ES/L et ses cosinus directeurs. Ces valeurs serviront à consti- 
tuer sa matrice de rigidité (sous-programme ’Matrig') et à calculer 
les efforts normaux à partir des déplacements. D’autres choix peu- 
vent être retenus en fonction des priorités accordées par l’utilisa- 
teur aux critères temps d’exécution et encombrement mémoire. 
Nous avons opté pour un compromis entre les deux solutions. La 
définition des cosinus directeurs est précisée (fig. 25). 





Cx = Lx/L Cx = Lx/L 
Cy = Ly/L Cy = Ly/L 
Cz = Lz/L 


Fig. 25. - Définition des cosinus directeurs. 


L’organigramme du sous-programme ‘Format’ est donné figure 
26. 
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liLiais(1,4) 
JisLiais(1,2) 
Kk=Liais(1,3) 
Lisliais(1,4) 









Cos(4)=X(Jj)-X(11) 
Cos(2)=Y(Jj)-Y(1i) 
Long=Cos(1)"2+Cos(2)*2 











Cos(3)=2(Jj)-Z(1i) 
Long=Long+tCos(3)*2 


STOP 


Long=S@R (Long) 
R(1,1)=E(L1)XGect(Kk)/Long 


RCI,J+4)=Cos(J)/Long 


Fig. 26. - Organigramme du sous-programme "Format. 


Listing du sous-programme Format’ 


3340 1 KP ENNNMNHNNMHHHNMNNNNHNNHHNNNMHHNNNNNNNNNNNHNNNNRNNNNNNNNNNEERONNNNNNRNNNNNRX 
3350 SUK Format(1) 

3360 1! HHOPONHNHHNNHHHHMHNNNNHNNEHNMHENNNNHNREENNNENNNNNHRRHNNRHRRNNX 
3370 ! 

3380 ! Appelé par Prog. Princ. 

3390 ! 

3400 OPTION BASE 

3410 ! 

3420 ! Calcul de la rigidité et de l'inclinaison de la barre 

3434 ! 

3440 COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelem,Nmodul ,Nsect 

3450 COM Titre$[80],Charge$(5)(80] 


3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3596 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3661 
3670 
3680 
3690 
3700 
3710 

780 
3730 
3740 
3750 
3760 
3770 


3.6 - MATRICE DE RIGIDITÉ D'UNE BARRE DANS LE PLAN 


On distingue la matrice de rigidité dans le système d’axes 
local, et la matrice de rigidité dans le système d’axes global. Dans le 
système d’axes local, c’est une relation qui relie les efforts et les 
déplacements suivant l’axe de la barre, et dont les coefficients ne 
dépendent que de la section, de la longueur et du module d’élas- 


COM X(200),*(200),2(200) ,E(10) ,Sect(30) 
COM INTEGER Liais(400,4) ,Isup(200) 

l 

CO“ Ndlb,Nn,Ntot 

COM INTEGER Iappui(50) ,Maxa(60i) 

COM A(6000) ,R(400,14) 

1 


DIM Cos(3) 

[ 

lisliais(1,i) 
JjLiais(1,2) 
Kk=Liais(1,3) 
LizLisis(1,4) 

| 

Cos(1)=X(Jj)-X(1i) 
Cos(2)=Y(Jj)-Y(11) 
Long=Cos(1)"2+Cos(2)*2 

IF Ndl=2 THEN 3660 
Cos(3)=2(Jj)-Z(1i) 
Long=Long+Cos(3)*2 
Long=S@R (Long) 

IF Long()0 THEN 3710 

PRINT "LONGUEUR NULLE POUR LA BARRE ";l 
STOP 

l 
R(I,1)=E(LL)XSect(Kk)/Long 
} 


FOR J=4 TO Ndl 
R(1,J+{)=Cos(J)/Long 

NEXT J 

SUREXIT 

SUBEND 


ticité du matériau. 


Afin de pouvoir écrire l’équilibre de chaque nœud de la struc- 
ture, il est nécessaire d’exprimer ces matrices de rigidité dans un 
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système d’axes unique, le système global. On utilise les relations 
entre les grandeurs mesurées suivant l’axe de la barre et dans le 
système global, exprimées figure 28 ci-dessous sous forme matri- 
cielle où l’on a posé : 


Cx=(Xj— XL ; Cy=(Yj—Yi/L 


ôi ôj 





Fig. 27. - Matrice de rigidité dans le système local. 


Il est important de noter que Cx et Cy sont des grandeurs 
algébriques qui dépendent de l’orientation de la barre : à est le 
nœud origine et j, le nœud extrémité. 





Fig. 28. - Changement de repère pour les structures planes. 
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On obtient ainsi la matrice de rigidité d’une barre pour une 
ossature plane. 





Fig. 29. - Matrice de rigidité d’une barre (ossature plane). 


3.7 - MATRICE DE RIGIDITÉ D'UNE BARRE DANS L'ESPACE 


Pour une barre dans l’espace, les relations entre les grandeurs 
mesurées dans le système local et celles mesurées dans le système 
global s’expriment par les relations (fig. 30). 


La matrice de rigidité de la barre dans l’espace est alors (fig. 31) : 


Ces matrices sont calculées par le sous-programme 'Matrig', que 
la structure soit plane ou spatiale, à partir des valeurs déterminées 
par le sous-programme ‘Format’ et rangées dans le tableau "R'. L’or- 
ganigramme du sous-programme ‘Matrig’ est donné figure 32. 
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Fig. 30. - Changement de repère pour une ossature spatiale. 






œ ce 
GC 


Fig. 31. - Matrice de rigidité d'une barre (ossature spatiale) 








ES 
L 
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Cste=R(1,1) 
Rig({)=R(I,2)*2 
Rig(2)=R(I,3)"2 

Rig(3)=R(I,3)KR(1,2) 












Jmax=21 
Rig(4)=R(I,4)"2 
Rig(S)=R(I,4)XR(I,3) 
Rig(6)=R(I,4)XR(1,2) 
Rig(7)=Rig(1) 
Rig(8)=-Rig(b) 
Rig(9)=-Rig(3) 
Rig(10)=-Rig(4i) 
Rig(ii)=Rig(e) 
Rig(i2)=Rig(3) 
Kig(13)=-Rig(5) 
Rig({4)=-Rig(2) 
Rig(15)=-Rig(3) 
Rig({6)=Rig(4) 
Rig(17)=R1g(5) 
Rig(i8)=R1g(6) 
Rig(19)=-Rig(4) 
Rig(20)=-Rig(S) 
Rig{21)=-Rig(é) 





















Jnax=10 
Rig(4)=Rig(i) 
Rig(5)=-Rig(3) 
Rig(b)=-Rig(i) 
Rig(7)=Rigte) 
Rig(8)=Riq(3) 
Rig(9)=-Rig(a) 
Rig(10)=-Rig(3) 



























Rig(J)=kig(J)xCste 





Fig. 32. - Organigramme ‘Matrig' 


Listing de 'Matrig' : 


4610! HNNNNNONNNNNNNNNNNNNNNONNNONNNNNNENOUNNNNNNENNNNNNNX 
4620 SUE Matrig(R(K),Rig(X) ,Ndl,1) 
4630! HPHPNNNNNNPNNNNONNNNNNNNNNENNNNNPNNNNTUNNNNNNNNNNNNNNNX 


4640 ! 


4660 ! 
4670 OPTION HASE { 
4680 ! 


4650 ! Appelé par Assemb 
! 


4690 ! Calcul de la matrice de rigidité de la barre sous la forme 


4700 ! 
4710 ©! Rig(i) Rig(3) 
4720 ! Rig(2) 
4730 ! 
4746 ! 


4750 Cste=R(I,1) 


4760 Rig(i)=R(I,2)*2 
4770 Rig(2)=R(I,3)"E 


Rig(b) 
Rig(5) 
Rig(4) . 


4780 Rig{3)=R(I,2)KR(I,3) 


4790 IF Ndl=3 THEN 4890 


4800 Jmax=i0 

AR1Q Rig(4)=Rig(i) 
4826 Rig(5)=-Rig(3) 
4930 Rig(b)=-Rig(1) 
4840 KRig(7)=Rig(e) 
4850 Rig(8)=Rig(3) 
4860 Rig(7)=-Rig(2) 


470 Fig(10)=-Rig(3) 


4466 GOTO 5090 
4890 Jmax=2 


4900 Rig(4)=R(I,4)*2 


4510 Rig(S)=R(I,4)#R(I,3) 
4926 R1g(6)=R(I,4)KR(I,? 


4930 Rig(7)=Rig({) 
4940 Rig(&)=-Rig(é) 
4950 Rig(9)=-Rig(3) 


4960 Rig(i0)=-Rig(i) 


4970 Rig(ii)=Rig(e) 
4980 Rig(i2)=Rig(3) 


4590 Rig({3)=-Rig(5) 
5000 Rigti4)=-Rig(2) 
S040 Rig(15)=-Rig(3) 


5928 Rigtib)=Rigi4) 
5030 Rig(17)=Rig(5) 
5040 Rig({8)=Rigtb) 


5050 Rig(i9)=-Rig(4) 

5060 Rig(20)=-Rig(5) 

5070 Rig(21)=-Rigtb) 
l 


Sat 


5090 FOR J=4 TO Jnax 


5100 Rig(J)=Rig(J)XCste 


S440 NEXT J 
5420 SUBEXIT 
5130 SUBEND 
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3.8 - DÉFINITION DES CHARGES 


Le second groupe de données, relatives aux chargements de la 
structure, est défini ici dans le sous-programme ’Second’. 


Quatre types de charges sont prévus : 


1. Charges arrivant directement sur les nœuds. C’est le plus 
simple. Il suffit de placer directement les deux (ou trois) compo- 
santes au bon endroit dans le second membre du système linéaire. 
La variable ‘Indice’ donne l’adresse du premier coefficient du second 
membre en cours de chargement (les seconds membres sont placés 
à la suite des coefficients de la matrice dans le même tableau ’A'). 


2. Charge répartie provenant du poids propre de la barre. L’ac- 
tion de la pesanteur s'exerce suivant l’axe Y pour les structures 
planes comme pour les structures spatiales. Ce cas de charge est 
généré automatiquement à partir de la donnée de la densité du 
matériau. Cette densité est supposée être unique pour toute la 
structure et les charges nodales sont calculées à partir de la section 
des barres et des longueurs de celles-ci. On est obligé de conserver 
le poids de la barre en mémoire pour le calcul des efforts normaux 
dans les barres, ainsi que l’illustre la décomposition suivante : 





Fig. 33. - Décomposition de la charge provenant du poids propre. 


Les chargements 2 et 3 sont auto-équilibrés et, par conséquent, 
sont sans influence sur les déplacements des nœuds et les réactions 
d’appuis. Par contre, il faut tenir compte de la variation d’effort 
normal dans la barre due au chargement 2 lors du calcul des efforts 
normaux à partir des déplacements des nœuds. 


56 


3. Température. L'action de la température est assimilée à 
un effort axial provoquant le même allongement de la barre (fig. 
34). 


ESaAt.Cy 


F = ESaAt 






77_ÿj ESaAt.Cx 


7 — ESaAt.Cy 


Fig. 34. - Forces nodales x 





provenant de l'effet thermique. 


Indice=Maxa(Nn+{)-Nn 





Ic=0 
MAT D=(0) Renise à 0 du tableau 
des déplacements imposés 
Indice=Indice+Nn 


| 











CALL Genererdepl 


Fig. 35. - Organigramme du sous-programme ‘Second’ 
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Inoeud, 
F(4),F(2)? 


Inoeud, 
F(4),F(2),F(3)? 


li=Indice+Ndlk(Inoeud-1)-1 












ACTi+#K)=A(Li+K)+F(K) 


















Densité du matériau 





Tizliais(l,1) 
Jj=Liais(I,2) 
Kk=Liais(1,3) 






Long=(X(Jj)-K(1i))" 24 
CU HD-Yi))*2 






Long=Long+(Z(Jj)-Z(11))*2 
Long=SQR (Long) 
Po=Longksect(Kk)#Dens/2 


li=Indice+(1i-1)KNd1+i 
Jj=lndice+(Jj-1)4Nd1+1 






Longueur de la barre 









Adresse degré de liberté 
correspondant à Ÿ 











Cumul dans le second 
membre 


ACII)=A(Ti)-Po 
ACTjI=AUIj)-Po 


R(L,1+ 9)=Po 





Sauvegarde Fo 
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Tenperature 
ÉCRÉPRREREEEEEPRREEES 





Ni,Ne, Alpha, Tenper? 


CALL Generer tenp 












F(4),F(2),F(3)7 














Rep$='Oui’ 
In=(N-L)#Nd1 





Jj={,Ndl 





Autre noeud ? 





Comme dans la réalité il s’agit d’une déformation et non d’un 
état de contrainte, il faut sauvegarder ces efforts normaux équiva- 
lents de façon à en tenir compte au moment du calcul des efforts 
réels dans les barres. 


4. Les déplacements d’appuis. On introduit le numéro du nœud 
(qui doit être obligatoirement un appui) et on calcule les charges 
nodales équivalentes conformément au sous-programme 'Genererdepl' 
exposé au paragraphe 2.8. 


L’organigramme du sous-programme ‘Second’ est reproduit 
figure 35. 


Listing du sous-programme ‘Second’ : 


SbA 
5650 
5660 
Sé7i 
5680 
5690 
5700 
S710 
S7e0 
5730 
5740 
5750 
5760 
5770 
5780 
5790 
5800 
5610 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
5490 
5900 
5910 
5920 
5930 
5940 
595 
59360 
5970 


1 OOONOONONNNNNNNNENENNNNNNNNNNENONNEEREIERIEEENINNCOONOCNNNNNENERX 
SUE Second 

1 BONHOMME 
Î 

! Appelé par Prog. Princ. 

! Appelle : Generer temp, Genererdepl 

! 

OPTION BASE 

! 

! Introduction du chargement, reproduction des charges 
let constituiion du (ou des) seconds membres du système 
| 

COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelen,Nmodul ,Nsect 

CON Titre$[801,Charges$(5)1801 

COM X(200),Y(200),2(200) ,E(40) ,Sect{30) 

CON INTEGER Liais(400,4) ,Isup(&00) 

| 

COM Nd1b,Nn,Ntot 

COM INTEGER Lappui(50) ,Maxa(604) 

COM A(6000) ,R(400,14) 

j 


DIM F(3),Coef(3),D(600) 
1 
Indice=Maxa(Nn+i)-Nn 
FOR I1=1 TO Nchar 
Ic=0 ! Paramètre indiquant les déplacements d'appuis 
MAT D=(0) 
Indice=lndice+Nn ! Indice du fer Coef. du vecteur second membre 
PRINT LIN(3) | 
PRINT "CAS DE CHARGE N° ";I 
PRINT " 
PRINT 
INPUT " Titre du cas de charge ?",Charges$(I) 
PRINT Charges#{I) 
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60 


5980 
5990 
6000 
6t{0 
6020 
6030 
6040 
6050 
6060 
6076 
680 
6090 
6100 
biit 
6i20 
6130 
6140 
6150 
bi60 
bi70 
6180 
6190 
6200 
62410 
beet 
6230 
beal 
6250 
6&bÙ 
be7( 
6280 
6290 
6300 
6310 
6320 
6330 
634( 
6350 
6360 
6370 
6380 
6390 
6400 
6410 
64aû 
6430 
6440 
6450 
6460 
6470 
6480 
6490 
6500 
6510 


C 


À: 


ro 


PRINT 

Index=5 

DISP “ NATURE DU CHARGEMENT ?" 
DISF " 4 : Charges aux noeuds" 
DISP " & : Poids propre" 

DISP " 3 : Température" 

DISP " 4 . Déplacement d'appui" 
INPUT " S : FIN" ,Index 

ON Index GOTO Ci,C2,C3,C4 ,Nexti 
| 


| Charges aux noeuds 
| 
PRINT "CHARGES AUX NOEUDS" 
PRINT " N° noeud Fx Fy Fi” 
DISP "La description du chargement est achevée par ia donnée” 
DISP "de la charge sur ie dernier noeud" 
IF Nd1=2 THEN 6190 
INPUT " N° noeud, Fx, Fy, F1" ,Inoeud,F(1),F(2),F(3) 
PRINT USING 6170; Inoeud,F(1),F(2),F(4) 
IMAGE X,DDD,2X,MDDD .DDD,2%X,MDDD.DDD,2X,MHDDD .DDD 
GOTO 6220 
INPUT ® N° noeud, Fx, Fy",Inoeud,F(4),F(2) 
PRINT USING b210;Inoeud,F(4),F(2) 
IMAGE X,DDD,2X,MDDD .DDD ,2X,MDDD .DDD 
li=Indice+Ndlk(Inceud-1)-1 
FOR K={ TO Ndl 
ACTi+K)=ACIi+K)+F(K) 
NEXT K 
IF Inceud=Nnoeuds THEN 5990 
GOTO 6120 ! Introduction de La charge sur le noeud suivant 
l 
| Génération poids propre 
[l 
PRINT "POIDS PROPRE (Toutes les barres sont prises en compte)" 
INPUT "Densité ?",Dens 
PRINT " Densité = “;Dens 
FOR L={ TO Nbarres 
Tisbiais(L,1) 
Jj=Liais(L,2) 
Kk=Liais(L,4) 
Long=(X(Jj)-X(11))"8+(Y(Jj)-Y(11))*E 
IF Nd1=2 THEN 6410 
Long=Long+(Z(Jj)-2Z(11))*2 
Long=S@R (Long) 
Po=LongkSect(Kk)kDens/2 
li=Indice+(Ti-{)XNd1+i 
Jj=indicet(Jj-1)4Nd1+i 
ACII)=ACTi)-P 0 
AJ jr=At I j)-Po. 
| 
l| Sauvegarde du poids propre pour le calcul 


! des efforts normaux et des réactions d'appuis. 
| 


R(L,149)=P6 


6520 NEXT L 

6530 GOTO 5996 

6540 ! 

6550 C3: ! Température 

6560 | 

6571 PRINT "TEMPERATURE" 

6580 | 

6570 DISP "N° première barre, N° dernière barre, Alpha et Température" 
6600 DISP “d’un groupe de barres (Cas de charge terminé obligatoirement" 
6610 DISP "avec le numéro de la dernière barre de la,structure" 
6bël INPUT "Ni, N2, Alpha, Temper ?",N1,N2,Alpha,Tenper 

6630 PRINT Ni,N2,Alpha,Temper 

6640 CALL Generertemp(N4,N2,Alpha,Temper ,Indice,1) 

6650 IF N2=Nbarres THEN 5990 

6660 GOTO 6620 | Lecture du groupe de barres suivant 

6670 | 

6680 C4: ! Déplacements d'appuis 

6690 ! 

4700 PRINT "DEPLACEMENTS IMPOSES DES APPUIS" 

6710 | 

6720 INPUT “ N° du Noeud ?",N 

6730 IF Isup(N)=0 THEN Erreur 

6740 Ic=ic+i 

6750 Nombre=lsup(n) 

6760 CALL Decode(Nombre ,Nd1,Coef(*)) 

6770 IF Ndl=2 THEN 6810 

6780 INPUT "DeltaX, DeltaY, DeltaZ ?",F(4),F(2),F(3) 

6790 PRINT "Noeud N° ",N," Dx = ";F(4)," Dy = ";F(2)," Dr = ";F(3) 
6800 GOTO 6830 

6ai0 INPUT "Deltax, Deltaÿ ?",F(4),F(2) 

6820 PRINT "Noeud N° ";N," Dx = ";F(14)," Dy = ";F(2) 

6830 Rep$="Qui" 

6840 In=(N-1)#Nd1 

6850 | 

6860 ! Seuls les déplacements imposés correspondants à des degrés de 
ser ! liberté bloqués sont pris en compte et placés dans le vecteur D 
6890 FOR Jj=1 TO Ndl 

6900 IF Coef(Jj)=0 THEN 6920 
"6910 D(ln+Jj)=F(Jj) 

6920 NEXT Jj 

6930 INPUT "Autre déplacement d'appui (Qui) ?",Rep$ 

6940 IF Rep$="Oui" THEN 6720 

6950 GOTO 5990 

6960 Erreur: DISP "Le noeud n’est pas un appui" 

6970 GOTO 6780 

6780 Nexti:! 

6996 ! 

7000 ! Modification du second membre en fonction des déplacements d’appuis 
7010 ! 

7020 CALL Genererdepl(Indice,Ic,Nn,Nd1b,DCK),ACX) ,Iappui(k),Maxa(#)) 
7030 ! 

7040 NEXT I 

7050 SUBEXIT 

7060 SUBEND 
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3.9 - CHARGES NODALES D'ORIGINE THERMIQUE 


Les charges d’origine thermique sont définies par la donnée d’un 
groupe de barres soumises à la même variation de température et 
possédant le même coefficient de dilatation. Le sous-programme 
‘Generertemp' génère alors les forces nodales équivalentes dans le 
système d’axes global à partir des caractéristiques des barres rangées 
dans le tableau ‘R'. 













Ji=liais(I,i) 
Jj=Liais(],2) 
Long=(X(Jj)-X(11))"2#(Y(Jj)-Y(T1))*2 








Long=Long+iZ(Jj)-Z(11))"2 


Long=SGR (Long) 
Eff=R(I,1)XAlphakTenper XLong 


R(I,Ichar+4)=Eff 


Longueur 
de la barre 





Sauvegarde pour 
calcul des efforts 


LisIndice+(Ii-1)XNd1-{ 
Mn=Indice+(Mn-1)XNd1-1 







Cste=EfFRR(I,J+1) Cunul dans le 
second mes#bre 
A(LI+J)=A(L1+I)-Cste 


AlMmtJ)=A(Mnt])+Cste 





Fig. 36. - Organigramme du sous-programme 'Generertemp’. 
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De plus, il faut conserver la valeur de cet effort normal pro- 
duisant la même déformation de la barre, ES aAf et servant à déter- 
miner les charges nodales. En effet, puisque cet effort n’existe pas, 
et que dans la réalité on a une déformation, lors du calcul des efforts 
normaux à partir des déplacements et des matrices de rigidité il fau- 
dra le soustraire. Comme le programme est dimensionné avec 5 de 
charges, ces efforts normaux équivalents sont rangés dans le tableau 
"R', colonnes 5 à 9. 


La description des charges d’origine thermique est achevée par 
la donnée d’un groupe de barres incluant la dernière barre de la 
structure, même si celle-ci n’est pas soumise à un effet thermique. 


L’organigramme du sous-programme ’Generertemp’ est donné 
figure 36. 


Listing du sous-programme ’Generertemp' : 


7070 1 KKHOONNNNNNNNNNNANNMNANNNNNNNNNNNNNNRNNNNNNNNNNNNNNEERHIERNRRNEEIRIRPERX 
7080 SUE Generer temp(Ni,N2 Alpha, Temper ,Indice,Ichar) 

7090 1 KMNMNNNNNNNNNNMHNNNNNNNNNNNNNNNNNONONENNNNNEEONENENNEEENNENNX 
7409 ! 

7110 ! Appelé par Second 


7480 ! 

7430 OPTION HASE i 

7440 ! 

7450 ! Genère les charges nodales dues à la température 
7160 ! 


7470 COM Nnoeuds,Nbarres ,Nchar ,Nd1 ,Nnelen,Nmodul ,Nsect 
7180 COM Titre$1663,Charge$(5)(801 

7490 COM X(200),Y(200),2(200),E(10) ,Sect(30) 

7260 CON INTEGER Liais(400,4) ,Isup(200) 

7210 ! 

7220 COM Ndlb,Nn,Ntot 

7230 COM INTEGER Iappui(50),Maxa(604) 


7240 COM A(6000) ,R(400,14) 

7250 FOR I=N1 TO N2 

7260 Jisliais(1,1) 

7270 Jj=Liais(1,2) 

7280 Long=(X(Jj)-X(11))"2+(Y(J3)-Y(11))*2 
7296 IF Nd1=3 THEN Long=Long+(Z(Jj)-2(1i))"2 
7300 Long=SQR (Long) 

7310 EFF=RIT,1)KAlphaXTemper {Long 

7320 R(I,Ichar+4)=Eff 

7330 ! 

7340 | Caicul des compsantes des forces dues à la 
7350 | température et cumul 

7360 ! 
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7370 
7380 
7390 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
7470 


Li=Indice+(1i-1)#XNd1-1 
Mu=Indice+(Jj-1)XNd1-1 
1 
FOR J=f TO Ndi 
Cste=EffxkR(I,J+1) 
AtL1+J)=A(L1+J)-Cste 
AiMmt])=A(MmtT)+Cste 
NEXT J 
NEXT I 
SUBEXIT 
SUBEND 


3.10 - IMPRESSION DES DÉPLACEMENTS 


Le sous-programme ‘Depl' imprime les déplacements après 
résolution du système linéaire. Il s’agit, en fait, d’extraire du tableau 
"A' les coefficients correspondant au second membre 
la procédure de résolution du système linéaire, les coefficients du 
second membre, initialement les charges nodales, sont maintenant 


les composantes des déplacements des nœuds. 


La seule particularité de cette séquence consiste dans la pré- 
sentation des résultats, suivant que la structure est plane ou spatiale. 
L’organigramme ne présentant aucun. intérêt, nous ne donnons que 


le listing reproduit ci-dessous : 


Listing du sous-programme ‘Depl' : 


9460 
9470 
74e 
9490 
9500 
9546 
9520 
9530 
9541 
9550 
9560 
9570 
9580 
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LLC E EEE CRE CTI TR ER TT TENTE TES TT TT TT TT TETE TT TITI TT TT TT TT 
SUE Depl(Indice,Nd1 ,Nnoeuds ,A(X)) 


Le TETE TT TT ETES LE TETITIT EEE ET TT LES EITTTT ETES ETTT 
1 


| Appelé par Prog. Princ. 


OPTION BASE 1 
| 


| Impression des déplacements 
| 


PRINT LIN(4) 
PRINT 'DEPLACEMENTS DES NŒUDS" 
PRINT 


: à la fin de 


9590 
3600 
961ü 
9629 
9630 
9640 
9650 
9660 
9670 
9680 
9690 
9700 
9740 
9720 
9730 
9740 
9750 
9760 
9770 
9780 
9790 
9800 


IF Ndl=2 THEN 9620 


PRINT " N° Noeud Dx Dy Dr" 
GOTO 9630 

PRINT " N° Noeud Dx Dy" 

PRINT 


FOR 1=i TO Nnoeuds 
lis=Iindice+(I-1)4Ndi 
IF Ndl=2 THEN 9760 
1 


| Ossature spatiale 

! 
PRINT USING 9740;1,A(11),A(I1+1) ,A(Ii+2) 

IMAGE 3X,DDD,SX,MD.DDDDDDE ,2X,MD.DDDDDDE ,2X, MD. DDDDDDE 
GOTO 9780 

| 


| Ossature plane 
l 


PRINT USING 9770,1,A(I1),A(li+4{) 

IMAGE 3X,DDD,5X,MD.DDDDDDE ,2X,MD.DDDDDDE 
NEXT I 
SUREXIT 
SUBEND 


3.11- EFFORTS ET RÉACTIONS D'APPUIS 


Les efforts normaux aux extrémités des barres sont calculés 
à partir des déplacements des nœuds. Les réactions d’appui sont 
obtenues ensuite par cumul des actions des barres reliées à un nœud 
appui. L’organigramme du sous-programme ‘Effn' est donné figure 37. 
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lisLiais(1,4) 
J j=Liais(1,2) 
K=Indice+(I1-1)4Nd1 
L=indice+(Jj-1)#Nd1 
Cos({)=R(1,2) 
Cos(2)=R(1,3) 
Dep(1)=A(L)-A(K) 
Dep(2)=A(L+1)-A(K+4) 















Cos(4)=R(1,4) 
Dep(3)=A(L+2)-A(K+2) 







Allong=0 


Allong=Allong+Dep(J)kCos(J) 


F(2)=DXR(T,4)-RC1,Ichar +4) 
F(4)=-F(2) 
Fo=R(I,Ichar+9)KCos(2) 

Effi=F(2)-Fo 
Effj=F(2)+F0 







Jj=Liais(1,J) 





I, 1,15,Effi,Eft) 





Reac(Jj,K)=Reac(J3,K)+R(I,Ichar+ 9) 


Calcul de l'allongement 
de la barre 


Prise en compte 
du poids propre 












Fig. 37. - Organigramme ‘Effn’. 


Listing du sous-programme 'Effn’ : 


78410 
9820 
9830 
984 
7850 
9860 
9870 
9880 
9390 
9900 
9910 
9920 
9930 
9940 
9950 
9960 
9970 
9980 
9990 
16000 
10010 
10020 
10030 
10040 
10050 
106060 
10070 
10080 
10090 
10100 
16410 
104ec 
10130 
10146 
10150 
10160 
10170 
40180 
10190 
10200 
102310 
4022 
10230 
10240 
10250 
10360 
10270 
10280 
10290 
10300 


(RLLTT Te CETTE TETE ET TETE TT EEE EE TEEN ET EE TE TETE EEE EE LT T EEE ETS 3 | 
SUE Effn{Indice,Ichar ) 

[M LLECE EN EEE TETE TUE ET TENTE TTT SET LUTTE ETES ET LITTLE EEE TE LES 
| 


| Appelé par Prog. Princ. 
| 


OPTION RASE 1 
| 
! Calcul et impression des efforts dans les barres. 
! Calcul et impression des réactions d’appuis. 
Il 
COM Nnoeuds ,Nbarres ,Nchar ,Ndi ,Nnelen,Nmodul ,Nsect 
CO Titre$[801,Charges$(5)180] 
COM X(200),Y(200),Z(200) ,E(10) ,Set(30) 
COM INTEGER Liais(400,4) ,Isup (200) 
i 
COM Ndlb,Nn,Ntot 
CON INTEGER Iappui(50) ,Maxa(604) 
COM A(6000) ,R(400,14) 
1 
DIM Coef(3),Cos(4) ,Reac(200,3),F(2) ,Dep(4) 
PRINT LIN(I) 
PRNT'EFFORTS DANS LES RARRES {+ Traction)" 
PRENT 
PRINT " N° Earre N° Orig N° Extr Eff. Norm. Or EfF. Norm. Ex." 
PRINT 
FOR I={ TO Nbarres 
Ji=Lia1is(1,1) 
Jj=Lia1s(l,2) 
K=Indice+(I1-{)xNdl 
L=ndice+(Jj-1)XNd1 
Cos(i)=R(1,2) 
Cos{2)=R(I,3) 
Dep({)=A(L)-ACK) 
Dep(2)=A(L+1)-A(K+1) 
IF Ndl=2 THEN 10190 
Cos(3)=R(I,4) 
Dep(3)=A(L+E)-A(K+2) 
Allong=0 
FOR J={ TO Ndl 
Allong=Allong+tDep(J)XCos(J) 
NEXT J 
l 


L Effort normal dans la barre 
l 


F(2)=AllongxR(1,1)-R(I,Ichar+4) ! Force à l'extrémité 
F(4)=-F(2) ! Force à l’origine 
! 


! Variation de l'effort normal due au poids propre 
} 
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16310 
10320 
10330 
10340 
10350 
10360 
10370 
10380 
10390 
10400 
10410 
10420 
10430 
10440 
10450 
10460 
10470 
10480 
16490 
19500 
10510 
10520 
10530 
10540 
10550 
10560 
10570 
10580 
106590 
10600 
16616 
10620 
10630 
10640 
10650 
10660 
10670 
10686 
10690 
10700 
10740 
10720 
10730 
10740 
10750 
10760 


Fo=R(I,Ichar+9)kCos(2) 

Effi=F(2)-Fo ! Effort normal à l’origine 
Effj=F(2)+Fo ! Effort normal à l'extrémité 
| 


! Calcul des Réactions d'appui. 
l 
FOR J=i TO Nnelem 
Jj=Liais(1,J) 
IF Isup(Jj)=0 THEN 10540 
Nombre=Isup(Jj) 
CALL Decode(Nombre,Nd1 ,Coef(*)) 
FOR K=1 TO Ndl 
IF Coef(K)=0 THEN 10500 
Reac(Jj,K)=Reac(Jj,K)+F(J)KCos(K) 
| 


! Addition des charges nodales provenant du poids 
l propre (Direction Ÿ (K=2) seule concernée) 
l 


IF K=2 THEN Reac(Jj,K)=Reac(Jj,K)+R(I,Ichar+9) 
NEXT K 
NEXT J 
! 
| Impression des efforts normaux 
| 
PRINT USING 10560;1,1i,Jj,Effi,Eff; 
IMAGE 4X,DDD,8X,DDD,8X,DDD,2(5X,MDDDDDD . DDD) 


NEXT I 


! Impression des Réactions d'appuis. 


1 


PRINT LIN(4) 


PRINT'REACTIONS D'APPUIS" 

PRINT 

IF Ndl=2 THEN 10670 

PRINT * N° Noeud Rx Ry R1" 
GOTO 10680 

PRINT " N° Noeud Rx Ry" 

PRINT 


FOR 1=1 TO Nnoeuds 


IF Isup(1)=0 THEN 40770 

IF Nd1=3 THEN 10750 

PRINT USING 10730,1,Reac(1,1) ,Reac(1,2) ,Reac(1,3) 
IMAGE 3X,DDD,3X,3(2X MD. DDDDDDE) 

GOTO 10770 

PRINT USING 10730;1,Reac(1,1) ,Reac(1,2) 

IMAGE 3X,DDD,3X,2(2%X,MD.DDDDDDE ) 


10770 NEXT I 


16780 


[NI 


CVTT 


JUBLAA: 


10790 SUREND 


3.12 - EXEMPLE DE CALCUL D'UNE OSSATURE PLANE 


L'exemple suivant concerne une poutre à treillis. La définition 
géométrique ainsi que les chargements sont précisés figure 38 et dans 
la reproduction des données. 





Fig. 38. - Test structure plane. 


Les données sont reproduites ci-dessous : 


TITRE : TEST STRUCTURE PLANE 


Üssature plane axe Ÿ vertical 


Nonbre de noeuds 9 

Nombre de barres {5 

Nombre de cas de charges = 4 
Nombre de sections différentes 


nou 
ra 


Nombre de modules d’élasticité 

COORDONNEES DES NOEUDS 

N° X Y Code appui 
{ 0.000 6.000 ii 
è 2.000 2.000 û 
3 4.000 6.000 ( 
4 6.000 2.000 (l 
5 8.000 5.000 Î 
b 10.000 2.000 ( 
7 42.000 9.000 û 
& 14.000 2.000 û 
9 16.000 0.000 {i 
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DEFINITION DES LIAISONS 


N° Origine Extrémité N° Sect. N° Module 


COQ SJ NN Un UT Ci Ci pe © BE TJ JU Le 
Doc » Bb FO ro © Oo & -©O —] U1 0 
pe pe Dé pe En pe HA ER HE pe RE pa 
OS 


A un pe EE pe 
Ne ro © © 0 Jo UT BR D rar 


SECTIONS DES KARRES 

Section N° = .00È 
MODULES D'ELASTICITE 

Module N° 4 = 21000000 


CAS DE CHARGE N° 1 


Cherge noeud 3 P=-20 


CHARGES AUX NOEUDS 
N° noeud Fx Fy 
3 ü.006 -20.000 
3 0.000 ÿ.0600 


F2 


CAS DE CHARGE N° 2 


Poids propre 


POIDS PROPRE (Toutes les barres sont prises en compte) 
Densité = 8 


CAS DE CHARGE N° 3 


Elévation de température 20° 


TEMPERATURE 
1 15 6000125 20 


CAS DE CHARGE N° 4 


Déplacement appui S de -0.0i 


DÉPLACEMENTS IMPOSES DES APPUIS 
Noeud N° 5 Dx = Dy = - Ûi 


RESULTATS 


CAS DE CHARGE N° 


Charge noeud 3 P=-20 


DEPLACEMENTS DES NŒUDS 


N° Noeud Dx Dy 

{ 0.000000E+00 0.000000E+00 
ë 1.517159E-03 -2.726994E-03 
3 5.730911E-04 -4 880891E-03 
4 -1.938013E-04 -2.630091E-03 
S 9.S523810E-04 0 .000000E+00 
b 1.938013E-04 6.2454{4E-04 
7 3.792898E-04  6.699948E-04 
8 3.876026E-04  5.246408E-04 
Le] 


0.000000E+00  0.600000E+00 


EFFORTS DANS LES KHARRES (+ Traction) 


N° Harre N° Orig. N° Extr. Eff. Norm. Or. Eff. Norn. 


i Î 3 6.017 
è 3 5 3.983 
3 5 7 6.017 
4 7 9 -3.983 
5 è 4 -17.965 
0 4 6 4.070 
7 6 8 3.035 
8 Î è -12.703 
9 3 2 12.703 
40 3 4 15.581 
it 5 4 -15.581 
ie S 6 -1.439 
43 7 6 1.439 
14 7 8 -1.439 
15 9 8 1.439 


6.017 
3.983 
-6.017 
-3.983 
-17.965 
4.070 
2.035 
-12.703 
12.703 
15.581 
-15.581 
-1.439 
1.439 
-1.439 
1.439 
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REACTIONS D'APPUIS 


N° Noeud RY Ry 
Î 2.96S086E+00 8.982543 +00 
S U.000000E+00 1.203494E+01 


ç _ 


ro 


.-6S086E+00 -1.0174S7E+00 


CAS DE CHARGE N° 2 


Poids propre 


DEPLACEMENTS DES NGŒUDS 


N° Noeud Dx Dy 
{ 0.000000E+00 0.000000E+00 
è 1.442723E-05 -2.708276E-05S 
2 1.277948E-06 -4.248235E-05 
4 -4.079705E-06 -2.79S981E-05 
5 -4.000000E-17  O.000000E+00 
6 4.07970SE-06 -2.795981E-05 
7 -1.277348E-06 -4.24823SE-05 
8 -1.112723E-0S -2.708276E-05 
$ 0.000000E+00 0.000000E+00 


EFFORTS DANS LES KHARRES (+ Traction) 


N° Karre N° Orig. N° Extr. Eff. Norm. Or. Eff. Norn. 


{ 1 3 .t13 .013 
2 3 S -.013 043 
3 S 7 -.013 -.013 
4 7 9 013 .D13 
S è 4 -.160 -.160 
6 4 6 .086 086 
7 6 8 -.160 -.160 
8 { è -.184 2,152 
9 3 a .042 .074 
40 3 4 .080 .4i2 
ii 5 4 -.267 -.235 
ie 5 b -.267 -.235 
13 7 6 .080 .442 
14 7 8 .042 .074 
15 9 8 -.184 -.152 


REACT 


N° Noeud 


1 
5 
9 


IONS D'APPUIS 


Rx 


4.050454E-01 
0.000000E+00 
-1.050454E-0i 


CAS DE CHARGE N° 3 


ZSIZSZSTSESSSSZZE= : 


Elévation de température 20° 


DEPLA 


N° Noeud 


Î 
è 
3 
4 
S 
6 
7 
8 
9 


CEMENTS 
Dx 


0.000000E+00 
-b.860346E-04 
.034914E-04 
-9.301728E-05 
.400000E-15 
-301788E-05 
.D34914E-04 
.860346E-04 
-000000E+00 


ro 


e © Fr) © un 


EFFORTS DANS 


N° Barre 


0 ur = cor 


N° Orig. N° Extr. 


NII UE B Fr JU re 


Ry 


4.730942E-01 
4.638562E-01 
1.730942E-0{ 


DES NŒUDS 
Dy 


0.000000E+00 
1.398254E-03 
5.930173E-04 
1.194763E-03 
0.000000E+00 
1.194763E-03 
5.930173E-04 
1.398254E-03 
0.000000E+00 


LES EARRES 
Eff. Norn. Or. 
3 -12.637 
5 -8.363 
7 -8.363 
9 -12.637 
4 -4.273 
6 -8.547 
8 -4.273 
2 -3.022 
2 3.022 
4 -3.022 
4 3.022 
6 3.022 
ë -3.022 
8 3.022 
8 -3.022 


{ + Tra 


Eff. 


ction) 


Norm. 


-12.637 


-8.363 
-8.363 
12.637 
-4.273 
-8.547 
-4.273 
3.022 

3.022 
-3.022 


.(?2 


ï t 

Cu Ci Ci Cid cJ 
eee 
ro F9 ro F9 M 
ro ro ro ln 
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REACTIONS D'APPUIS 
N° Noeud Rx Ry 
1 1.477332E+01  2.136659E+00 
5 0.000000E+00 -4.273319E+00 
9 -1.477332E+01  2.136659E+00 
CAS DE CHARGE N° 4 
Déplacement appui 5 de -0.01 
DÉPLACEMENTS DES NŒUDS 
N° Noeud Dx Dy 
{ 0.000000E+00 0.000000E+00 
è 8.034944E-03 -2.754364E-03 
3 -5.087284E-04 -6.017457E-03 
4 1.017457E-03 -8.263093E-03 
5 -1.700000E-14 -{.000000E-02 
ë -1.017457E-03 -8.263093E-03 
7 5.087284E-04 -6.017457E-03 
8 -2.034914E-03 -2.754364E-03 
9 0.000000E+00 0.000000E+00 
EFFORTS DANS LES BARRES 
N° Barre N° Orig. N° Extr. Eff. Norm. Or. 
{ { 3 -5.342 
2 3 5 5.342 
3 5 7 5.348 
4 7 9 -5.342 
5 è 4 -10.683 
6 4 6 -21.367 
7 b 8 -10.683 
8 { 2 -7.554 
9 3 è 7.554 
40 3 4 -7.554 
ii S 4 7.554 
{2 S b 7.554 
13 7 b -7.554 
{4 7 8 7.554 
45 9 8 -7.554 
REACTIONS D'APPUIS 
N° Noeud Rx Ry 
Î 1.068330E+01 S.341648E+00 
5 0.000000E+00 -1.068330E+04 
9 -1.068330E+01 5.341648E+00 


{ + Traction) 
Eff. Norm. 


-5.342 
5.342 
5.342 

-5.342 

10.683 
-21.367 
-10.683 

-7.554 
7.554 

-7.554 
7.554 
7.554 

-7.554 
7.554 

-7.554 


Ex. 














NZ 


Fig. 39. - Déformation sous poids propre 


















































Fig. 42. - Déformation sous déplacement de l'appui 5 


3.13- EXEMPLE DE CALCUL D'UNE OSSATURE SPATIALE 


La charpente tridimensionnelle présentée ici à titre d’illustra- 
tion du calcul d’une ossature spatiale présente un grand nombre de 


7S 


données. Celles-ci ne seront donc reproduites que partiellement. Les 
résultats sont visualisés à l’aide de dessins, ce qui met en évidence la 
nécessité des sous-programmes de représentation graphique. 


La définition géométrique de la structure et des cas de charges 
est présentée ci-dessous. 


4.00 4.00 4.00 , 4.00 








VWNVN 


Z 
Ÿ 
AT TAT 
———+ 
4.00 | 4.00 | 4.00 4.00 | 


Fig. 43. - Définition de la structure spatiale 
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Fig. 44. - Perspective de la structure spatiale 


TITRE : TEST STRUCTURE SPATIALE 


Gssature spatiale axe Ÿ vertical 
Nombre de noeuds ai 


Nombre de barres 128 
Nombre de cas de charges = à 


Nombre de sections différentes = { 

Nombre de modules d'élasticité = 1 

COORDONNEES DES NOEUDS 

N° X Y Z Code appui 
i 0.060 0.000 0.000 {ii 
2 4.000 0.000 0.000 ( 
3 8.000 0.000 0.000 û 
4 42.000 0.000 5.000 û 
5 16.000 0.000 0.000 iii 
b 2.000 2.000 2.000 { 
7 6.000 2.000 2.000 (i 
8 10.000 2.000 2.000 û 
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9 44.000 2.000 è.000 (l 
{10 0.000 0.000 4.000 î 
{ii 4.000 0.000 4.000 (l 
42 8.000 0.000 4.600 ( 
13 12.000 0.000 4.006 (r 
{4 46.000 0.000 4.000 (] 
45 2.000 2.000 6.000 (] 
16 6.000 2.000 6.000 ( 
17 10.000 2.000 6.000 (n 
18 14.000 2.000 6.000 (] 
19 0.000 6.000 8.000 û 
&0 4.000 0.000 8.000 (] 
ai 8.000 0.000 8.000 û 
22 42.000 0.000 8.000 0 
23 16.000 6.000 8.000 (1 
24 2.000 2.000 10.000 û 
85 .006 ë.000 10.000 0 
8b 10.000 2.000 10.000 (l 
27 14.000 2.600 10.000 (rl 
28 û.000 0.000 12.000 ( 
29 4.000 0.000 12.000 (i 
30 8.000 0.000 12.000 (l 
3i 12.000 6.000 42.000 û 
42 46.000 6.000 12.000 Û 
33 2.000 2.006 44.000 (! 
34 6.000 2.00 14.000 (l 
5 10.000 8.090 14.000 ( 
36 14.000 2.000 14.000 (l 
37 0.000 0.000 16.000 iii 
36 4.000 0.000 16.000 û 
35 8.000 0.006 16.00€ (l 
40 12.000 t.000 16.000 ( 
4i 16.000 0.000 16.000 {ii 
DEFINITION DES LIAISONS 
N° Origine Extrémité N° Sect. N° Module 
{ { 2 { { ë 82 23 
ë 2 3 { Î 13 28 89 
3 3 4 { { {4 29 30 
4 4 5 Î { 45 30 3i 
5 10 if { { 16 Ji 32 
6 ii {2 { i 17 37 38 
7 {2 13 î { 18 38 39 
8 13 14 Î îÎ 49 39 ac 
9 19 et { { ë 40 4i 
10 el ei { Î ei { 10 
11 21 2e Î { 22 10 {19 
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O0 


RE pe be pe px be pe 
MR pe EE pe De pe EE pa EX px 
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ee pe PE pe ee pe HE pa EE pe HE pe De pe EE pe EE pa te pe PE pa De pe De pe pe EE pe Et pa be pe be pa He pe He pe Er pa He pe HE pe EE pe De pre pe 


Domi 


pe 
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{20 
iei 
{22 
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M pee pe ee pe 2 
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SECTIONS DES HARRES 
Section N° { = .00e 
MODULES D'ELASTICITE 


Module N° 1 = 21000000 


CAS DE CHARGE N° 


Poids propre 


POIDS PROPRE (Toutes les barres sont prises en compte) 
Densité = 8 


CAS LE CHARGE N° 2 


Rotation globale autour des appuis { et 37 


DEPLACEMENTS IMPOSES DES APPUIS 


Noeud N° Of Dx = © Dy = .üi 
Noeud N° 37 Dx = 0 Dy = -.0i 
RESULTATS 
CAS DE CHARGE N° 
Poids propre 
REACTIONS D'AFPUIS 
N° Noeud Rx Ry R1 
4 8.4676S5E+00  1.910810E+00  2.467655E+00 
$ -2.467655E+00  1.910810E+00  2.467655E+00 
37 2.467655E+00  1.910810E+00 -2.467655E+00 
4i -2.46765SE+00  1.910810E+00 -2.467655E+00 


D1 
Dz 


LS 
ee COS we rune 


DONNE ro 


res 
[=] 





Fig. 45. - Déformation sous poids propre. 


CAS DE CHARGE N° 2 


Rotation globale autour des appuis { et 37 


DEPLACEMENTS 


Dx 


00000GE+00 
.438669E-04 
.22706SE-04 
.438669E-04 
.G00000E+00 
.349702E-04 
.374197E-04 
.374197E-04 
-349702E-04 
-324963E-04 


DES 


nTroEeS ou Jr 


NŒUDS 
Du Dz 
.000000€-02  0.000000E+00 
.793817E-03 -1.331963E-04 
.000000€-03  {.652000E-13 
.206183E-03  1.331963E-04 
.000060E+00  0.000000E+00 
.868865E-03  2.084970E-03 
.819411E-03  {.S4SS46E-03 
.680589E-03  9.544542E-(4 
.3413S4E-04  4.150298E-04 
.293817E-03 -2.438669E-64 


REACTIONS D'APPUIS 


N° Noeud Rx Ry R1 
{ 3.414137E+00 8.535343E-01 3.414137E+00 
5 3.414137E+00 -8.535343E-01 -3.414137E+00 
37 -3.414137E+00 -8.535343E-01  3.414137E+00 


ai -3.414137E+00 8.535343E-01 -3.414137E+00 





Fig. 46. - Déformation sous déplacements d'appuis. 
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LES POUTRES 
CONTINUES 


4.1 - GÉNÉRALITÉS 


Une poutre continue est une poutre droite reposant sur des 
appuis et soumise à des charges s’exerçant perpendiculairement à 
la fibre moyenne. Le problème que l’on se propose de résoudre 
consiste à déterminer les déplacements, les réactions d’appuis et les 
sollicitations dans la poutre sous l’action de charges extérieures 
(fig. 47). 


Pb 





Fig. 47. - Poutre continue chargée perpendiculairement à sa fibre moyenne. 


La théorie du calcul des structures élastiques débouche sur deux 
approches différentes : 


1. La méthode des forces, qui consiste à effectuer des coupures 
dans la poutre de manière à la rendre isostatique, puis à déterminer 
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les forces nécessaires pour rétablir la continuité. Cette méthode, 
dans le cas des poutres continues, se ramène à une série d’équations 
dites ‘équations des trois moments’. 


2. La méthode des déplacements, qui prend comme inconnues 
les déplacements d’un certain nombre de points de la poutre et dé- 
termine les valeurs de ces déplacements de façon à réaliser l’équilibre 
de l’ensemble. 


Comme nous avons axé l’ouvrage sur la méthode des déplace- 
ments, c’est celle que nous allons utiliser. De plus, cette méthode 
permet plus facilement de prendre en compte les diverses particula- 
rités de la structure, même si dans certains cas particuliers la méthode 
des forces est plus efficace (lignes d’influence). 


4.2 - DÉFINITION DES DONNÉES ET NOTATIONS 


Le fait qu’une poutre continue est constituée d’une succession 
d'éléments linéaires impose logiquement une numérotation naturelle 
des appuis et des travées. Nous supposerons que tous les tronçons 
sont de caractéristiques constantes. Cette restriction se justifie pour 
deux raisons : il est possible d’introduire des points intermédiaires 
de façon à suivre au mieux les variations des caractéristiques, ou bien 
il sera possible de créer des séquences de calcul spécifiques tenant 
compte de la variation de l’inertie et de la section. 


Le nombre de nœuds est au moins égal au nombre d’appuis de 
la poutré continue, mais il peut être supérieur. Le nombre de tronçons 
sera égal au nombre de nœuds moins un (fig. 48). 


1D2 © 3 O 4 566 ©  ! 
Re — 


Fig. 48. - Numérotation des nœuds et des barres. 


Afin de n’apporter aucune modification aux sous-programmes 
généraux, nous utiliserons le tableau des connexions Liais’. Chaque 
tronçon sera caractérisé par sa longueur et ses caractéristiques : iner- 
tie, module d’élasticité, éventuellement section réduite permettant 
de tenir compte des déformations dues à l’effort tranchant. 
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Le déplacement de chaque nœud est défini par sa flèche (dépla- 
cement vertical) et sa rotation. Il y a donc deux degrés de liberté par 
nœud. Si les deux composantes sont nulles, on est en présence d’un 
encastrement rigide. Si seule la flèche est nulle, on est en présence 
d’un appui simple. Le cas où seule la rotation est nulle est plus rare, 
mais sera envisagé. 


Afin de rendre le programme plus général, nous avons tenu 
compte de la possibilité d'introduire des articulations aux extrémités 
des barres : soit on libère le moment fléchissant, soit on libère 
l'effort tranchant. Cela est précisé par un code défini figure 49. 


Liais (1, 5) = 0 


RER 
md Liais (1, 5) = 1000 
Ÿo À Liais (1, 5) = 100 
N PR Liais (1, 5) = 10 

à —À Liais (I, 5) = 1 


Fig. 49. - Convention utilisée pour la définition des rotules. 


Le listing du programme principal, où les principales modifica- 
tions portent sur la définition des variables, est reproduit ci-dessous : 


Listing du programme ‘Poutr! : 


1 û RE DERREEEEEENEENEEEEEEEENNEEEENENEREENEEEENE 
30 !* CALCUL DES POUTRES CONTINUES * 
40 1+ X 
50  1+X PAR LA METHODE DES DEPLACEMENTS x 
60 1% * 
70 OX  XXXXXX XtX * X HO HO X 
80 1x * X OX XX * * * X * 
90 !X * X + X * * k * X * 
400 UK XKKKXX X * * * LS COLE SE: 
40 1x XX * X * * * X  X  * 
420 !X * X OX + * * * X * 
430 OX * xt Étteel * * XX 
440 1% X 
150 LE HHHHHNNNHHNHNHHMNNNEHUNRNHHNNNNNNNNENNNHNNNNNONNNONNNRNENX 
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160 
170 
180 
190 
200 
210 
220 


230 


e50 
260 
270 
280 
290 
300 
310 
320 


490 
500 
510 
sel 
530 
540 
550 


570 
580 
590 
600 
640 
620 
630 
640 
650 
660 


! 
| ELOC COMMUN 
OPTION EASE 1 
PRINTER IS 7,1 
; 
Définition des variables 


4°/ VARIARLES DEFINISSANT LE PROBLEME 


[ 
! 
[ 
; 
! Nnoeuds........... Nombre de noeuds 
| Nbarres........... Nombre de liaisons 
LNchañ sic Nombre de cas de charges 
PNA roues Nombre de degrés de liberté par noeud 
| Ndl = 2 Poutre continue 
[ 
| Nnelem............ Nombre de noeuds d’un élément 
! Nnelen = 2 pour une barre (2 noeuds) 
| Nnodul............ Nombre de modules d’élasticité 
l'Nsect. sis hs Nombre de sections différentes 
Titres 27m us Titre du probleme 
Charges$.......... Titres des cas de charge 
NET aa ares Abscisse du noeud I 
Isup(1) .......... Code appui : Nombre de 2 chiffres = poutre continue 


00 Noeud libre 


10 Noeud bloqué suivant Ÿ 
4 Noeud bloqué en rotation 
{1 Encastrement rigide 


CRUE EL 


Liais(1,4)........ Nunéro du noeud origine de la barre 
1 Liais(1,2)........ Numéro du noeud extrémité la barre 
Liais(1,3)........ Numéro des caractéristiques géométriques : 
Section réduite et Inertie 
Liais(I,4)........ Numéro des caractéristiques physiques : 
; Module de cisaillement G et module d'YOUNG € 
Liais(1,5)........ Code articulation aux éxtrémités des barres 


0000 Aucune articulation 

1000 Articulation de cisaillement à l’origine 
0100 Rotule à l’origine 

0040 Articulation de cisaillement à l'éxtrémité 
! 0001 Rotule à l’éxtrémité 

COM Nnoeuds ,Nbarres,Nchar ,Nd1 ,Nnelen;Nmodul ,Nsect 

COM Titre$(80),Charge$(5)[80] 

COM X(100),E(S,2),Sect(15,2) 

COM INTEGER Liais(200,5),Isup(100) 

! 


! 2°/ VARIABLES AUXILIAIRES EN COMMHUN 


MOU OH OU 


L'NdIb res Nombre de degrés de liberté fixés ‘appuis) 

LAND ds à crie enets NoeudskNdl : Dimension du système linéaire 

TENTE 25e ones Nombre de coefficients de la matrice À 

1! Lappuit(I) ........ Numéros des degrés de liberté fixés 

l Maxa(I)...........Adresses des termes diagonaux de la matrice 


670 
680 
690 
700 


720 
730 
740 
750 
760 
770 
780 
790 
800 
8it 
826 
830 
840 
850 
860 
870 
880 
870 
900 
gi0 
9e 
930 
940 
950 
9eû 
970 
980 
990 
1000 
1010 
4020 
1030 
1040 
1050 
1060 
4070 
1080 
1090 
1100 
4110 
{120 
4190 
1140 
4150 
1160 
1170 
1180 


MAGIE ame Tableau des coefficients de la matrice et des 
| seconds membres 
LRO SR8) 7 reed e Coefficients de la matrice de rigidité de la barre I 


| de { à 4, de S à 24, pour chaque cas de charge, les 
| efforts aux extrémités des barres. 

COM Ndib,Nn,Ntot 

COM INTEGER Iappui(50) ,Maxa(301) 

COM A(5150) ,R(200 ,24) 

l 


| Initialisation 
! 

MAT Maxa=(0) ! 
MAT A=(0) ! 
MAT Isup=(0) ! CAR VARIABLES EN COMMUN 
MAT Liais=(0) ! 
MAT R=(0) ! 
Nsect=i 
Nnodul=i 
Nnelen=2 

[ 


! Lecture des données générales et des coordonnées des noeuds 
l 


! CALL Lect 
1 


! Calcul des indices des degrés de liberté bloqués 
l 


CALL Indapp(Nnoeuwds,Nd1 ,Ndlb,Isup(x) ,Iappui(*)) 
! 


1 Reproduction des données relatives à La structure 
1 


CALL Imprin 
j 


FOR 1=i TO Nbarres 
1 


| Formation de la matrice de rigidite de la liaison I 
1 


CALL Format (1) 
1 


!{ Calcul de la hauteur des colonnes correspondantes 
[ 


CALL Hcol(I,Ndl,Nnelen,Liais(k) ,Maxa(x)) 
Il 


NEXT I 
1 


1 Calcul des indices des termes diagonaux 
1 


CALL Indiag(Nn,Ntot,Maxa(#)) 
l 


| Assemblage de la matrice globale 
} 


CALL Assemb 
! 
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4190 
4200 
12410 
{220 
4230 
1240 
4250 
1260 
1270 
1280 
4290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
i4el 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
4590 
1600 
4610 
1620 
1630 
1640 
1650 


| Lecture des charges et création du ou des seconds membres 
1 


CALL Second 
Ï 


1 Modification de la matrice en fonction des appuis 
1 


CALL Modif (Nn,Nd1b,ACK) ,Maxa (X) ,lappui (#)) 
! 


| Résolution du système linéaire 
l 


CALL Solve(Nn,Ntot,Nchar ,ACX) ,Maxa(X)) 

! 

L Pour chaque cas de charge, impression des déplacements des noeuds, 
! calcul et impression des efforts aux extrémités des barres et des 

| réactions d'appuis. 

! 

Indice=Maxa(Nn+i)-Nn ! permet de récuperer les déplacements de chaque 
! cas de charge dans le tableau À. 

PRINT LIN(2) 

PRINT “RESULTATS "* 


PRINT 


FOR I=i TO Nchar 
Ê 
PRINT 
PRINT "CAS DE CHARGE N° ",1 
PRINT "============2==2=2=2" 
PRINT Charge$(1) 
1 


Indice=Indice+Nn 
l 


! Impression des déplacements 
] 


CALL Depl(Indice,Nd1 ,Nnoeuds,A(X)) 
1 


! Calcul et impression des efforts et des 
| réactions d’appuis. 
1 


CALL Efforts(Indice, I) 
} 


PRINT LIN(2) 
NEXT I 
DISP “FIN PROBLEME. FAIRE CONTINUE POUR UN AUTRE PROBLEME * 
PAUSE 
GOTO 780 


END 
| 


4.3- LECTURE DES DONNÉES 


Les données définissant la structure sont particulièrement 
simples. Il convient de ne pas oublier d’initialiser la variable ‘Ndl' 
à 2 (deux degrés de liberté par nœuds). Le fait de ne pas donner une 
des caractéristiques, physique ou géométrique, concernant la flexion 
ou la distorsion conduit à considérer la barre comme sans déforma- 
tions vis-à-vis de cet effet. Nous n’avons pas envisagés le cas où un 
tronçon de barre se comporterait comme un solide parfaitement 
rigide. La définition des connexions est générée automatiquement. 
Le listing du sous-programme ’Lect’ est reproduit ci-dessous : 


Listing du sous-programme ’Lect' : 


1660 1 HONOR 
4670 SUB Lect 
1680 1! KHHHNNNHHMNMNMHMNHNNNHNNNNHNHRMNNNHMHNNNNNNNNNNNRNNNRNRNNNNNNNONONNNHNTOINNNX 
1690 ! 
4700 ! Appelé par Prog. Princ. 

] 


4740 ! 

4720 OPTION HASE 

1730. ! 

4740 ! Lecture des données relatives à l’ossature 
4750 ! 


1760 COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelen,Nnodul ,Nsect 
4770 COM Titre$[801,Charge$(5)[801 

1780 COM X(100),E(S,2) ,Sect(i5,2) 

1790 COM INTEGER Liais(200,9) ,Isup(100) 

1800 ! 

4840 COM Ndlb,Nn,Ntot 

4820 COM INTEGER Iappui(50) ,Maxa(301) 

1830 COM A(S150) ,R(200,24) 


1840 ! 
4850 INPUT "Titre",Titre$ 


1860 Ndl=2 

4870 INPUT "Nombre de barres" ,Nbarres 

1880 Nnoeuds=Nbarres+i 

1890 Nn=NnoeudskNdl 

4900 INPUT “Nombre de cas de charge" ,Nchar 

4910 INPUT "Nombre de Sections (1)”,Nsect 

4920 INPUT "Nombre de Modules E (1)",Nmodul 


1930 ! 

1940 ! Lecture des abscisses des noeuds 
4950 ! 

1960 FOR 1=i TO Nnoeuds 

1970 DISP " X(",1;,") 2° 


1980 INPUT X(1) 
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1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
ei10 
2120 
2430 
2440 
2154 
2460 
2170 
2180 
2190 
2e00 
2ei0 


eeet 


ro TO 
+ CG 
ee 


rJ 
u1 
CRE 


se 


ro ro PO ro FO ro FO ro F2 
+ > DO 
= 


S ci FO ro MO ro 


e 


INPUT "Code Appui ? (0)=Noeud libre",Isup(1) 
NEXT I 
j 


| Définition des liaisons 
| 
FOR 1=4 TO Nbarres 
Liais(1,4)=1 
Liais(1,2)=1#1 
Liais(I,3)={ 
Liais(1,4)=1 
DISP “HARRE N° ";l 
DISP "N° Section (4) N° Module (4) FAIRE CONTINUE Si i et°{" 
INPUT Liais(1,3),Liais(I,4) 
DISP "Code articulation ? (FAIRE CONTINUE SI PAS D'ARTICULATION) * 
INPUT Liais(1,5) 
NEXT I 


! Lecture des sections 

! 

FOR I=i TO Nsect 
DISP “Section réduite, Inertie N° ";l 
INPUT Sect(1,1),Sect(1,2) 

NEXT I 

1 


Lecture des modules 

! 

FOR 1=1 TO Nnodul 
DISP "Module G, module E N° ";l 
INPUT E(I,1),E(1,2) 

NEXT I 

| 

SUBEXIT 

SUKEND 


4.4 - REPRODUCTION DES DONNÉES 
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Le listing concernant le sous-programme de reproduction 
des données, particularisé essentiellement par des labels adéquats, 
est reproduit ci-après : 


Listing du sous-programme Imprim” : 


2770 
2780 
2790 
2800 
28410 
2820 
2830 
2840 
2850 
2860 
2876 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3040 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3400 
3140 
3420 
3430 
3140 
3450 
60 
3470 
3180 
3490 
3200 
35) 
3240 
3220 
3230 


1 epnnnNNNHHNNNHHHHHHHNNNHHNNNNHHHHNNHHNNNHNNRRMHNNNNHNIHNNNNNANNONHNNHENNRX 
SUB Imprin 

Le 0oNNNNNHHNNHHNHHIHNRHMHHNNNNHAHRNHHENNNNNNNRIEREMHONNAIHONNNNNERREONNNENX 
i 


| Appelé par Prog. Princ. 
! 

OPTION BASE 

! 


| Reproduction des données générales de la structure 
| 

CON Mnveuds ,Noarres ,Nchar ,Nül ,Nnelen,Nmodui ,Nsect 
COM Titre$(801,Charge$(5)180] 

COM X(400),E(5,2) ,Sect(15,2) 

COM INTEGER Liais(200,5),Isup(100) 

| 


COM Ndlb,Nn,Ntot 

COM INTEGER Iappui(50) ,Haxa(304) 
COM A(S150) ,R(200 ,24) 

l 


PRINT “TITRE : "ATitres 


PRINT 

PRINT " Poutre continue, axe Ÿ vertical" 
PRINT 

PRINT "Nombre de noeuds = “;Nnoeuds 
PRINT "Nombre de barres = “;Nbarres 


PRINT "Nombre de cas de charges = ";Nchar 
PRINT "Nombre de sections différentes = ";Nsect 
PRINT "Nombre de modules d’élasticité = ";Nmodul 


PRINT 

PRINT "ABSCISSES DES NOEUDS" 

PRINT 

PRINT * N° X Code appui" 
PRINT 


FOR 1={ TO Nnoeuds 
IF Nd1=2 THEN PRINT USING 3120;1,X(1),Isup{1) 
IMAGE DDD,3X,MDDD . DDD,14X,DDD 
NEXT I 
PRINT 
PRINT "DEFINITION DES LIAISONS" 
PRINT 
PRINT " N° Origine Extrémité N° Sect. N° Module Code articulation" 
PRINT 
FOR 1=1 TO Nbarres 
PRINT USING 3240;1,Liais(1,4),Liais(1,2),Liais(1,3),Liais(l,4),Liais(] 


IMAGE DDD,SX,DDD,6X,DDD,8X,DD,8X,DD,8X,DDDD 


NEXT I 
PRINT 
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3240 PRINT "CARACTERISTIQUES DES HARRES" 
3250 PRINT 
3260 FOR I=i TO Nsect 


3270 PRINT “Section réduite N° ";l;" = "isect(1,1) 
3280 PRINT “inertie N° ";l;" = “;Sect(I1,2) 


3290 PRINT 
3300 NEXT I 

3310 PRINT 

3320 PRINT "MODULES D'ELASTICITE G et E" 
3330 PRINT 

3340 FOR 1=1 TO Nnodul 

3350 PRINT “Module G N° "1," 
3360 PRINT “Module E N° ";1;" 
3370 PRINT 

3380 NEXT I 

3390 SUBEXIT 

3400  SUBEND 


“;EC1,1) 
*E(1,8) 


4.5 - MATRICE DE RIGIDITÉ D'UN ÉLÉMENT DE POUTRE EN FLEXION 


La relation entre les forces appliquées aux extrémités d’une 
poutre et les déplacements de celles-ci peut s’obtenir rapidement 
de la manière exposée ci-après. 

Soit un tronçon de poutre dont l’extrémité gauche est i et 
l’extrémité droite j. Tiet Mi l’effort tranchant et le couple appliqués 
en i, Tj et Mj l'effort tranchant et le couple appliqués en j. Les dé- 
placements des nœuds à et j sont respectivement Vi, Gi, Vj, et Oj. 
Ces grandeurs sont positives sur la figure 50. 


Ti Tj 


Nu 


Oj 
Vi 






Mi Mj 
Fig. 50. - Définition des forces et des déplacements. 


Les forces appliquées aux extrémités de la barre doivent satis- 
faire aux équations.d’équilibre : 
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Ti+Tj=0 
Mi+Mj+TiL=0 ou Mi+Mj-TiL=0 


jou encore, sous forme matricielle : 





6 EI 


En introduisant le paramètre sans dimensions : y = L'Gs 
F: 


et en inversant on obtient : 





Si on fixe maintenant Vj — 6j — O0, on trouve de la même 
manière : 





d’où la matrice de rigidité de l’élément de poutre en flexion : 





Fig. 52. - Matrice de rigidité d'un élément de poutre en flexion. 


Il convient maintenant de construire les matrices de rigidité 
dans le cas où, à l’une des extrémité, il y a une rotule (fig. 53). 


Ti Tj | 





Fig. 53. - Cas d’une 


Mi 
poutre avec rotule à gauche. DA Re Mj 


Dans ce cas un couple Mi est sans effets. De même une rotation 
Oi. Les seules grandeurs à considérer sont les forces Ti, Tj, Mj et les 
déplacements correspondants Vi, Vj et 0j. Les conditions d’équilibre 
sont alors : 
Ti +Ti=0 
TiL+Mj=0 
En fixant chaque fois deux des déplacements, on obtient facile- 


ment les expressions des forces en fonction du troisième. D’où la 
matrice de rigidité : 








Fig. 54. - Matrice de rigidité avec rotule à gauche. 
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De même, dans le cas d’une rotule à droite, on trouve : 





Fig. 55. - Matrice de rigidité avec rotule à droite. 


Si maintenant on a une rotule vis-à-vis de l’effort tranchant à 
gauche (fig. 56), l’élément de poutre est en flexion pure : 





Fig. 56. - Élément en flexion pure. 


On obtient dans ce cas : 





Fig. 57. - Matrice de rigidité d'un élément de poutre en flexion pure. 


Le sous-programme ’Format', dont l’organigramme est donné 
figure 58, détermine les différents coefficients de la matrice de rigidi- 
té et les range dans le tableau ’R'. 
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Kk=Liais(1,3) 
LlsLiais(1,4) 
Long=X(I+i)-X(I) , 
Gamma=( 









Den=E(Ll,£)kSect(Kk,1) 






Kigidité 
distorsionelile 


Gamna=6xE(L1,2)kSect(Kk,2)/Den/Long'à 









Barre avec 
articulation 


Aux=E(L1,2)kSect(Kk,2)/(1+2XGamma) 
RUI,1)=18XAuxx/Long"3 
R(I,8)=R(I,1)XLong/2 

R(1,3)=4K({+Gamna/2)XAux/Long 

R(T,4)=2X(1-Gamma)Xâux/Long 
















Rotule 
Effort tranchant 


Liais(1,5) = 
1090 OR { 









Rotule 
Moment 


CODE INCORRECT : STOP 






RL, L)=RCI,2)=0 
R(I,3)=E(L1,2)KSect(Kk,2)/Long 
RCI,4)=-R(I,3) 









Aux=JKE(L1,2)XSect(Kk,2)/(1+2XGamma)/Long*3 
R(I,2)=AuxkLong 
R(I,3)=AuxkLong*2 
R(I,4)=0 






Fig. 58. - Organigramme du sous-programme Format’. 


Listing du sous-programme ‘Format’ : 


3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3540 
3520 
3530 
3540 
3550 
3560 
3570 
3586 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3490 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3830 


LOL CE TU DO DT TOUTE CE TER CUT RS LESC T CT ST TR LUTTE LL TT TT TETE TT 
SUE Format(1) 
Laon 
1 


! Appelé par Prog. Princ. 
[ 


OPTION BASE 1 
l 


| Calcul des coefficients de la matrice de rigidité I 
| 

COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelen,Nmodul ,Nsect 
CON Titre$[801,Charge$(5)(80] 

COM X(100),E(5,2) ,Sect(15,2) 

COM INTEGER Liais(200,5),Isup(100) 

1 

COM Nd1b,Nn,Ntot 

COM INTEGER Iappui(50) ,Maxa(301) 

CON A(S150) ,R(200 ,24) 

l 

lisLisis(I,1) 

Jj=Liais(1,2) 

Kk=Liais(1,3) 

Li=Liais(I,4) 

l 

Long=X(I+1)-X(1) 

Ganna=0 

Den=E(Ll,i)kSect(Kk,i) 

IF Den=06 THEN 3700 
Gamna=6KE(L1,2)kSect(Kk,2)/Den/Long*2 

IF Liais(1,5)()0 THEN Artic ! Harre avec articulations 
Aux=E(L1,2)k6ect(Kk,2)/(1+2XGanmma) 
R(I,1)=12*Aux/Long"3 

R(I,2)=R(I,4)KLong/2 
R(1,3)=4KX(1+Gamma/2)XAux/Long 
R(1,4)=2K(1-Gamna)KAux/Long 

SUHEXIT 

| 
frtic: | 

IF (Liais(1,5)=1000) OR (Liais(1,52=10) THEN R{ 
IF (Liais(1,5)=100) OR (Liais(I,5)=4) THEN Ra 
PRINT “CODE LIAISON INCORRECT" 

STOP 
Ri: ! 

R(I,4)=R(1,2)=0 

R(I,3)=E(L1,2)kSect(Kk,2)/Long 

R(I,4)=-R(T,3) 

: SUBEXIT 
Re: ! 
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3890 Aux=3kE(L1,2)XSect{Kk,2)/({+Gamna/2)/Long*3 
3906 R{I,1)=Aux 

3910 R(I,2)=AuxkLong 

3920 R(I,3)=AuxXlong"£ 

3930 R(1,4)=0 

3940 SUKEXIT 

3950  SUREND 


4.6- FORMATION DE LA MATRICE DE RIGIDITÉ 
SOUS FORME STANDARD 


Le sous-programme ‘Matrig', dont l’organigramme est donné 
figure 59, reconstitue la matrice de rigidité à partir des valeurs 
précédemment calculées, sous sa forme prévue pour la formation 
de la matrice globale. 




















Rig(i)=R(I,1) 
Rig(2)=R(I,3) 
Rig(3)=R(I,2 
Rig(4)=R(I,4) 
Rig(S)=-R(I,2) 
Rig(6)=-R(I,1) 
Rig(7)=R(I,3) 
Rig(4)=-R(I,2) 
Rig(9)=R(1,4) 
Rig(40)=R(I,2) 








Rig(a)=Q 

Rig(3)=0 
“Rig(5)=0 
Rig(9)=0 


Fig. 59. - Organigramme du sous-programme ‘Matrig_. 
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Listing du sous-programme 'Matrig' 


4800 LAON 
4810 SUK Matrig(R(X),Rig(x),I,Code) 

4820 ! RDPHOPHHOPONOPENONNENNNNNONEENENENPENEOEEDENNENENENENEX 
4830 ! 

4840 ©! Appelé par Assemb 

4850 ! 

4860 OPTION HASE 1 

4870 ! 

4880 ! Calcul de la matrice de rigidité de la barre sous la forme 

4690 ! 

4900 ! Rig(i) Rig(3) Rig(b) 


4910 ! Rig(2) Rig(S) 
4920 ! Rig(4) 
4930 ! 


4940 Rig(i)=R(I,1) 

4950 Rig(2)=R({I,3) 

4960 Rig(3)=R(1,2) 

4970 Rig(4)=R(I,1) 

4980 Rig(5)=-R(1,2) 

4990 Rig(b)=-R(I,1) 

5000 Rig(7)=R(I,3) 

5010 Rig(8)=-R{I,2 

5020 Rig(9)=R(1,4) 

5030 Rig(i0)=R(1,2) 

5040 IF Code=0 THEN SUBEXIT 

5050 IF Code=100 THEN Gauche 

S060 IF Code=1 THEN Droite 

5070 SUREXIT ! Rotule au cisaillement 
5080 Gauche: ! Rotule à l’origine 
5090 Rig(2)=Rig(3)=Rig(5)=Rig(9)=0 
5100 SUBEXIT 

S410 Droite: ! Rotule à l'extrémité 
S120 Rig(7)=Rig(8)=Rig(9)=Rig(10)=0 
S130 SUBEXIT 

5140 SUBEND 


4.7 - CHARGES EN TRAVÉES 


Ces charges sont soit concentrées, soit réparties. On décompose 
ce type de chargement en charges nodales, égales aux réactions de 
la poutre bi-encastrée, changées de signe, et d’un système auto- 
équilibré, ainsi que l’illustre la figure 60. 


99 


En R; ua Ra É Ra 


Fig. 60. - Décomposition des charges en travées. 


Ce type de décomposition fut déjà rencontré pour les ossatures 
lors de la prise en compte de l’action du poids propre. Le système 
de charges auto-équilibré correspondant à un tronçon de poutre 
bi-encastrée donne des déplacements et des réactions d’appui nulles 
dans la structure. Il convient néanmoins de sauvegarder les charges 
nodales équivalentes pour le calcul des efforts aux extrémités des 
barres (ces valeurs, 4 par élément de poutre, sont conservées dans 
le tableau ’R'). 


Nous allons décrire la présentation utilisée pour le calcul des 
réactions d’appuis d’une poutre bi-encastrée, sous l’action de forces 
ou de couples concentrés puis répartis. 


Soient R1 et R2 la réaction verticale et le moment d’encastre- 
ment à gauche, R3 et R4 la réaction verticale et le moment d’encas- 
trement à droite. Soient Fy et Cz la force et le couple appliqués à 
une distance ‘a’ de l’extrémité gauche. Un calcul classique donne, 
sous forme matricielle, les valeurs suivantes pour les réactions d’appui 
dues à la force Fy et au couple Cz (fig. 61). 


CE 


| L2G+m —L#(7+1) —3L| 2L/(2+7) |-L(r+1) 
DEEE 


DER CN EST ETS 












Fig. 61. - Réactions d'appui pour une force Fy et un couple Cz. 
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Le vecteur formé par les puissances de ‘a’ dépend de la position 
de la charge concentrée, alors que les coefficients de la matrice ne 
sont fonction que de la portée et des caractéristiques du tronçon de 
poutre. Il ne seront donc pas recalculés si le chargement concerne 
toujours le même tronçon. 


Pour une charge variant linéairement (fig. 62), entre les abscisses 
a’ et ‘b'" où les composantes de la densité de charge sont respective- 
ment Pa, Ca, Pb, et Cb, l'expression de P entre a et b est donnée par : 


x Pb 





Fig. 62. - Densité de charge entre a et b. 


On remplacera les puissances de ‘a’, a" par(b"*! —ag"*1)/(n+ 1) 
et on fera la somme de la partie constante et la partie linéairement 
variable. 


Dans le programme, on a posé : 


Coef1=(Ph—Pa)/(b—a) 
Coef2=(Ch—Ca)/(b—a) 
Coef3=(b*Pa—a*Pb)/(b—a) 
Coef4=(b+Ca—axCb)/(b—a) 


Dans le cas où on est en présence d’une articulation à l’extré- 
mité de la poutre, il faut utiliser les matrices (fig. 63 et 64). 
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Fig. 63. - Rotule à gauche. 














Fig. 64. - Rotule à droite. 


Pour un élément avec rotule, vis-à-vis de l’effort tranchant on 
utilise les matrices figures 65 et 66. 
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Fig. 65. - Effort tranchant nul à gauche. 














R3 


F,>0 





Effort tranchant nul à droite. 


Fig. 66. 
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4.8 - DÉFINITION DES CHARGEMENTS 


Quatre types de chargement ont été prévus : 


— Les charges directement appliquées aux nœuds. C’est évi- 
dement le cas le plus simple et les valeurs algébriques de ces charges 
sont directement cumulées dans le second membre. 


— Les charges appliquées sur les travées. Celles-ci peuvent être 
concentrées ou réparties, ce qui correspond aux valeurs ‘Index'=2 
et 'Index'=3. 

— Les déplacements d’appuis. La séquence correspondante 
est identique à celle reproduite pour les ossatures, et l’organigramme 
n’en sera pas donné. 


Indice=Maxa(Nn+i)-Ni 















1c=( 
Itrav=0 


Déplacements d'appuis 
Travée en cours de 
chargement 

Remise à Ü0 du tableau 
des déplacements imposés 


MAT D=(0) 


Indice=Indice+Nn 





Index ? 





CALL Genererdepl 


Fig. 67. 
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L’organigramme général du sous-programme ‘Second’ est donné 
figure 67. 









Calcul des coefficients 
de la matrice transfer 


MAT Tr=0 
Long=X(1+1)-X(1) 
Longä=Long*2 
Long3=Long*3 
Kk=Liais(li,3) 
LisLiais(li,4) 
Code=Liais(1i,5) 
Ganna=0 
Den=E(Kk,1)kSect(L1,1) 









Rigidité de distorsion 





Ganna=b*E(Kk,2)XSect(L1,2)/Den/Long2 






Code=10 






Code incorrect 
STOP 











Normal | Voir 
Listing 


Den=(2+Gamma)XLong3 
Den=(1+24Gamna)kLonga 












Tr€K,L)=TrK,L)/Den 
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Aa,Fy,Ci? 


oui 
| Aa ,Pa,Ca? 
cb? 


Eb,Fb,Cb? 





F(K)=Aa"(4-K) 










Coefi=(Pb-Pa)/(Bb-Aa) 
Coef2=(Cb-Ca)/(Kb-ña) 
Coef3=(EbkP a-AakPb)/(Eb-Aa) 
Coef4=(HbkCa-AaxCb)/(Hb-Aa 


C(L)=(Hb*(b-L)-Aa"(8-L))/(6-L) 


F(LI=C(L)KCoefi+C(L+i)kCoefs 
F(L+4)=C(L)kCoef2+C(L+i)kCoef4 


MAT Fnod=Trxf 
K=IX4 Sauvegarde des 


us) forces aux extremités 
RCT1,K+L)rR(I1,K+L) de là poutre 
+FiL) 


K=Indice+Nalk(I1-1)-1 
Que) Cumul dans le 
second membre 


A(K#L)=AÎK+L)-Fnodil) 


@) 





FCK+4)=F(K)XCz 
F(KI=F(KIKF y 
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5650 
5660 
5670 
5680 
5690 
5700 
5740 
5720 
5730 

740 
5750 
5760 
5770 
5780 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
S860 
5870 
5860 


5890 


5900 
5910 






ACL+K)=ACLEK)+F (K) 


Inoeud=Nnoeuds 





ML TCCCER TT EE ETC ETS TETE TT TT PTT TT STE TT TT TT TT TEST TT TT TT TT IT TEE TS 
SUE Second 

Lx NNHHNNNNNHNNHNHNNRNHRMNNNNHNNENNENNNNERNNONNNHRNONHENNNRRENNNONENIENIONENX 
l 

| Appelé par Prog. Princ. 

! Appelle : Genererdepl 

! 


OPTION BASE { 

| 

! Introductiun du chargement, reproduction des charges 
let constitution du (ou des) seconds membres du système 
1 

COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelem,Nmodul ,Nsect 

COM Titre$[80] ,Charges$(5)(80) 

COM X(100),E(5,2),Sect(15,2) 

CON INTEGER Liais(200,5),Isup(100) 

| 

COM Nd1b,Nn,Ntot 

COM INTEGER Iappui(50) ,Maxa(301) 

COM A(S150) ,R(200 ,24) 


DIM F(8) ,Coef(3),D(300) ,Tr(4,8),Fnod(4) 
l 


Indice=Maxa(Nn+i)-Nn 

FOR 1={ TO Nchar 
Ic=0 | Paramètre indiquant les déplacements d’appuis 
Itrav=Q 
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5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6096 
6100 
b110 
6120 
6130 
6140 
650 
6160 
6170 
6180 
6190 
6200 
beit 
62et 
6230 
6240 
6250 
b260 
6270 
6280 
6290 
6300 
6310 
b3e0 
633 

6340 
6350 
6360 
6370 
6380 
6390 
6400 
6410 
b4et 
6430 


Ci: 


eo 
ra 


MAT D=(0) 

Indice=Indice+Nn ! Indice du ier Coef. du vecteur second membre 
PRINT LIN(3) 

PRINT "CAS DE CHARGE N° ";l 

PRINT L 

PRINT 

INPUT " Titre du cas de charge ?",Charges#$(1) 
PRINT Charges#$(I) 

PRINT 

Index=5 

DISP " NATURE DU CHARGEMENT ?" 

DISP " 4 : Charges concentrée aux noeuds” 
DISP " à : Charges concentrées en travées" 
DISP " 3 : Charges linéaires en travées" 
DISP " 4 : Déplacement d'appui" 

INPUT * S : FIN", Index 

ON Index GOTO Ci,C2,C2,C4 ,Nexti 

[ 


|! Charges concentrées aux noeuds 
i 
PRINT "CHARGES AUX NOEUDS" 
PRINT " N° noeud Fy Ci" 
DISP "La description du chargement est achevée par la donnée" 
DISP "de la charge sur le dernier noeud” 
INPUT "N° noeud, Fy, Cz",Inoeud,F(1),F(2) 
PRINT USING 6180 ;Inoeud ,F(1),F(2) 
IMAGE X,DDD,2X,MDDD.DDD ,2X ,MDDD .DDD 
L=Indice+Ndlk(Inceud-4)-1 
FOR K=1 TO Ndl 
ACL4K)=ACL+K)+F(K) 
NEXT K 
IF Inoeud=Nnoeuds THEN 6010 
GOTO 6440 ! Introduction de La charge sur le noeud suivant 
| 


| Charges en travée 

| 

INPUT " NUMERO DE LA TRAVEE ?",Ii 

PRINT LIN(1) 

PRINT " TRAVEE N° ";li 

IF lis=ltrav THEN Suite ! Ne pas recalculer les coef de Tr 
Itrav=li 

i 


! Calcul des coefficients de la matrice transfer 
1 


MAT Tr=(0) 

| 
Long=X(1i+i)-X(11) 
Long2=Long*2 
Long3=Long2kLong 
Kk=Liais(1i,3) 
Li=liais(1i,4) 
Code=Liais(li,5) 


6440 
6450 
b4oû 
6470 
6480 
6490 
6500 
6510 
6520 
6530 
6540 
6550 
6560 
6576 
6580 
6590 
6600 


6640 


b620 
6630 
b4l 
6650 
6bbl 
bb70 
6680 
6690 
6700 
6740 
6720 
6730 
b740 
6750 
6760 
6770 
6780 
6790 
6800 
6810 
6820 
6830 
6840 
6850 
6860 
6870 
6880 
6890 
6900 
6910 
6920 
6930 
6940 
6950 


Ganna=0 
Den=E(Kk,i)#Sect(Ll,i) 

IF Den: THEN 6460 
Ganma=6kE(Kk,2)kSect(Ll,2)/Den/Long*e 
IF Code=0 THEN Normal 

IF Code=1000 THEN Tgauche 
1F Code=i100 THEN Mgauche 
IF Code=10 THEN Tdroit 

IF Code=i THEN Mdroit 
DISP "CODE INCORKECT" 
STOP 


Normal :! 


Den=(1+2XGamma)kLong3 
Tr(4,1)=-2 
Tr(4,2)=34Long 
Tr(1,3)=2XGanmakLongé 
Tr(1,4)=-(4+2XGamma)XLongs 
Tr(4,6)=-6 
Tr(14,7)=6xLong 
Tr(1,8)=2XGanmaklonge 
Tr(2,1)=-Long 
Tr(2,2)=(2+Ganmna)XLonge 
Tr(2,3)=-(i+Ganna)XLongs 
Tr(2,6)=-3#Long 
Tr(2,7)=2XTr(2,2) 
Tr(2,8)=Tr(2,3) 
Tr(3,1)=2 

Tr(3,2)=-Tr (4,2) 
Fr(3,3)=-Tr (4,3) 
Tr(3,6)=6 

Tr(3,7)=-Tr (4,7) 
Tr(3,8)=-Tr(1,8) 
Tr(4,1)=-Long 
Tr(4,2)=(1-Gamma)KLongé 
1r(4,3)=GammakLongs 
Tr(4,6)=Tr (2,6) 
Tr(4,7)=2XTr (4,2) 
Tr(4,8)=Tr(4,3) 

GOTO Divise 


Tgauche:! 


Tr(2,8)=.5/Long 
Tr(2,3)=-1 
Tr(2,4)=Long/2 
Tr(2,7)=4/Long 
Tr(2,8)=Tr(3,4)-1 
Tr(4,2)=-Tr(2,2) 
Tr(4,4)=Tr(2,4) 
Tr(4,7)=-Tr(2,7) 
GOTO Suite 


Mgauche:! 


Den=(2+Gamma)kLongz 
Tr(i,1)=-1 
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6960 
6970 
6780 
6990 
7000 
7040 
7020 
7030 
7040 
7050 
7060 
7070 
7080 
7090 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
7170 
7166 
7490 
7e00 
72410 
7eel 
7230 
7246 
7350 
7èbû 
7270 
7280 
7290 
7300 
7310 
7320 
7530 
7340 
7350 
7360 
7370 
7380 
7390 
7460 
7410 
7420 
7430 
7440 
7450 
746 
7474 
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Tr(1,3)=(3+Gamma)kLonge 
Tr(1,4)=-(2+6amna)XLongs 
Tr(i,6)=-3 
Tr(4,8)=Tr(1,3) 
Tr(3,4)=4 
Tr(3,3)=-Tr(41,3) 
Tr(4,6)=3 

Tr(3,8)=-Tr (1,8) 
Tr(4,1)=-Long 
Tr(4,3)=Long3 
Tr(4,6)=-3XLong 
Tr(4,8)=Long3 

GOTO Divise 


Tdroit:! 


Tr(1,4)=Tr(2,3)=1r(2,8)=-1 
Tr(2,2)=.5/Long 
Tr(2,7)=i/Long 
Tr(4,8)=-Tr(2,2) 
Tr(4,7)=-Tr(2,7) 

COTO Suite 


Mdroit:! 


Den=(2+6Gamma)XkLong3 
Fr44,4)=-i 
Tr(4,8)=3XLong 
Tr(1,3)=GammakLongz 
Tr(4,4)=-(2+Gamma)kLongs 
Tr(4,6)=-3 
Tr(4,7)=6XLong 
Tr(4,8)=Tr(41,3) 
Tr(2,4)=-Long 
Tr(2,2)=3KLonge 
Tr(8,3)=-2#Long3 
Tr(2,6)=-Tr(1,2) 
Tr(2,7)=6*Longe 
Tr(2,8)=Tr(2,3) 
Tr(3,1)=1 
Tr(3,2)=-Tr(1,2) 
Tr(3,3)=-Tr (41,3) 
Tr(3,6)=3 
Tr(3,7)=-7r (1,7) 
Tr(3,8)=Tr(3,3) 


Divise:! 


FOR K=1 TO 4 
FOR L={ TO 8 
.Tr(K,L)=Tr(K,L)/Den 
NEXT L 
NEXT K 


Suite:! 


IF Index=3 THEN Repartie 

INPUT " Abscisse de la charge, Fy, C1 ?",Aa,Fy,Cz 
PRINT " Abscisse = "Aa," Fy = ";Fy," C1 = ";Cz 
FOR K=1 TO 4 


7480 
7490 
7500 
7510 
7520 
7530 
754 
7550 
7560 
7570 
7580 


7590 
7600 


7640 
7620 
7630 
7640 
7650 
7660 
7670 
7680 
7691 
77u0 
7740 
7720 
174 
T7AÙ 
7750 
7780 
777 
7780 
7790 
7800 
7810 
7820 
7830 
7840 
7850 
7860 
7870 
7880 
7890 
7506 
7940 
7920 
7940 
7940 
7950 
7960 
7970 
7986 
7990 


F(K)=Aa" (4-K) 
FUK+4)=F (K)XCz 
F(KI2F CIOKFy 
NEXT K 
GOTO Produit 


Repartie:! 


DISP " Abscisse origine, Densité de forces, Densité de couples ?" 
INPUT Aa ,Pa,Ca 

PRINT "Abscisse origine............. 
PRINT "Densité de forces origine...= ";Pa 
PRINT "Densité de couples origine...= ";Ca 


DISP " Abscisse extrémité, Densité de forces, Densité de couples ?" 
INPUT Eb,Pb,Cb 


PRINT "Abscisse extrémité........... "Bb 
PRINT "Densité de forces extrémité..= “;Pb 
PRINT "Densité de couples extrémité = ";Cb 
PRINT 
IF Eb)Aa THEN 7680 
PRINT "ERREUR" 
GOTO Repartie 
Coefi=(Pb-Pa)/(Eb-Aa) 
Coef2=(Cb-Ca)/(Kb-Aa) 
Coef3=(HbXFa-AakPb)/(Eb-Aa ) 
Coef4={EbkCa-AakCb)/(Bb-Aa) 
FOR L=4 TO S 
C(L)=(Kb*(6-L)-Aa"(6-L))/(b-L) 
NEXT L 
FOR L=i TO 4 
F(L)=C(L)XCoef{+C(L+{4)4Coef3 
F(L+4)=C(L)kCoef2+C(L+{)kCoef4 
NEXT L 


";Aa 


Produit:! 


C4: 


MAT Fnod=TràF 
[ 


Sauvegarde 

[ 

K=IX4 

FOR L={ TO 4 
R{Ii,K#L)=R(Ti,K+#L)+Fnod(l) 

NEXT L 

! 

| Cumul dans le second membre 

K=Indice+Ndlt(1i-1)-1 

FOR L={ TO 4 
A(K+L)=A(K+L)-Fnod(l} 

NEXT L 

GOTO 6040 

1 


! Déplacements d'’appuis 
l 


PRINT "DEPLACEMENTS IMPOSES DES APPUIS" 
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8000 
80410 
8020 
8030 
8040 
8050 
8060 
8070 
8080 
8090 
8100 
8110 
8120 
830 
8140 
8150 
8160 
8170 
8180 
8490 
8200 


l 
INPUT " N° du Noeud ?",N 
IF Isup(N)=0 THEN Erreur 
Ic=lc+i 
Nombre=lsup(N) 
CALL Decode (Nombre ;Nd1 ,Coef (k)) 
INPUT "DY, Teta ?2",F(4),F(2) 
PRINT "Noeud N° ";N," Dy = ";F(4)," Teta = ";F(2) 
Rep$="Qui" 
In=(N-£)XNd1 
l 
| Seuls les déplacements imposés correspondants à des degrés de 
|! liberté bloqués sont pris en compte et placés dans le vecteur D 
| 
FOR Jj=i TO Ndl 
IF Coef(Jj)=0 THEN 8170 
D(In+Jj)=F(J;j) 
NEXT Jj 
INPUT "Autre déplacement d'appui (Qui) ?",Rep$ 
IF Rep#="Oui" THEN 8040 
GOTO 6010 


8210 Erreur: DISP "Le noeud n’est pas un appui" 


8220 


GOTO 8010 


8230 Nexti:! 


8240 
8250 
8260 
8270 
8280 
8290 
8300 
8310 


! Modification du second membre en fonction des déplacements d'appuis 
CALL Genererdepl(Indice,Ic,Nn,Ndlb,DCX) ,ACK) ,Lappui(x),Maxa(x)) 
NEXT I! 


SUREXIT 
SUBEND 


4.9 - IMPRESSION DES DÉPLACEMENTS 


Cette séquence est simplement particularisée par des labels 
adéquats. Il convient de rappeler le point important suivant : si une 


barre est articulée à une de ses extrémités, le déplacement relatif au 


nœud concerné se rapporte à l’extrémité de la barre adjacente. 


Comme il revient au même de considérer l’articulation pour la barre 


de droite ou de gauche (mais surtout pas pour les deux!), on peut 
trouver deux déplacements différents pour le nœud. Ce fait ne 
doit pas étonner. 
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Le listing du sous-programme ‘Depl' : 


10260 1 KE HOONHNNNNNNNNNNNNNNNNNNNNNNANNNNNNNNNNNNNNNNCNNNNNNNONNNNN 
10270 SUE Depl(Indice,Nd1 ,Nnoeuds ,A(X)) 
10280 1 HopHHNNNNNNHNNNNNNNHHNNNNHRMNNNNNNNNNNEEHNNENEINNNNNNEEEX 
10290 ! 
40300 ! Appelé par Prog. Princ. 

i 


10340 ! 

10320 OPTION RASE { 

10330 : 

10340 ! Impression des déplacements 
10350 ! 


10360 PRINT LIN(1) 

40370 PRINT 'DEPLACEMENTS DES NOEUDS" 
10380 PRINT 

10390 PRINT " N° Noeud Dy Teta" 
10409 PRINT 

10410 FOR I=1 TO Nnoeuds 

10420 li=Indice+(I-1)XNd1 

10430 PRINT USING 40440;1,A(11),A(Ii+i) 

10440 IMAGE 3X,DDD,5X,MD.DDDDDDE ,8X,MD.DDDDDDE 
10450 NEXT I 

10460 SUBEXIT 

40470 SUREND 


4.10 - EFFORTS ET RÉACTIONS D'APPUI 


Les efforts aux extrémités des barres sont calculés à partir des 
matrices de rigidité (le sous-programme "Matrig’ est appelé de nou- 
veau) et des déplacements. Il convient de rajouter les réactions d’ap- 
pui pour le tronçon supposé parfaitement encastré qui ont été 
sauvegardées dans le tableau ‘R’. La séquence de calcul des réactions 
d'appui de la poutre continue est identique à celle décrite pour les 
ossatures. 


L’organigramme du sous-programme 'Efforts' est donné figure 
68. 
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li=liais(l,1) 
Jj=Liais(1,2) 
Code=Liais(l,5) 






CALL Matrig Formation de la 
K=Indice+Ndlk(I-1) 


matrice de rigidité 





F(AI=RIQUI)KACKIHRIUIIKAKELIHRIGUH)IKALK HZ) RIAD) KACKEI)HRET , Li) 
F(2)=R1g(3)KA(K)+Rig(2)KACK+4)+R1QCS)KACK+Z)+R1g(P)KACKEIIHRT ,L+2) 
F(3)=R1g(6)KACK)+RIg(S)KACK+L)+R1QADKAK+2)+RIgBIKACKEIIHR(T ,L+3) 
F(4)=R1g(40)KA(K)+R1g(7)XACK+L)+RIGCBIKAK+2I+RIg7)KACKEI)HRUT ,L+4) 


Jj=liais(1,J) 


Nombre=Isup(J; 
CALL Decode 


NATE 
: oui 
er l 


Reac(Jj,K)=Reac(Jj,K)+F(2K(J-4),K) 


0 









M 1,11,3j,F(4),F{2),F(3),F (4) 





Fig. 68. - Organigramme du sous-programme 'Efforts. 
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Listing du sous-programme 'Efforts’ : 


10480 
10490 
10500 
105410 
10520 
10530 
10540 
10550 
10560 
10570 
10580 
10590 
10600 
10610 
10620 
106630 
10640 
10650 
10660 
10670 
10660 
10690 
10700 
10710 


16720 
10730 
10740 
10750 
10760 
10770 
10780 
10790 
10800 
10810 
10820 
10830 
10840 
10850 
10860 
10870 
10880 
10890 
10900 
10910 
10920 
10930 
10940 
10950 


MR LLLOTE EEE TES TT TETE TT TETE TETE TES TETE TT TT ET TITI TETE TTTTTEETEE I 
SUE Efforts(Indice,Ichar) 
LPO 
! 
! Appelé par Prog. Princ. 
Leu 
OPTION HASE 
| 
| Calcul et impression des efforts dans les barres 
| 
COM Nnoeuds ,Nbarres ,Nchar ,Nd1 ,Nnelem,Nmodul ,Nsect 
CON Titre$[80],Charges$(S)(80] 
COM X(400) ,E(5S,2) ,Sect(15,2) 
CON INTEGER Liais(200,5),Isup(100) 
| 
COM Ndib,Nn,Ntot 
COM INTEGER Iappui(50),Maxa(301) 
COM A(5150) ,R(200,24) 
! 
DIM Rig(24),Coef(2),F(4) ,Reac(50,8 
PRINT LIN(1) 
PRONT'EFFORTS DANS LES BARRES" 
PRINT 
PRINT " N° Karre N° Orig. N° Extr. T Orig. M Orig. T Extr. 
H Extr." 
PRINT 
FOR 1=1 TO Nbarres 
li=Liais(1,1) 
Jj=Liais(I,2) 
Code=Liais(1,5) 
l 


CALL Matrig(R(W) ,Rig(X), I, Code) 
1 


K=Indice+Ndlk(I-1) 

L=Ichark4 
F(LI=RIgCLIKACKI4RIIZ)KACK#LIHRIQHIKAK+2)HRIGLO)KACK+ZI HR ,L+1) 
F(2)=Rig(3)XA(K)+Rig(2)KACK+L)+RIg(SIKACK+2)+RIIPIKACKE+I)HR(I ,L+2) 
F(D=RIg(b)KACK)+R1g(S)KA(K+I)+RIg(A)KACK+2)+RIG(BIKACK+IIHR(I ,L+3) 
F(4)=Rig(i0)KA(K)+Rig(9)KA(K+L)+Rig(B)KA(K+2)+Rig(7)KA(K+II4R(I,L+4) 
1 


| Calcul des Réactions d’appuis. 
j 
FOR J=1 TO Nnelen 
Jj=Liais(1,J) 
IF Isup(Jj)=0 THEN 10980 
Nombre=Isup(J;j) 
CALL Decode(Nonbre ,Nd1,Coef (X)) 
FOR K=1 TO Nd1 
IF Coef(K)=0 THEN 10979 
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10960 Reac(Jj,K)=Reac(Jj,K)+F(2K(J-1)#K) 


10970 NEXT K 

10980 NEXT J 

10990 | 

31000 ! Impression des efforts 
44010 | 


11020 PRINT USING 41030;1,1i,1j,F(4),F(2),F(3),F(4) 
11030 IMAGE 4X,DDD,68X,DDD,BX,DDD,4( X,MDDDDDD .DDD,X) 
44040 NEXT I 

11050 ! 

11060 ! Impression des Réactions d'appuis 

11070 ! 

11980 PRINT LIN({) 

44090 PRINT 'REACTIONS D'’'APPUIS" 

41100 PRINT 

44410 PRINT " N° Noeud Ry Mz° 

11120 PRINT 

15430 FOR I=1 TO Nnoeuds 

41140 IF Isup(1)=0 THEN 14170 

44150 PRINT USING 11160 ;1,Reac(1,i),Reac(I,2) 

41160 IMAGE 3X,DDD,2(2X,MD.DDDDDDE ) 

44170 NEXT I 


11180 SUBEXIT 
11190 SUREND 


4.11 - EXEMPLE DE CALCUL DE POUTRE CONTINUE 


La poutre est définie par la figure 69. 





Fig. 69. - Exemple poutre continue 
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TITRE : Essai poutre continue 3 Appuis 


Poutre continue, axe Ÿ vertical 


Nombre de noeuds = {1 
Nombre de barres = 10 

Nombre de cas de charges = 
Nombre de sections différentes 
Nombre de modules d’élasticité 


ARSCISSES DES NOEUDS 


N° X. Code appui 
{ 0.000 10 
a 1.000 û 
3 2.000 û 
4 3.000 (l 
5 4.000 û 
6 5.000 10 
7 6.000 (1 
8 7.000 0 
9 8.000 0 

10 9.000 û 
ai 10.000 ii 


DEFINITION DES LIAISONS 


N° Origine Extrémité N° Sect. N° Module Code articulation 


Î { ë Î { û 
a è 3 { { (l 
3 3 4 { { û 
4 4 5 Î { û 
5 5 6 { { û 
6 b 7 { î û 
7 7 8 i { û 
8 8 9 i { (t 
9 9 10 { { û 
40 10 ii { 1 Î 
CARACTERISTIQUES DES KARRES 
Section réduite N° 4 = .000S 


Inertie N° 41 = .0005 


MODULES D’ELASTICITE G et E 


Module G N° 


Î 800000 
Module E N° 


2000000 
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TRAVEE N° 
Abscisse = .75 Fy = -10 Cz = 
TRAVEE N° 10 
Abscisse = .25 Fy = -10 C1 
RESULTATS 
CAS DE CHARGE N° 1 
Test synétrie 
DEPLACEMENTS DES NŒUDS 
N° Noeud Dy Teta 
{ 0.000000E+00 -3.545673E-04 
2 5.699038E-03  2.15961SE-03 
3 S.410817E-03  b.39663SE-04 
4 3.929327E-03 -2.269234E-04 
S 1.907933E-03 -4.401442E-04 
6 0.000000E+00 -8.000000E-15 
7 4.907933E-03  4.401442E-04 
8 3.929327E-03  2.269231E-04 
9 5.410817E-03 -b.39663SE-04 
10 5.699038E-03 -2.159615E-03 
qi 0.000000E+00  0.000000E+00 
EFFORTS DANS LES KARRES 
N° Earre N° Orig. N° Extr. T Orig. M Orig. 
{ { 2 410.653 0.000 
è a 3 .653 1.847 
3 3 4 .653 1.193 
4 4 S .b53 .540 
5 5 6 .653 -.113 
6 b 7 -.653 -.767 
7 7 8 -.653 -.113 
8 8 9 -.653 .540 
9 9 A0 -.653 1.193 
10 10 ii —.653 1.847 
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CAS DE CHARGE N° 


Test symétrie 


1 


10 


-10 


ï 


Extr. 


.653 
.653 
.653 
.bS3 
653 
653 
.b53 
.653 
.653 
40. 


653 


M Extr. 


-1.847 
-1.193 
-.540 
143 
.767 
443 
-.540 
-1.193 
-1.847 
0.000 


REACTIONS D'APPUI 
N° Noeud Ry Mz 
4 1.065337E+01  0.000000E+00 


b -1.306731E+00 0.000000E+00 
41 1.065337E+01 0.000000€+00 
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ÉLASTICITÉ 
PLANE 


5.1 - BASES THÉORIQUES 


Le problème est résolu par la connaissance des deux compo- 
santes du déplacement d’un point, u (x, y) et » (x, y). On déduit 
alors les déformations par : 


_ ou 
Ts 

_Ôy 
LT 

_0u Ôv 
T3 + 


Les déformations et les contraintes sont reliées par la loi de 
HOOKE : 
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pour un état de contrainte plane (0, =0, €, ()0) et par 


. E(1-—) 





pour un état de déformation plane (0, (0 et e, =0). 
Cette relation peut être représentée matriciellement par : 
o=[Dle 


5.2- APPROXIMATION DES FONCTIONS DE DÉPLACEMENT 


La connaissance exacte des composantes du déplacement en 
tout point étant impossible, on cherchera à déterminer leurs valeurs 
’en un certain nombre de points’. A cette fin, on découpe le domaine 
étudié en éléments triangulaires, et on cherche à interpoler les 
déplacements à l’intérieur du triangle à partir des déplacements 
des sommets. Le plus simple est de supposer que les fonctions 
déplacements u(x,y) et v(x,y) sont des fonctions linéaires à l’inté- 
rieur de chaque triangle élémentaire : 

u (x, y) = Got à x+ @ y 


vx, Y) = &3 + Ga x + Os y 


u (x, y} et v (x, y) prennent les valeurs U;, V;, U,, V;, Uz et Vz aux 
sommets d’un triangle i, j, k. On peut ainsi déterminer les constantes 


QXo 1 A2 A3 Aa ds 


en fonction des coordonnées des sommets i, j, k du triangle. 
On obtient pour u et y les expressions suivantes : 


122 


1 
u (x, y) s [@;+b;x+ Gy)U;+ (a; +b;x+c;y)U;+ (ax +bxx+cx y)Ur] 


1 
v(x, y) later a MN +(a;+ bjx+c;y)V;+ (ax + bxx + cx y) Vx] 


GE XYR XKYS : XVI XIYk dk Xi}; Xj Yi 
bi=y;— yk bj=yK— Yi bx=yi- y; 
Ci Xk— Xj CG Xi Xk Ck Xi Xi 


À = Surface du triangle ijk. 


Ces notations seront également utilisées lors de l’étude du pro- 
blème de torsion. 


5.3 - MATRICE DE RIGIDITÉ D'UN ÉLÉMENT TRIANGULAIRE 


Par application des opérateurs différentiels reliant les déforma- 
tions aux déplacements, on obtient la relation suivante : 





que l’on écrit matriciellement sous la forme 
e=[B]6ô 


Cherchons la relation entre les forces appliquées aux nœuds à, j, 
k d’un triangle élémentaire et les déplacements de ses sommets. 


On considère un déplacement virtuel 8* des sommets à, j, k. Il 
en résulte des déformations virtuelles 


e* = [B] ô* 
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Le théorème des travaux virtuels permet d’écrire : 


16*F = [ le* o dv 


y 
— L'intégration porte sur le volume de l’élément 
— !ô est un vecteur ligne : { Ui Vi Uj V; Ux Vx| 
— Fest un vecteur colonne dont les composantes sont : 
j {Ex Fri F,; Fy; Fxk Fyx } 


En utilisant les relations précédentes, on obtient : 


ô*F = | ‘6*'B D B6 dv 


y 


et par conséquent 
F = | ‘BDBdv | ô 


y 


Comme les coefficients des matrices sous le signe somme 
sont constants, l'intégration consiste à multiplier ce produit par 
le volume de l’élément. On pose 


K='BDB(Ae) 


matrice de rigidité de l’élément triangulaire dont les déformations 
sont approchées par des fonctions linéaires. 


Sous forme explicite, les coefficients de la matrice de rigidité 
de l’élément triangulaire sont donnés, figure 70 : 


+ l oies vbic;+ ab bibx 41 ac vbicx + 5" cibk 
12e = I=v lv 
‘ F Eu | a EI 


bib; oc 


= a. 
€ 
v 2 1+v 
bi be 'ejck verbe + hier b; be + cjek vojbe + bjek + 7 bkck 
mr Fra rh ne M to 7 


Fig. 70. - Matrice de rigidité en contrainte plane. 























vbic+} cb) cicj+ 
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5.4 - PROGRAMME PRINCIPAL 


Dans le programme principal, nous avons introduit la variable 
’Nat$’ qui définit la nature du problème : contrainte plan ou défor- 
mation plane. Cinq cas de charge ont été prévus. Le tableau 'E' 
définit le module de YOUNG et le coefficient de POISSON du maté- 
riau. Le tableau "Epais', l’épaisseur et la densité du matériau, afin 
de pouvoir générer éventuellement un chargement équivalent au 
poids propre. Le listing du programme principal est reproduit ci- 
dessous : 


{0 1 OHONNNNNNENHANNHNNNNNNNNNNNNNNHNNRNNNNNRHRNNNNHINNNNNNNNHNNNEX 
20 ! * 
CONTRAINTE PLANE OÙ DEFORMATION PLANE 


PAR LA METHODE DES DEPLACEMENTS 


* XX 0000 XXX PONRX 
x * X * X  * x 

*X * * 

* HONX * 

* * x * 
* * *X L 
HOONERX HHERX * 


Kerr x 


; 
* * 
* * 
* * 
* * 
* * 
UX * 
EE: * 
Li * 
* * 
* * 
* * 
* * 
* *X 


XPH 


en 
LA 
1 


| BLOC COMMUN 
180 OPTION BASE À 
490  PRINTER IS 7,1 ,MIDTH(130) 


200 ! : 

240 ! Définition des variables 

eë0 ©! 

230 ! 1°/ VARIABLES DEFINISSANT LE PROBLEME 

ë40 ! 

250 ! Nat$.............. Nature du probleme 

260 ! Nat$="C" Contrainte plane 

270 | Nat$="D" Déformation plane 

880 ! Nnoeuds........... Nombre de noeuds 

290 ! Nelem............. Nonbre d'éléments 

300 ! Nnelem............ Nombre de noeuds par élément 

340 ! Nchar............. Nombre de cas de charges 

Set. “UNdE.ssieuhu.n Nombre de degrés de liberté par noeud 
330 ! Ndl = & 

340 ! Nnodul............ Nombre de modules d’élasticité et de coef. de POISSON 
350 ! Nepais............ Nombre d'épaisseurs diffèrentes 

360 ! Titre$............ Titre du probleme 

370 ! Chargest.......... Titres des cas de charge 

380 1! XCI), YCI)........ Coordonnées du noeud I 
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390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
540 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
740 
720 
73 
740 
750 
760 
770 
780 
790 
800 
810 
826 
830 
840 
850 
860 
870 
880 
890 
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Isup(1) .......... Code appui 

00 Noeud libre 

10 Noeud bloqué suivant X 
4 Noeud bloqué suivant Ÿ 


Autre combinaison : 14 


un 


| 
| 
Ë 
j 
! 


| Liais(1,1)........ Numéro du i er noeud de l'élément 

! Liais(I,2)........ Numéro du 2 ème noeud de l’élénent 

| Liais(1,3)........ Numéro du 3 ème noeud de l'élément 

|! Liais(l,4)........ Numéro du module d’élasticité E et coef POISSON 
| Liais(1,5)........ Numéro de l'épaisseur et densité 


COM Nnoeuds ,Nelem,Nchar ,Nd1,Nnelen,Nnodul ,Nepais 
COM Titre$[801,Charge#$(S)1801,Nat$ 

COM X(400),Y(400),E(S,2) ,Epais(15,2) 

COM INTEGER Liais(200,5),Isup(400) 

1 


| 2°/ VARIABLES AUXILIAIRES EN COMMUN 


Ü Ndib.............. Nombre de degrés de liberté fixés (appuis) 

DNNian oups énsteréeusst NoeudskNdl : Dimension du système linéaire 
YONtot.............. Nombre de coefficients de la matrice À 

D Ipp(i)............ ={ Si le Ième cas de charge comprend le poids propre 
! =0 Si non 

PLappui(T) ........ Numéros des degrés de liberté fixés 

| Maxa(I)........... Adresses des termes diagonaux de la matrice 

PACED Sie ar Tableau des coefficients de la matrice et des 

! seconds membres 

RCD ES D. Leu de 1 à 10 Car. des éléments 


| de {1 à 15 les déformations dues à ia température 
COM Ndib,kn,Ntot, Ipp(5) 

CON INTEGER Iappui(S0) ,Maxa(301) 

COM A(S150) ,R(E00, 15) 

| 

| Initialisation 
l 

MAT Liais=(0) ! 
MAT Isup=(0) ! 
MAT Maxa=(0)  ! 
MAT A=(0) ! 
MAT R=(0) ! 
MAT Ipp=(0} ! 
Nepais={ 

Nmodul=i 

Nnelen=3 | CAR ELEMENTS TRIANGULAIRES 
Ndi=2 ! DEUX ddl PAR NOEUDS 

! 


\ Lecture des données générales et des coordonnées des noeuds 
i 


CALL Lect 
} 


! Calcul des indices des degrés de liberté bloqués 
! 


CAR VARIABLES EN COMMUN 


900 

910 

get 

930 

940 

950 

960 

970 

980 

990 

1000 
4010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
{140 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
4210 
1280 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1330 
1390 
1400 
1410 


CALL Indapp(Nnoeuds,Nd1 ,Ndlb ,Isup(x) ,Iappui(x)) 
! 


L Reproduction des données relatives à la structure 
1 


CALL Imprin 
[ 


FOR 1=4 TO Nelen 
1 


! Formation de la matrice de rigidite de la liaison I 
[ 


CALL Format(1) 
I 


| Calcul de la hauteur des colonnes correspondantes 
[ 


CALL Hcol(I,Ndl,Nnelem,Liais(*) ,Maxa(x)) 
NEXT 1 
fais 
l Calcul des indices des termes diagonaux 
| 
CALL Indiag{Nn,Ntot,Maxa(x)) 
l 


| Assemblage de la matrice globale 
! 


CALL Assemb 
l 


! Lecture des charges et création du ou des seconds membres 
1 


CALL Second 
l 


Modification de la matrice en fonction des appuis 
! 

CALL Modif (Nn,Nd1b,ACK) ,Maxa(k) ,Lappui(xk)) 

! 


| Résolution du système linéaire 
| 


CALL Solve(Nn,Ntot,Nchar ,ACK) ,Maxa(X)) 
! 
| Pour chaque cas de charge, impression des déplacements des noeuds, 
! calcul et impression des contraintes et des réactions d'appuis. 
| 
Indice=Maxa(Nn+i)-Nn ! permet de récuperer les déplacements de chaque 
| cas de charge dans le tableau À. 
PRINT LIN(2) 
PRINT "RESULTATS" 
PRINT 
FOR 1=4 TO Nchar 
! 
PRINT 
PRINT "CAS DE CHARGE N° ";f 


PRINT Charge$(1) 
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1420 | 


1439 Tlndics=Indice+Nn 

1440 ! 

1450 | Impression des déplacements 

1460 ! 

4470 CALL Depl(lndice,Nd1 ,Nnoeuds ,A(X)) 
1480 ! 

1490 ! Impression des contraintes et des réactions d'appuis 
1500 Û 

1510 CALL Effnilndice,1) 

1520 ! 

1530 PRINT LIN(2) 

1540 NEXT I 

1550 ! 


1560 DISP "FIN PROBLEME. FAIRE CONTINUE POUR UN AUTRE PROBLEME" 
1570 PAUSE 

1580 GOTO 730 

1590 END 


5.5- LECTURE DES DONNÉES 


Les problèmes d’élasticité traités par la méthode des éléments 
finis nécessitent un découpage du domaine en un grand nombre 
d'éléments. Il n’est pas pensable d’introduire la définition des élé- 
ments et les coordonnées des nœuds à l’aide du programme ‘Lect', 
qui est reproduit ici pour des raisons de structure logique du pro- 
gramme. Un pis-aller consiste à utiliser un sous-programme dans le- 
quel les données sont en DATA’ et les variables à initialiser, en 
COM’. En fait, pour un programme destiné à être utilisé fréquem- 
ment par des personnes différentes, il convient d’attacher une grande 
importance à la méthode de description du maillage. Le listing 
du sous-programme ‘Lect', reproduit ci-après, ne fait que définir 
les données nécessaires. 


1600 ! 

1610 1 HBHHNNNNHNMNNNNNNNNONNNNNTENNNNHEPNNNNNNNNHONNNNNNNNNREPHNNONNNNONNNNNNNEEONX 
4620 SUE Lect 

1630 1 RHONE 
1640 ! 

4650 ! Appelé par Prog. Princ. 

1660 ! 

4670 OPTION HASE i 

1680 ! 

4690 ! Lecture des données relatives à l’ossature 

4700 ! 

1740 COM Nnoeuds,Nelem,Nchar ,Nd1 ,Nnelen,Nnodul ,Nepais 
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1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
atet 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2400 
2110 
2420 
2130 
2440 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 


COM Titre$[80],Charge#(5)(801) ,Nat$ 
COM X(400),Y(100),E(S,2) ,Epais(15,2) 
COM INTEGER Liais(200,5) ,Isup(100) 
1 
COM Nd1b,Nn,Ntot,Ipp(5) 
COM INTEGER Lappui(50) ,Maxa(304) 
COM A(S450) ,R(200,15) 
[ 
INPUT "Titre",Titre$ 
Nat4$="C" 
INPUT "Contrainte plane (C) ou déformation plane (D) ?",Nat$ 
INPUT "Nombre de Noeuds" ,Nnoeuds 
INPUT "Nombre d'éléments" ,Nelen 
INPUT "Nombre de cas de charge" ,Nchar 
INPUT "Nombre de d’épaisseurs (1)",Nepais 
INPUT “Nombre de Modules E (1)",Nmodul 
Nn=NnoeudskNd1 
i 
! Lecture des coordonnées des noeuds 
i 
FOR 1=1 TO Nnoeuds 
INPUT "X, Ÿ, ",XC1) ,Y(1) 
INPUT "Code Appui ? (0)=Noeud libre",Isup(I) 
NEXT I 
1 


| Définition des liaisons 

| 

FOR 1=1 TO Nelen 
DISP "Noeud I, Noeud J, Noeud K ?" 
INPUT Liais(1,1),Liais(1,2) ,Liais(1,3) 
Liais(1,4)={ 
Liais(1,5)=1 
DISP "N° Epaisseur (4) N° Module (4) FAIRE CONTINUE Si 4 et 1" 
INPUT Liais(1,4),Liais(1,5) 

NEXT I 

[ 


' Lecture des sections 

( 

FOR 1=1 TO Nepais 
DISP “Epaisseur, Densité N° ";1 
INPUT Epais(1,1) ,Epais(1,2 

NEXT 1 

! 


\ Lecture des modules _ 
1 
FOR I=1 TO Nnodul 
DISP “Module N° ";l 
INPUT EUI,1) 
DISP "Coef. POISSON N° ";I 
INPUT E(I,2) 
NEXT I 
! 
SUBEXIT 
SUBEND 
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5.6 - REPRODUCTION DES DONNÉES 
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Si la reproduction des données est toujours souhaitable, il faudra 
prévoir en plus des contrôles graphiques reproduisant le maillage. 
Le listing du sous-programme ‘’Imprim' est donné ci-dessous : 


2710 
e7e0 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
280t 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3080 
3030 
3040 
3050 
3660 
3070 
3180 
3090 
3100 
3110 
3120 
3436 
3140 


1 POONENNNNNNNNNNNNNNNNONEPNNNNNNNNNONNRNENNNENNNENRNENNINENX 
SUE Imprin 

1 00H OX 
l 


! Appelé par Prog. Princ. 
l 


OPTION BASE 
l 


l Reproduction des données générales de la structure 

| 

CON Nnoeuds ,Nelen,Nchar ,Ndi ,Nnelem,Nmodul ,Nepais 

CON Titre#$(801,Charge$(5)[80],Nat$ 

COM X(100) ,Y(400),E(S,2) ,Epais(15,2) 

COM INTEGER Liais(200,5),Isup(100) 

l 

COM Ndlb,Nn,Ntot,Ipp(5) 

COM INTEGER Iagoui(S0) ,Maxa{301) 

COM A(S150) ,R(200,15) 

j 

PRINT "TITRE : "&Titres 

PRINT | 

IF Nat$="C" THEN PRINT ” Probleme de contrainte plane." 
IF Nat$="D" THEN PRINT " Probleme de déformation plane." 
PRINT 

PRINT "Nombre de noeuds = ";Nnoeuds 

PRINT "Nombre d’éléments= ";Nelen 

PRINT "Nombre de cas de charges = ";Nchar 

PRINT "Nombre d’épaisseurs et densités différentes = ";Nepais 
PRINT "Nombre de modules d’élasticité = ";Nmodul 


PRINT 

PRINT "COORDONNEES DES NOEUDS" 

PRINT 

PRINT " N° X Y Code appui" 
PRINT 


FOR 1={ TO’ Nnoeuds 
PRINT USING 3070,1,X(1),Y(1),Isuptl) 
IMAGE DDD,3X,MDDD.DDD,5X,MDDD.DDD,44X,DDD 
NEXT I 
PRINT 
PRINT "DEFINITION DES ELEMENTS" 
PRINT 
PRINT "Elément I J K N° Sect. N° Module" 
PRINT 
FOR 1=i TO Nelem 


3450 
,5) 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3240 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 


5.7 - MATRICE 'B' ET MATRICE D'ÉLASTICITÉ 'D' 


PRINT USING 3460;1,Liais(1,4),Liais(1,2),Liais(1,3),Liais(1,4),Liais(l 


IMAGE DDD,SX,DDD,6X,DDD,8X,DD,8X,DD,8X,DD 
NEXT I 
PRINT 
PRINT "EPAISSEURS ET DENSITE DES ELEMENTS" 
PRINT 
FOR I=i TO Nepais 
PRINT "Epaisseur N° ";1," 
PRINT "Densité N° ";l," 
PRINT 
NEXT I 
PRINT 
PRINT "MODULES D'ELASTICITE"* 
PRINT 
FOR 1=i TO Nmodul 
PRINT “Module E N° ";l;" 
PRINT “Coef POISSON N° ";1;" 
NEXT I 
SUBEXIT 
SUBEND 


";Epais(I,4) 
“;Epais(1,8) 


“E(1,4) 
*E(1,2) 


CL 


Les coefficients de la matrice ’B' sont calculés à partir des 


coordonnées des nœuds des éléments, et rangés dans la mémoire 
auxiliaire "R'. Les coefficients de la matrice d’élasticité sont calculés 
à partir des caractéristiques physiques de l’élément (module de 
YOUNG, coefficient de POISSON) et de la nature du problème 


traité : contrainte plane ou déformation plane. 


L’organigramme du sous-programme ‘Format’ est reproduit 
figure 71. 
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TisLiais(1,1) 
Jj=Liais(l, 2) 
Kk=Liais(1,3) 
LiLiais(1,4) 
Ma=Liais(1,5) 





Xi=4{1i) 
Yi=Y(15) 
X3=X(J3) 
Y3=Ÿ(J 3) 
Xk=XIKK) 
Yk=YCKK) 


Delta=(Xj-X1)k(Yk-Y12-(Yj-Y1)XK(Xk-X 1) 





RIT, 4)=(Ya-Yk)/Delta 
R(T,2)={Yk-Y1)/Delta 
R(1,3)=(Y1-Y3)/Delta 
R(1,4)={Xk-X7)/Delta 
R(1,5)=(Xi-Xk)/Delta 
R(I,6)=(Xj-X1)/Delta 
R(1,7)=Delta/2 

Eyoung=E(Li,i) 

Xmu=E(L1,2) 













la matrice ‘ER’ 


R(I,7)=R(1,7)XEpais(Mn,1) 
R(1,8)=Eyoung/(i-Xmu”2 











F{1,9)=Eyoung/(1-Xmu) 

(A+Xmu)/(1-2XXMu ) 
R(T,7)=R(I,8)XXmu/ 
(4-Xnu) 


R(I,10)=Eyoung/(1+Xmu)/2 





RCI,9)=RÇI,8)XXmU 


Fig. 71. - Organigramme de Format’. 


Calcul et sauvegarde 
des coefficients de 


Surface de l'élément 


Listing du sous-programme ‘Format’ : 


3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
36410 
3620 
3630 
364 
3650 
3660 
3670 
3680 
3690 
3700 
3740 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 


1 HNONMHNNNMNNHNNNRHNNNNNNNNNRRNNRNNNRIHNNNNHNRNOONINENNNINNNNEONNNNNNX 
SUB Format(1) 

1 HBnNBBHNNHAINNHNMHHHNMHHNNNHNRANNNNHANHNNNHNHONNNNNNNRNHRONNRRNNNRNERNNNONNENNN 
l 


\ Appelé par Prog. Princ. 
[ 


OPTION EASE 

| 

! Calcul des éléments pour la matrice de rigidité 
l 

COM Nnoeuds ,Nelem,Nchar ,Nd1l ,Nnelen,Nnodul ,Nepais 
CON Titre$[801,Charge$(5)[801,Nat$ 
COM X(400),Y(400) ,E(5,2) ,Epais(15,2) 
COM INTEGER Liais(200,5) ,Isup(100) 

l 

COM Nd1b,Nn,Ntot,Ipp(5) 

COM INTEGER Iappui(50) ,Maxa(301) 

COM A(5150) ,R(200,45) 

| 

lisliais(1,1) 

Jj=Liais(l,2) 

Kk=Liais(1,3) 

LisLiais(I,4) 

Mu=Liais(T,5) 

Xi=X(11) 

Yi=Y(1i) 

Xj=X(J j) 

Yj=Y(j) 

Xk=X(KK) 

Yk=YCKK) 
Delta=(Xj-X1)X(YkK-Y1)-CYj-YI2KCXk-X1) 
IF Delta)0 THEN 3690 

PRINT "ERREUR ELEMENT : Surface négative ou nulle ";l 
STOP 

R(1,43=(Yj-Yk)/Delta 
R(I,2)=(Yk-Y1)/Delta 
R(1,3)=(Yi-Yj)/Delta 
R(1,4)=(Xk-Xj)/Delta 
R(I,5)=(Xi-Xk)/Delta 
R(I,6)=(Xj-X1)/Delta 

R(1,7)=Delta/2 

Eyoung=E{Ll,1) 

Mu=E(L1l,2) 

IF Nat$="C" THEN 3880 

1 


| Déformation plane 
il 


R(I,8)=Eyoungk(i-Mu)/(1+Mu)/(1-2XMu) 
R(1,9)=R(1,8)XMu/(1-Mu) 
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3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 


5.8 - MATRICE DE RIGIDITÉ DE L'ÉLÉMENT 


La matrice est calculée globalement à l’aide des fonctions 
matricielles. On reconstitue ’B' et D’ à partir des valeurs sauvegar- 
dées dans ’R’ ; on calcule la matrice de rigidité ’Rr' que l’on met 
sous la forme prévue pour son cumul dans la matrice globale. 


L’organigramme du sous-programme ‘’Matrig’ est donné figure 


T2 
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GOTO 3910 
! 


! Contrainte plane 

| 
R(I,7)=R(1,7)KEpais(Mn,1) 
R(1,8)=Eyoung/(1-MukMu) 
RCI,9)=RCI,8)XMu 
R(1,10)=Eyoung/(1+Mu)/2 

1 


SUBEXIT 
SUREND 









MAT H=(0) 
MAT D=(0) 









H(4,4)2E(3,0)=R(1,4) 
E(4,3)=(3,4)=R(1,2) 
R(4,5)=H(3,6)=R(I,3) 
KH(2,8)=H(3,1)=R(1,4) 
R(2,4)=E(3,38)=R(I,5) 
B(2,6)=B(3,5)=R(1,6) 












D(4,4)=D(2,2)=R(I,8) 
D(4,2)=D(2,4)=R(1,9) 
D(3,3)=R(1,16) 







MAT Tb=TRN(E) 
MAT Aux=DXk 
MAT Rr=TbXAux 


Vol=R(1,7) 
L=Q 











LEL# 
FigéL)=Rr (I ,K)#Vol 





Formation de 
la matrice ‘EH’ 


Formation de la 
matrice ‘D’ 


Matrice transposée 


Volume de l'élément 


Fig. 72. - Organigramme de ‘Matrig’. 
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Le listing correspondant est reproduit ci-dessous : 


4780 
4790 
4800 
4810 
4820 
4830 
4640 
4850 
4860 
4870 
4890 
4896 
49700 
4710 
4900 
4930 
49740 
4950 
4960 
4970 
4980 
4996 
5000 
5016 
Soät 
S030 


5040 


505 
5060 
5070 
5080 
5090 
5100 
5410 
5126 
5130 
5140 
5150 
Siet 
5170 
5160 
5196 
5200 
5210 
Seal 


À HOOONNNHHNNNNNNNNNNNNNNNNNMNENNNNNNNHNNNNNNNNNNNNEEX 
SUR Matrig(R(X) ,Rig(X) ,Ndl,Nnelem, I) 

À KHONHHOPHNNNNPHNPPNNNNNNNENNUNNNNRNHNNNRNNHDNNHNNUNHNNPNNNENNNE 
[l 


Appelé par Assemb 
} 


OPTION ASE 1 
[ 


! Calcul de la matrice de rigidité de la barre sous la forme 
! 

! Rig(i) Rig(3) Rig(é) .... 

! Rig{2} Rig(5) .... 


Rig(4) ..… 

} 

DIM H(3,6),D(3,3),Tb(6,3) ,Aux(3,6) ,Rr(b,6) 
MAT E=(0) 

MAT D=(0) 


l 
E({,4)=8(3,2)=R(1,4) 
B(4,3)=E(3,4)=R(1,2) 
E(4,5)=E(3,6)=R(1,3) 
H(2,2)=R(3,4)=R(I, 4) 
E(2,4)=H(3,3)=R(1,5) 
H(2,6)=K(3,5)=R(1,6) 
1 


D(4,4=D(2,2)=R(1,8) 
Dt,2)=D(2,4)=R(1,9) 
D(3,3)=R(I,40) 

l 


HAT Tb=TRNCE) 
MAT Aux=DXE 
MAT Rr=TbkAux 
! 
Vol=R(I,7) 
! 
L=Q 
FOR J=1 TO NdlkNnelen 
FOR K=J TO 4 STEP -1 
LEL+{ 
Kig{L)=Rr(I,K)XVol 
NEXT K 
NEXT J 
| 
SUBEXIT 
SUKEND 


5.9 - DÉFINITION DES CHARGES 


Comme pour les ossatures, on a prévu des charges sur les nœuds, 
le poids propre, les effets thermiques et les déplacements d’appuis. 


Les différences avec le programme d’ossatures étant minimes, 
nous reproduisons ci-dessous le listing du sous-programme ‘Second’ 
sans en faire les organigrammes. 


5730 
574( 
575 
5760 
5770 
5780 
5790 
5800 
5810 
Sal 
5830 
5840 
5850 
S86i 
5871 
_ 5880 
5890 
5900 
5910 
Syel 
5930 
5740 
5956 
5960 
5970 
5980 
5990 
6090 
6010 
6020 
6030 
6040 
605 
6060 
6070 
et8t 
6070 
6100 
6116 


1 0onmoNnpHNNNHNNNNUMNHNHNRHNNNMNNMNNNNONNNNNNNONNENNNNNNNRNNNONENNNNNNENNIRX 
SUK Second 

{2 000HNPHHNANNNNNNNIEOENNNMNNNONIHNNNNNNNNNNNNHNNONNENNINNNNOOIENNNNNOOEOOOINEONK 
| 

| Appelé par Prog. Princ. 

l Appelle : Generertemp, Genererdepl 

l 


OPTION HASE 1 
! 
1 Introduction du chargement, reproduction des charges 
let constitution du (ou des) seconds membres du système 
l 
COM Nnoeuds ,Nelem,Nchar ,Nd1 ,Nnelen,Nmodul ,Nepais 
COM Titre$1801,Charges$(5) (80) ,Nat$ 
COM X(400) ,Y(100),E(S,2) ,Epais(i5,2) 
COM INTEGER Liais(200,5) ,Isup(100) 
! 
COM Nd1b,Nn,Ntot,Ipp(5) 
COM INTEGER Iappui(50),Maxa(304) 
COM A(S150) ,R(200 ,15) 
! 
DIM F(2),Caef(2),D(200) 
} 
Indice=Maxa(Nn+i)-Nn 
FOR 1=1 TO Nchar 
Ic=0 | Paramètre indiquant les déplacements d'appuis 
MAT D=(0) 
Indice=IndicetNn ! Indice du fer Coef. du vecteur second membre 
PRINT LIN(3) 
PRINT "CAS DE CHARGE N° ";! 


PRINT 

INBUT " Titre du cas de charge ?",Charges#(l) 
PRINT Charges#$(I) 

PRINT 

Index=5 

DISP " NATURE DU CHARGEMENT ?" 

DISP " 4 : Charges aux noeuds" 

DISP " & : Poids propre" 


137 


bi20 
6130 
6140 
bi50 
bibl 
6170 
6186 
6190 
6eût 
6210 
beat 
6230 
634( 
6250 
b2bl 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 
6360 
b370 
6380 
6390 
6400 
64i0 
6480 
6430 
6440 
645 
6465 
6470 
6480 
6490 
8500 
6510 
8520 
6330 
6540 
6550 
6560 
6570 
6560 
6590 
6eût 
6610 
bbël 
6630 
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Ci: 


Ce: 


DISP " 3 : Température" 

DISP " 4 : Déplacement d'appui" 
INPUT " S : FIN" ,Index 

ON Index GOTO Ci,C2,C3,C4 ,Nexti 
1 


! Charges aux noeuds 
| 
PRINT "CHARGES AUX NOEUDS" 
PRINT " N° noeud Fx Fy" 
DISP "La description du chargement est achevée par la donnée" 
DISP "de la charge sur le dernier noeud" 
INPUT " N° noeud, Fx, Fy",Inoeud,F(1),F(2) 
PRINT USING 6250; Inceud,F{4),F(2) 
IMAGE X,DDD,2X,MD.DDDDE ,2X ,MD.DDDDE 
Ji=Indice+Ndlk(Inoeud-1)-1 
FOR K=1 TO Ndl 
A(TitK)=A(Ii+K)+F(K) 
NEXT K 
IF Inoeud=Nnoeuds THEN 6080 
GOTO 6240 ! Introduction de la charge sur le noeud suivant 
l 
! Génération poids propre 
I 
PRINT "POIDS PROPRE (Tous les éléments sont pris en compte)" 
Ipp(T)=i 
FOR L=i TO Nelen 
li=Liais(L,1) 
Jj=Liais(L,2) 
Kk=Liais(Ll,3) 
LisLiais(Ll,5) 
l 
Dens=Epais(Ll,2) 
Po=R(L,7)XDens/3 ! RIL,7)= Volume de l'élément 
l 


li=indicet(li-i)xNd1+i 
Jj=indicet(Jj-1)XNd1+i 
Kk=Indicet(Kk-1)KNdl1+i 
A(li)=4(1i)-Po 
A(Jj)=A(Tj)-Po 
ACKkK)=A(KkK)-Po 

NEXT L 

COTO 6080 

[ 


| Température 
1 


PRINT "TEMPERATURE" 

j 

DISP "N° premier élément, N° dernier élément, Alpha et Température" 
DISF "d'un groupe d’éléments(Cas de charge terminé obligatoirement" 
DISP "avec le numéro du dernier élément de la,structure" 

INPUT Ni, Ne, Alpha, Tenper ?",N4,N8,Alpha,Temper 

PRINT Ni,N2,Alpha,Temper 


664( CALL Generertemp(ni,N2,Alpha,Tenper ,Indice,1) 


665 IF N2=Nelen THEN 6080 

666% GOTO 6620 ! Lecture du groupe d'éléments suivant 
6670 ! 

6680 C4: 1! Déplacements d'appuis 

6690 | 

6700 PRINT "DEPLACEMENTS IMPOSES DES APPUIS"* 

6710 ! 

6720 INPUT " N° du Noeud ?",N 

6736 IF Isup(N)=0 THEN Erreur 

6740 Ic=lc+i 

675û Nombre=Isup(N) 

6760 CALL DecodeiNombre ,Nd1,Coef (4)) 

6770 INPUT "DeltaX, DeltaY ?",F(1),F(2) 

6780 PRINT "Noeud N° ",N," Dx = ";F(4)," Dy = ";F(2) 
6790 Kep$="Oui" 

6800 In=(N-1)KNd1 

6810 | 

6Sel ! Seuls les déplacements imposés correspondants à des degrés de 
6830 ! liberté bloqués sont pris en compte et placés dans le vecteur D 
6840 ! 

6850 FOR Jj=i TO Ndl 

6860 IF Coef(Jj)=0 THEN 6880 

6870 D(IntJj)=F(J 3j) 

6880 NEXT Jj 

6890 INPUT "Autre déplacement d'appui (Oui) ?",Rep$ 
6960 IF Rep$="Gui" THEN 6720 

69710 GOTO 6080 

6920 Erreur: DISP "Le noeud n’est pas un appui” 

6930 GOTO 6720 

6940 Nexti:! 

6950 ! 


6960 ! Modification du second membre en fonction des déplacements d’appuis 
6970 ! 

6980 CALL Genererdepl(Indice,Ic,Nn,Nd1b,D(X),ACK) ,Iappui(k) ,Maxa(X)) 

6990 ! 

7000 NEXT-I 

7040 SUREXIT 

7020 SUBEND 


5,10 - CHARGES NODALES ÉQUIVALENTES AUX EFFETS THERMIQUES 


Les effets thermiques sont des déformations imposées. On 
montre que, dans le cas du problème de contrainte plane, on a : 
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et, dans le cas de déformation plane : 


Ex At 
€y | =(1+v) |aAt 
VYxy 0 


On en déduit les contraintes donnant les mêmes déformations : 
o=[Dle 


et, par application du théorème des travaux virtuels, on obtient les 
charges nodales équivalentes : 


F = Î BD edy 
y 


. En tenant compte des coefficients des matrices ’B' et 'D', on 
trouve : 





en contrainte plane. 


L’organigramme du sous-programme ’Generertemp' est repro- 
duit figure 73. 
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Jisligis(l,1) 
Tiliais(1,2) 
Kk=L13i5(1,3) 

Coef=AlphakTemper 










Coef=Coefx(itE(Liais(1,4),2)) 













R(I,Ichar+1(6)=Coef Sauvegarde 
Coef=Coefk(R(I,8)+R{1,9))KR(1,7 
lislndice+(li-{)#Ndl Adresses 
J j=indicet(Jj-i)XNdl second 
Kk=Indicet(Kk-4)XNd1 nenbre 





ACID=ACIL)4Coef#R (I, 4) 
ACIIH{)=A(Ti#{)+Coef#R(I, 4) 
ACT )=A(T j)+CoefkR(I ,2) 
ACTiHL)=AU jt) +Coef*R (1,5) 
AUKI=A(KID+CoeFXR (I ,3) 
AUKKHLI=ACKK HI +CoeF*R(T, 6) 


Fig. 73. - Organigramme de Generertemp”. 
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Le listing en est donné ci-dessous : 


7030 
7040 
7050 
7060 
7070 
7080 
7090 
7100 
7440 
7420 
7130 
7140 
7450 
7460 
7470 
7480 
7190 
7200 
7240 
7220 
7230 
7240 
7850 
7260 
727 
7280 
7890 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 
7410 


1 oNHPHNNNHNNNHHONNNNNNNNNENENNRANNONNNINNNNNNEENNENONNIONNINX 
SUE Generertemp(Ni,N2,4lpha,Temper ,Indice, Ichar } 

1 oponppHNNNNNRNNNNHNNNNNNNNNHONHNNNNERNNNNHINNNRHNNNNNNNERNNNEINNNTANNNIENX 
1 


| Appelé par Second 
| 


OPTION BASE 1 
[l 


! Genère les charges nodales dues à la température 

l 

COM Nnoeuds ,Nelem,Nchar ,Nd1 ,Nnelem,Nnodul ,Nepais 

COM Titre$(80]1,Charge$(S)[80j,Natt 

COM X(100) ,Y(400),E(S,2) ,Epais(15,2) 

COM INTEGER Liais(200,5),Isup(100) 

1 

COM Nd1b,Nn,Ntot,Ipp(5) 

COM INTEGER Iappui(S0),Maxa(301) 

COM A(5150) ,R(200,45) 

| 

FOR I=N{ TO NE 
li=Liais(I,1) 
Jj=Liais(l,2) 
Kk=Liais(1,3) 
Coef=AlphakTemper 
IF Nat$="C" THEN 7290 
Coef=Coefx(itE(Liais(1,4),2)) 
R(I,Ichar+{0)=Coef 
Coef=Coefx(R(I,8)+R(I,9))KR(1,7) 
li=Indice+(1i-1)XNd1 
Jj=indice+(Jj-1)KNd1 
Kk=Indice+(Kk-1)#Nd1 
A(Ii)=A(Ii)+Coef#R(I,1) 
A(Ii+{)=A(Iiti)+CoefxR(I,4) 
A(Jj)=A(Tj)+CoefkR(I,2 
A(Jj+1)=A(Tj+1)+CoefkR(I,5) 
ACKKI=A(KK)+CoefkRCI,3) 

NEXT I 

SUBEXIT 

SUREND 


5.11 - IMPRESSION DES DÉPLACEMENTS 


Le listing correspondant est reproduit ci-après : 
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9360 1! KoHoNHHNNNNHNNNNPNNNNNNNHNNHIHNNHNNMIHNNNNHNNHNNNNHIENNNNNENNNERX 
9370 SUB Depl(Indice,Nd1 ,Nnoeuds ,A(K)) 
9380 ! XHHNHNNNNNNNMNNNHNNMHNNNNNNNNNHANNNONNHINNHNNNNNNNNNHENNNNEREMNNNRE 
9390 ! 
9400 ! Appelé par Prog. Princ. 

1 


g4a10 ! 

9420 OPTION BASE À 

9430 ! 

9440 ! Impression des déplacements 
9450 ! 


9460 PRINT LIN(i) 

9470 PRINT'DEPLACEMENTS DES NOEUDS" 
9480 PRINT 

9490 PRINT " N° Noeud Dx Dy" 

9500 PRINT 

9510 FOR 1=1 TO Nnoeuds 


9520 li=lndice+(I-1)KkNd1 

9530 PRINT USING 9540;1,A(T1),A(Ii+i) 

9540 IMAGE 3X,DDD,5X,MD.DDDDDDE ,2X,MD .DDDDDDE 
9550 NEXT I 


9560 SUHEXIT 
9570 SUBEND 


5.12 - CALCUL DES CONTRAINTES ET DES RÉACTIONS D'APPUI 


Une boucle sur les éléments permet de calcul les contraintes à 
partir des déplacements. Il faut d’abord reconstruire la matrice "B, 
calculer les déformations, soustraire les éventuelles déformations 
d’origine thermique et, à l’aide de la matrice d’élasticité de l’élément, 
calculer les contraintes. On calcule également les forces nodales 
équivalentes afin de déterminer les réactions d’appuis. 


Outre ces résultats, il est particulièrement intéressant de calculer 
les contraintes principales ainsi que les directions principales. 


L’organigramme du sous-programme ‘Effn' est donné figure 74. 
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Tislndice+(Liais(1,1)-£)xNdl 
Jj=lndice+(Liais(l,2)-1)XNd1 
Kk=Indice+(Liais(1,3)-1)KkNd1 






U(H)=ACT1) 
U(2)=A(Ti+4) 
U(3=ACJ j) 
U(A4)=ACT #4) 
U(S)=ACKK) 
U(B)=ACKK+4) 








EC4,L)=E(3,2)=R(I,4) 
E(1,3)=E(3,4)=R(I,2) 
B(4,5)=EK(3,6)=R(1,3) 
H(2,2)=B(3,4)=R(1, 4) 
R(3,4)=B(3,3)=R(1,9) 
B(2,6)=K(3,5)=R(1,6) 








DC, L)=D(2,2)=R(1,8) 
D(4,2)=D(2,4)=R (1,9) 
Di3,2)=R(I,10) 





MAT S=HX0 
Cste=Ril ,Ichar+10) 
S(4IES(L)-RCT ,Ichar+i0) 
6(2)=8(2)-R(I,Ichar +10) 
MAT Sigma=DxS 
HAT F=TbkSigns 





Sx=Sigma(i) 
Sy=Sigma(e) 
Txy=Sigma(3) 









Fiè=Fii+90 


S1=(Sx+6y)/2-(Sy-Sx)/2KC0S(24F14)+TxyXSINC2KF 14) 
S2=(Sx+Sy)/2-(Sy-Sx)/2HCUS(2XF IS) +TXyXSIN(2F 10) 


Adresses des 
déplacements 
de l'élément 


U vecteur des 
déplacements 
de l'élément 


Formation de 
la matrice ’K’ 


Formation de 
la matrice ‘D’ 


Déformations 
imposées 


Contraintes 
Forces nodales 














Smax=oè 
Snin=Si 
Fi=Fii 













Reac(Jj,K)=Reac(Ij,K)+ 
FCSkCT-L)+K5KRCT,7) 


Tpp(Ichar)=0 ge 
D oui 


Li=Liais (1,5) 
Po=R{l,7;#Epais(Ll,2)/3 
Reac(Jj,K)=Reac(T],K)+Po 











.1,liais(,f),Liais(f, 


HI,Reac(1,1),Reac(l,2) 












1 (T,2),Liais(1,4) 
Bigna(t),61gn8(2) ,5igna(4),Smax Snin,Fi F 





Pas de poids propre 
dans le chargement, 







Fig. 74. - Organigramme 


du sous-programme 'Effn' 
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Le listing du sous-programme 'Effn' : 


9580 1 Kb ONNNEHHNNHNPHMNNNNNNMPNNPNNNNNHNNNNNNNNRHENENNNNONNONNNNNNNEONNEK 
9590 SUR Effn(Indice,Ichar) 
9600 1 KKMNNNNNMNNNMHNNRMANNNNNNNNNEMNNNNNNNNHNNNNONNNNNINOOUNNENNX 
9610 ! 
9620 ! Appelé par Prog. Princ. 
9630 ! 
9640 OPTION BASE 1 
9650 DEG 
9660 ! 
9670 ! Calcul et impression des contraintes dans le éléments 
9680 ! Calcul et impression des réactions d'appuis. 
9690 ! 
9700 COM Nnoeuds ,Nelem,Nchar ,Nd1 ,Nnelem,Nmodul ,Nepais 
97140 COM Titre$[801,Chargest$(S)[80],Nat$ 
9720 COM X(100),Y(400),E(S,2) ,Epais(15,2) 
9730 COM INTEGER Liais(200,S) ,Isup(100) 
9740 ! 
9750 COM Ndib,Nn,Ntot,Ipp(5) 
9760 COM INTEGER Iappui(50),Maxa(304) 
9770 COM A(S150) ,R(200,15) 
9780 ! 
‘ 9790 DIM Coef(e),Sigma(3),B(3,6),Tb(6,3) ,U(6),D(3,4),8(3),F (6) ,Reac(100,2 
9800 MAT Reac=(0) 
9810 ! 
9820 1! Sigma vecteur contraintes 
9830 ! S vecteur déformation 
9840 ! U vecteur déplacements 
9850 ! F vecteur forces nodales 
9860 ! K matrice reliant les déplacements aux déformations 
9870 ! Tb Matrice transposée de K 
9880 ! D matrice d’'élasticité 
9890 ! 
9900 PRINT LIN(4) 
9940 PRINT'CONTRAINTES DANS LES ELEMENTS" 
9920 PRINT 


9930 PRINT " N° Elen I! J K 5x Sy Txy 
Sin Smax Teta” 

9940 PRINT 

9950 FOR 1=1 TO Nelen 

9960 lislndice+(Liais(I,1)-4)KNd1 

9970 Jj=indice+(Liais(I,2)-4)#Nd1 

9980 Kk=Indicet(Liais(I,3)-1)KNd1 

9990 | 

10000 | Déplacements des noeuds de l'élément I 

40040 | 


10020 U(4)=A(T1) 
40030 U(2)=A(Titi) 
10040 UC3)=A(T;j) 
10050 U(4)=A(J j+4) 
10060 U(S)=ACKK) 
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10070. 


10080 
10090 
10100 
10410 
40120 
10130 
10140 
10450 
40160 
10170 
10180 
10190 
10209 
10210 
10220 
10230 
10240 
10250 
10260 
10270 
10280 
10290 
10300 
10310 
10320 
10330 
10340 
10350 
10360 
10370 
10380 
10390 
10400 
10410 
10420 
10430 
10440 
10450 
16460 
10470 
10480 
10490 
10500 
10510 
10520 
10530 
10540 
10550 
10560 
10570 
16580 


U(b)=A(Kk+1) 
[ 


! Matrice K 

! 
B(4,1)=K(3,2)=R(1,1) 
E(1,3)=K(3,4)=R(I,2) 
E(1,5)=K(3,6)=R(I,3) 
K(2,2)=8(3,1)=R(I,4) 
E(2,4)=B(3,3)=R(1,5) 
K(2,6)=B(3,5)=R(I,6) 
1 


MAT Tb=TRNCE) 
l 


\ Matrice d'élasticité 
! 

D(1,1)=D(2,2)=R(1,8) 
D(1,2)=D(2,1)=R(1,9) 
D(3,3)=R(I,10) 

1 


Calcul des déformations 

HAT S=BXU 

Prise en compte des déformations imposées 
S()=8(LI-RCI, Ichar+10) 
S(2)=5(2)-R(I,Ichar+ti0) 

| 


! Calcul des contraintes 
1 


MAT Sigma=DKS 
l 


| Calcul des forces aux noeuds (pour les réactions d’appuis) 
[ 


MAT F=TbtSigna 
1 


Sx=Sigma(i) 
Sy=Sigma(e) 
Txy=Signa(3) 
IF AKS(Sy-Sx)€.000601 THEN 10490 


F14=.SXATN(2XTxy/(Sx-5y)) 
GOTO 10500 


Fii=45 

Fi2=F11+90 

S1=(Sx+Sy)/2-(Sy-Sx) /2KCOS(2kF11)+TxyKSIN(2KF 11) 
S2=(5x+Sy)/2-(5y-Sx)/2KCO0S(2XF12)+TxyXSIN(2YF 12) 
IF S2)=51 THEN 10580 

Snax=5{ 

Snin=52 

Fi=Fie 

GOTO Suite 

IF S2=61 THEN 10630 
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10590 Snax=62 
10600 Snin=5{ 
16610 Fi=Fif 

19620 GOTO Suite 
10630 Snax=Smin=Si 


10640 Fi=0 

10650 Suite:! 

10660 | 

10670 | Calcul des réactions d'appuis 
10680 ! 

10690 FOR J=i TO Nnélen 

10700 Jj=Liais(1,3) 

10710 IF Isup(Jj)=0 THEN 10860 

10720 Nombre=Isup(J;j) 

10730 CALL Decode(Nombre ,Nd1 ,Coef (X)) 
10740 FOR K=1 TO Ndl 

10750 IF Coef(K)=0 THEN 10850 
10760 Reac(Jj,K)=Reac(Jj,K)+F(2K(J-1)4K)XR(I,7) 
40770 ! 

10780 | Prise en compte éventuelle du poids propre 
10790 Û 

10800 IF Ipp(Ichar)=0 THEN 10850 
10810 IF K=1 THEN 10850 

10820 LisLisis(1,5) 

10830 Po=R(I,7)XEpais(Ll,2)/3 
10840 Reac(Jj,K)=Reac(Jj,K)+Po 
10850 NEXT K 


10860 NEXT J 

10870 PRINT USING 10880;1,Liais(1,1) ,Liais(I,2),Liais(]1,3),Sigma({) ,Sigma(2) 
,Signa(3),Smin,Smax,Fi 

10880 IMAGE 4(3X,DDD) ,2X,5(2X,MD.DDDE) ,2X,MDDD .DDD 
10890 NEXT I 

10900 ! 

10910 ! Impression des Réactions d'appuis. 

10920 

10930 PRINT LIN({) 

40940 PRINT 'REACTIONS D'’'APPUIS" 
40950 PRINT 

10960 PRINT " N° Noeud Rx Ry" 

10970 PRINT 

10980 FOR I=i TO Nnoeuds 

10990 IF Isup(1)=0 THEN 11020 

11000 PRINT USING 14040;1,Reac(1,1),Reac(I,2) 
14016 IMAGE 3X,DDD,3X,2(2X,MD.DDDDDDE) 

11020 NEXT I 

41030 SUBEXIT 

11040 SUREND 
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5.13 - EXEMPLE 


Le présent programme a servi pour calculer la répartition des 
contraintes dans un chapiteau. La modélisation comportait 313 
nœuds et 535 éléments. Il a fallu modifier le programme en suppri- 
mant tous les commentaires et la mémoire auxiliaire 'R'. Il est 
intéressant de noter que le nombre total de coefficients de la matrice 
du système linéaire est théoriquement 391876. La méthode de range- 
ment par ‘ligne-de-ciel' n’en conservait que 16567, soit 4.22 %. Le 
chapiteau était en béton et les aciers tendus dans la partie supérieure 
sont pris en compte en donnant des caractéristiques convenables à la 
deuxième bande de triangles (fig. 75). Les données ne sont pas repro- 
duites, et nous présentons uniquement le dessin des contraintes 
principales de compression et des contraintes principales de traction 
(fig. 76 et 77). 
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Fig. 75: - Modélisation. 
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Fig. 76. - Contraintes principales de compression. 

















Fig. 77. Contraintes principales de traction. 
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INERTIE 
DE TORSION 


6.1 - BASES THÉORIQUES 


SAINT-VENANT a montré que le problème de la torsion pure 
se résout à l’aide d’une fonction auxiliaire ® par les équations 
suivantes : 

d* 9? d0 
ee 
à x Ô y dz 








d0 
c-2/| pdx dy =GK— 
dz 


Dans ces relations, G est le module d’élasticité transversal, 
d0/dz l'angle de rotation d’une section par unité de longueur, K 
l’inertie de torsion et C le moment de torsion. 


La fonction Ÿ permet de relier ces grandeurs et d’exprimer les 


composantes de la contrainte de cisaillement (fig. 78). On montre 
de plus que Ÿ doit s’annuler sur le contour de la section 
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Y 
Tx =— 


Tx 79 x 


Fig. 78. - Définition des contraintes. 


La flèche prise par une membrane tendue avec une tension 
constante g sur une ouverture ayant la forme de la section et soumise 
à une pression p obéit à l'équation aux dérivées partielles : 


8?w d?w. p 





SES = — 
Ôx*. 0y? q 


D'où les analogies : 


W — % 
d0 
q dz 


6.2 - EQUILIBRE D'UN PETIT ÉLÉMENT TRIANGULAIRE DE MEMBRANE 


Considérons un petit élément triangulaire dont les sommets sont 
ï, j, k et qui est soumis à une tension q (fig. 79). Cette tension est 
statiquement équivalente à des forces appliquées au milieu des côtés 
ayant pour intensité : 


fkx =Oj—yi)q fry = Xi —x;) q 
fix = Ox — }j) q Ji = (x; — xx) q 


fix =Oi-yx)q fiy = (xx — Xi) q 





Fig. 79. - Forces équivaientes. 


On introduit les notations du paragraphe 5.2 
Les forces équivalentes aux sommets du triangle sont alors : 


l b; l Ci 
Fix na Fiy = 5 Ces + fiy)=Q = 


2 2 

l b; | C; 

F;x date) a Fjy = Ciy + fky) = 4 
l bk 1 Ck 

Fkx 75 xt fx) = Fey 73h + fiy)=4a 


Considérons maintenant le triangle à, j, k après application 
de la pression p (fig. 80). 


Soient W;, W;, Wz les déplacements verticaux des sommets i, j, k. 
Ces déplacements sont supposés être faibles. L'équation du plan 
passant par à’, j’, k' s'exprime par : 


W(z)=@+axt+æ y 
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Fig. 80. - Forces nodales sur un élément triangulaire. 


et les pentes suivant les axes de coordonnées par : 


Ô w Ô w 
— — —_— a 
ax / 0y & 
On montre (voir paragraphe 5.2) : 
; 1 b; b; bg W; 
RE X W; 
2A 
d) Ci €j Ck We 


La projection des forces appliquées au sommet i' sur l’axe Z, 
dans l’hypothèse des petits déplacements s’exprime par : 


Ô w Ô w 
F;=F,;, de an 
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Pour les trois sommets du triangle, on peut écrire : 








F; b; Ci b; b; bx W; 
q 

F; F9 b; Cj °5n Ci Cj Ck W; 

Fx bx Ck We 


En développant cette expression, on obtient la matrice de rigidi- 
té de l’élément triangulaire en équilibre de membrane uniformément 
tendue. 






Les charges nodales s’obtiennent en écrivant que la pression p 
multipliée par la surface du triangle se répartie également entre les 
trois sommets. D’où les charges nodales relatives à l’élément : 








En revenant aux notations utilisées dans le problème de torsion, 
on trouve pour un triangle : 





Il suffit d'exprimer que chaque nœud de la section est en équili- 
bre sous l’action des éléments incidents pour construire le système 
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d’équations linéaires dont la solution fournit les déplacements 
verticaux de la membrane. Il ne reste plus qu’à calculer le volume 
compris entre la membrane et le plan de la section pour obtenir l’i- 
nertie de torsion et les dérivées partielles pour obtenir les composan- 
tes de la contrainte totale de cisaillement. 


6.3 - LE PROGRAMME PRINCIPAL 


Le programme principal se trouve très peu modifié. Les points 
sur le contour de la section sont définis par Isup(l) =1. Il n’y a qu’un 
seul cas de charge (pression p sur la membrane) qui est généré auto- 
matiquement. 


Listing du programme principal : 


1 | DRAM Rennt 
è ! 

40 Lx CALCUL DE L'INERTIE DE TORSION ET * 
40 1x * 
50 Lx DISTRIBUTION DES CONTRAINTES DE CISAILLEMENT * 
60 UX * 
70 Lx XXX OPHPEK OO XOONX * * 

80 Lx X * * X _* X XX X  * 
90 Lx * *X * X * XX OX  * 
100 !* x x XX XXEXX X X *X * 
410 !X * * * * KO OX XX OX 
420 !* * * * X * XX x x 
130 !x xxx * * * KNX *X X _* 
440 !X * 
450 !Xx Lt 
160 1 KKONNNNANNNNNNNHNNNNHNNNNNNNNNNHNONNNNNNONHINNONONNENNNNX 
470 ! 

480. 1! ELOC COMMUN 


190 OPTION BASE { 
200  PRINTER IS 7,1,MIDTH(130) 


210 ! 

200 ! Définition des variables 

230 ©! 

240 ! 1°/ VARIABLES DEFINISSANT LE PROHLEME 

850 ! 

260 ! Nnoeuds........... Nombre de noeuds 

870 ! Nelem............. Nombre d'éléments 

280 ! Nnelem............ Nombre de noeuds par élément 

890 ! Nchar............. Nombre de cas de charges (=) 

300 ! Ndi............... Nombre de degrés de liberté par noeud 
340 ! Nnodul............ Nombre de modules d’élasticité et de coef. de POISSON 
480 1! Titre$............ Titre du probleme 
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330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
456 
460 
A70 
480 
496 
500 
510 
580 
530 
540 
55€ 
560 
570 
seu 
S9t 
600 
610 
CA 
630 
640 
650 
66 
670 
680 


690 


700 
710 
7e 
730 
74( 
7S0 
760 
770 
780 
790 
800 
810 
820 
830 


LOXUL YU). 554% Coordennées du noeud I 

! Isup(l) ..... ...Lode appui 

| = ( Noeud libre (Intérieur) 

i Â Noeud bloqué (Sur le contour extérieur) 


! Liais(1,1)........ Numéro du { er noeud de l'élément 
| Liais(1,2)........ Nunéro du & ème noeud de l'élément 
| Lia1s(1,43)........ Numéro du 3 ème noeud de l'élénent 
| Liais(1,4)........ Numéro du module d’élasticité G 


COM Nnceuds ,Nelen,Nchar ,Nd1 ,Nnelen,Nmodul 
COM Titre$(80] 

COM X(400) ,Y(400),G(5) 

COM INTEGER Liais(200,4) ,Isup(100) 

l 

! 2°/ VARIABLES AUXILIAIRES EN COMMUN 


DeNdIDE es seen es Nombre de degrés de liberté fixés (appuis) 

! Nm. ue ren Noeuds#Ndl : Dimension du système linéaire 
LANTOT a Eaemenan Nombre de coefficients de la matrice À 

li Tappui(l) ........ Numéros des degrés de liberté fixés 

| Maxa(l)........... Adresses des termes diagonaux de la matrice 
AUDE bas éme Tableau des coefficients de la matrice et des 
| seconds membres 

LRO Sa ete 6 pa ones de { à 6 Coef de la matrice rigidité 

l 


| R(1,7) surface de l'élément 
COM Nd1b,Nn,Ntot 
COM INTEGER Iappui(50) ,Maxa(204) 

COM A(S150) ,R(200,7) 

1 

! Initislisation 

} 

MAT Isup=(6) | 

MAT Maxa=(Q) 1! 

MAT A=(0) | CAR VARIAELES EN COMMUN 

MAT R=(U) | 
MAT Liais=(0) ! 

Nnodul=i 
Nnelem=3 

Nchar={ 
Ndl=i 

| 

| Lecture des Gonnées générales et des coordonnées des noeuds 
I 

| CALL Lect 

CALL Donnees 

1 


| Calcul des indices des degrés de liberté bloqués 
l 


CALL IndappNnoeuds ,Ndl ,Nd1b, Isup CH), lappui(#)) 


| Reproduction des données relatives à la structure 
| 
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840 
850 
860 
870 
880 


CALL Imprin 
| 


FOR 1=1 TO Nelen 
1 


! Formation de la matrice de rigidite de la liaison I 
1 


CALL Foraat(l) 
l 


1 Calcul de la hauteur des colonnes correspondantes 
! 


CALL Hcol(I,Ndl,Nnelen,Liais(k),Maxa(x)) 
l 


NEXT 1 


| Calcul des indices des termes diagonaux 
l 


CALL Indiag(Nn,Ntot,Maxa(X)) 
l 


i Assemblage de la matrice globale 
[ 


CALL Assenb 
l 


1 Lecture des charges et création du ou des seconds membres 
| 


CALL Second 
1 


| Modification de la matrice en fonction des appuis 
l 


CALL Modif (Nn,Ndlb,ACX) ,Maxa(X) , Iappui (X)) 


| Résolution du système linéaire 
i 


CALL Solve(Nn,Ntot,Nchar ,ACK) ,Maxa(X)) 


1 
Indice=Maxa(Nn+{)-Nn permet de récuperer les déplacements de chaque 
; cas de charge dans le tableau À. 
PRINT LIN(2) 
PRINT "RESULTATS" 
PRINT 
FOR I={ TO Nchar 
! 
Indice=Indice+Nn 
l 


! Impression des déplacenents 
! 


CALL Depl(indice,Nd1 ,Nnoeuds ,A(X)) 
! 


| Impression des contraintes et des réactions d'appuis 
[ 


CALL EffntIndice, I) 


1360 
1376 
1380 
1390 
41400 
4410 
1420 
1430 


PRINT LIN(2) 
NEXT I 
! 
DISP "FIN PROBLEME FAIRE CONTINUE FOUR UN AUTRE PROBLEME" 
PAUSE 
GOTO 630 


END 
| 


6.4 - LECTURE DES DONNÉES 


Le sous-programme 'Lect' permet de définir le maillage de la 


section. 


Un contrôle graphique sera nécessaire comme pour le pro- 


gramme d’élasticité plane. 


Listing du sous-programme Lect' : 


1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1540 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 


1 BppPHHANNRNRHRHNHNRNNNNNNNNNNNNNRNNNNNNNNRORONNRNRHINNNNNNNNNNNRROONNNNX 


SUB Lect 
L opHponNHNNHEHHHHHNNRHHNRPHANNHHNINRMHNNHNNHHNNNNHANHHPHONNRPANNNHRHRHNRNNHRONNNRX 
! 


! Appelé par Prog. Princ. 
1 


OPTION BASE À 
[ 


|! Lecture des données relatives à l’ossature 
1 
COM Nnoeuds ,Nelen,Nchar ,Nd1 ,Nnelen,Nmodul 
CON Titre$(80] 

COM X(100),Y(100),6(S) 

CON INTEGER Liais(200,4),Isup(100) 

[ 

COM Ndib,Nn,Ntot 

COM INTEGER Iappui(50) ,Maxa(204) 

COM A(5150) ,R(200,7) 

[ 

INPUT "“Titre",Titref$ 

INPUT “Nombre de Noeuds" ,Nnoeuds 

INPUT "Nombre d'éléments" ,Nelen 

Nchar=1 

Nn=Nnoeuds*Nd1 

INPUT “Nonbre de Modules G (1)",Nmodul 

1 


! Lecture des coordonnées des noeuds 
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1700 
1740 
4720 
1730 
4740 
1750 
1760 
1770 
1760 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
{920 
1930 
4940 


L. 
FOR 1=i TO Nnoeuds 

INPUT °X, YŸ, ",XCI),Y(I) 

INPUT "Code Appui ? (0)=Noeud libre",Isup(1) 
NEXT I 
1 


| Définition des liaisons 
| 
FOR I=i TO Nelen 
DISP "Noeud 1, Noeud J, Noeud K ?" 
INPUT Liais(1,4),Liais(1,2) ,Liais(I,3) 
Liais(1,4)=1 
DISP "N° Section (1) N° Module (4) FAIRE CONTINUE 51 i et 1" 
INPUT Liais(1,4) 
NEXT I 
l 


| Lecture des modules 

! 

FOR 1=i TO Nnodul 
DISP “Module N° ";l 
INPUT G(I) 

NEXT 1 

[ 

SUBEXIT 

SUBEND 


6.5- REPRODUCTION DES DONNÉES 


Le sous-programme est adapté à la nature du problème, mais 
ne saurait remplacer un cont-ôle graphique. Le listing est reproduit 


ci-dessous : 
2420 1! HHOBENNONNMHNNNNMHNNNNENNNHNENNNNNNNNNNNNNNNNENNNENNNNNENONNX 
2430 SUB Imprin 
2440 1 KB ONENNNPNONNNHNNNNNNHNNNNNNNNNNNNENNENNNNNNENNNNTNNTONNLE 
2450 ! 
2460 ! Appelé par Prog. Princ. 
2470 ! 
2480 OPTION BASE 1 
2490 ! 
2500 ! Reproduction des données générales de la structure 
2510 ! 
2520 COM Nnoeuds ,Nelem,Nchar ,Nd1 ,Nnelem,Nnodul 
2530 COM Titre$(80] 
2540 COM X(100),Y(400) ,G(5) 
2554 COM INTEGER Liais(200,4),Isup(100) 
2560 ! 
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2570 COM Ndib,Nn,Ntot 

2580 COM INTEGER Iappui(50) ,Maxa(204) 

2590 COM A(S150) ,R(200,7) 

2600 ! 

8610 PRINT "TITRE : "aTitreé 

2620 PRINT 

2640 PRINT "Nombre de noeuds = ";Nnoeuds 

2640 PRINT "Nombre d'éléments= ";Nelen 

2650 PRINT ."Nombre de modules d’élasticité = ";Nmodul 
2660 PRINT 

2670 PRINT “COORDONNEES DES NOEUDS" 

2680 PRINT 

2690 PRINT “ N° X Y Code appui” 
2700 PRINT 

2710 FOR 1=1 TO Nnoeuds 

8720 PRINT USING 2730;1,X(1),YCI) ,Isup(l) 

2730 IMAGE DDD,3X,MDDD.DDD,5X,MDDD.DDD,14X,DDD 
2740 NEXT I 

2750 PRINT 

2760 PRINT "DEFINITION DES ELEMENTS" 

2770 PRINT 


2780 PRINT "Elément ! J K N° Hodule" 

2790 PRINT 

2800 FOR 1=1 TO Nelen 

2810 PRINT USING 2820,1,Liais(1,4),Liais(1,2),Liais(1,3),Liais(I,4) 
2820 IMAGE DDD,5X,DDD,6X,DDD,8X,DD,8X,DD 

2830 NEXT I 

2840 PRINT 


2850 PRINT “MODULES D'ELASTICITE" 

2860 PRINT 

2870 FOR 1=i TO Nnodul | 
2880 PRINT "Module G N° ";1;" = ";G(1) 
2890 NEXT I 

2900 SUBEXIT 

2910 SUBEND 


6.6- FORMATION DE LA MATRICE DE RIGIDITÉ DE L'ÉLÉMENT 


Compte tenu de la simplicité du programme, les six coefficients 
de la demi-matrice sont calculés et rangés dans la mémoire auxiliaire 
’R'. Un contrôle est effectué sur la surface de l’élément, qui ne doit 
pas être ni négative ni nulle. 


L'’organigramme du sous-programme ‘Format’ est donné fi- 
gure 81. 
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N=liais(l.i) 
JisLiais(1,2 
Kk=Liais(1,3) 
Li=Liais(I,4)+ 













Hi=Y(Tj)-YCKk) 
Hj=Y(Kk)-Y(Ti) 
Ek=Y(11)-Y(J 5j) 
Ci=X(KKI-X(T3) 
Cj=X(1i)-X(Kk) 
CKk=X(Jj)-X(11) 


Lelta=(CkEj-CjXEk)/2 












Coef=4xDeltatG(Ll) 
R(I,1)=(Ki*2+C1"2)/Coef 
R(I,8)=(Kj"2+Cj"2)/Coef 

| R(T,H=(BikBj+CiXCj)/Coef | 
R(I,4)=(Hk"2+Ck"2)/Coef | 
RCI,5)=(BjXBk+C]XCk)/Coef | 
R(I,6)=(HiXEk+C1KCk)/Coef 
R(1,7)=Delta 








Fig. 81. - Organigramme du sous-programme Format’. 


Listing du sous-programme ‘Format’ : 
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2910 
evet 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 


1 opNNNNNNNNNNTNNNNNNNNNNNNNNNNNEANNONENENNNNNNNNANNNNNNIENXE 
SUR Format(I) 

MR LLECUESTITE STATUT ES LOST TT TUE TE TT TT ET TT TI TT TT DT TT TT TT TETE ET TI 
\ 


! fppelé par Prog. Princ. 
| 

OPTION BASE i 

| 


!| Calcul des éléments pour la matrice de rigidité 
! 

COM Nnoeuds ,Nelen,Nchar ,Nd1 ,Nnelem,Nmodul 

CON Titre$(80]1 


3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3410 
3i20 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
322 

3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 


COM X(400),Y(100) ,G(5) 
COM INTEGER Liais(200,4),Isup(100) 
} 


COM Ndlb,Nn,Ntot 

CON INTEGER Iappui(50),Maxa(2041) 
COM A(S150) ,R(200,7) 

l 


li=Liais(1,i) 
Jj=Liais(1,2) 
Kk=Liais(1,3) 
LiLiais(1,4) 
Bi=Y(Jj)-Y(Kk) 
Bj=Y(Kk)-Y(Ii) 
Bk=Y(11)-Y(Jj) 
Ci=X(Kk)-X(Jj) 
Cj=X(11)-X(Kk) 
Ck=X(Jj)-X(11) 
Delta=(CkkBj-Cj*#Bk)/2 

IF Delta)0 THEN 3240 
PRINT "ERREUR ELEMENT : Surface négative ou nulle 
STOP 

Coef=4kDeltakG(Ll) 
R(1,1)=(Bi"2+Ci"2)/Coef 
R(1,2)=(Bj"2+Cj"2)/Coef 
R(I,3)=(HikBj+Ci4Cj)/Coef 
R(I,4)=(Bk"2+Ck"2)/Coef 
R(1,5)=(EjkEk+Cj#Ck)/Coef 
R(I,6)=(BiKBk+CikCk)/Coef 
R(1,7)=Delta 

1 


SUBEXIT 
SUBEND 


6.7 - LE SOUS-PROGRAMME 'Matrig' 


‘1 


Il ne fait que récupérer les valeurs des coefficients de la matrice 
de rigidité dans ‘R', et il a été conservé afin de ne pas changer 
l’organisation générale du sous-programme ‘Assemb’. Son listing 
est reproduit ci-dessous : 


4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 


1 1009000000H0pPHNNNHNHHNNHHNNHHNNHHRHIONHNINNNHNNRHNNRHERNNNNHIHNNNRNHHNNNRNHENNNONNENX 


SUK Matrig(R(Kk) ,Rig(*) ,Nd1l,Nnelem,1) 


1 apopHONNHNHNNHNNHNNNNNHNNNRNNNNNENNNNNNNNNONNEIONNONIONENONNENEIONX 


| 
| Appelé par Assemb 
1 


OPTION BASE 4 
l 
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4270 
4280 
4290 
4300 
4310 
4320 
4330 
434( 
4350 
4360 
4370 
4380 


| Calcul de la matrice de rigidité de la barre sous la forme : 
l 
LRig(i) Rig(3) Rig(é) .... 
l Rigte) Rig(S) ... 
Û Rig(4) .... 
l 
FOR J=i TO b 

Rig(J)=R(1,J) 
NEXT J 
| 
SUREXIT 
SUBEND 


6.8 - GÉNÉRATION DU SECOND MEMBRE 


Le second membre est généré à l’aide du sous-programme ’Se- 
cond’ au cours d’une boucle sur les éléments. Le sous-programme 
’Genererdepl' est conservé afin d’annuler les coefficients second 
membre correspondant aux nœuds situés sur le contour. 


Listing du sous-programme ’Second' : 


4890 
4900 
4910 
A9el 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 


164 


RTC TTC TOO TE TT ET ETS TT TT TE CRT ET ET EE TE EE TT TEEN TEL TE CET EEE LE EL 2 1 
SUK Second 

OMR LOL TIR CUT CT TNT AT CT TT TT TT TT LT TITLE TT TT SET ELLES TEL LE TS 
! 

| Appelé par Prog. Princ. 

! Appelle : Genererdepl 

1 


OPTION BASE 4 
| 


| Formation du second membre 

! 
COM Nnoeuds ,Nelen,Nchar ,Nd1 ,Nnelen,Nnodul 
CON Titre$[80) 
COM X(100),Y(400),G(5) 

COM INTEGER Liais(200,4) ,Isup(100) 
| 
COM Nd1b,Nn,Ntot 

COM INTEGER Iappui(S0) ,Maxa(201) 
COM A(S150) ,R(200 ,7) 

! 


DIH D(200) 
l 


Indice=Maxa(Nn+i)-41 
[ 


Si30 
5140 
5150 
5160 
5470 
5180 
5190 
5200 
5210 
See 
5230 
5240 
5250 
5260 
5270 
5280 


FOR I=i TO Nelen 
li=Liais(1,1) 
Jj=Liais(I,2) 
Kk=Liais(I,3) 
Delta=2xR(1,7)/3 
A(Indice+li)=A(Indicetli)+Delta 
A(Indice+Jj)=A(Indice+]j)+Delta 
A(IndicetKk)=A(IndicetKk)+Delta 
NEXT I 
I 


! Modification second membre et sauvegarde charge nodale appuis 
1 


CALL Genererdepl(Indice,Nn,Nd1b,D(#),ACK),Lappui(k),Maxa(*)) 


! 
SUBEXIT 
SUBEND 


6.9 - CALCUL DES CONTRAINTES DE CISAILLEMENT 
ET DE L'INERTIE DE TORSION 


Les composantes des contraintes de cisaillement s'expriment 
en fonction des valeurs trouvées pour le déplacement de la membrane 


par : 





L’organigramme en est donné figure 82. 
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Volume menbrane 
Surface section 


li=indicetLia1s(l,J)-1 


F(J)=ACTi) Déplacements des noeuds 


lisliais(l,4) 
Jjliais(i,8) 
Kk=Liais(l,3) 
Delta=k(i,7)42 
T(1,1)=CX(Kk)-X(J3))/Lelta 
T{4,8)=(X(11)-X(Kk))/Deita 
T(1,3)=(X(Jj)-X(11))/Delta 
Té2,4)=-CYCTj)-YKk))/Delta 
Tié,2)=-{Y(Kk)-Y(11))/Delta 
T(2,2)=-(Y(11)-Y(J3))/Delta 





















MAT Txy=TXkF 


Vol=Voit(F({)+F(2)+F(3))XDeita/b 
Surf=Gurf+Delta/E 
Tmax=SQR(Txy (4) 24Txy(BINE 


lnertie de fürsion 
Surface 


Fig. 82. - Organigramme du sous-programme 'Effn’. 


Le listing est reproduit ci-dessous. 


7230 1 KHAN 
7240 SUK Effn(Indice,Ichar) 

7250! Xk0000nnHnHnNnHHHHHNHNNHNHNRHRRHIUNHNRHNNNNNRENNNNNIEINTNNRIHINNNNNNER 
7260 ! . 

7270 ! Appelé par Prog. Princ. 

7280 1! 

7290 OPTION BASE { 

7300 ! 

7340 ! Calcul et impression des contraintes dans le éléments 

7380 ! 

7330 COM Nnoeuds ,Nelem,Nchar ,Nd1,Nnelen,Nmodul 

7340 CON Titre$[80] 

7350 COM X(100),Y(100),6(5) 

7360 COM INTEGER Liais(200,4),Isup(400) 

7370 ! 

7380 COM Nülb,Nn,Ntot 

7390 COM INTEGER Iappui(S0) ,Maxa(201) 

7400 COM A(5150) ,R(200,7) 

7410 ! 

7420 DIM T(2,3),F(3),Txy(2) 

7430 DEG 

7440 PRINT LIN({) 

7450 PRINT'CONTRAINTES DE CISAILLEMENT" 
7460 PRINT 

7470 PRINT " N° Elen Tx Ty Tnax Fi" 
7480 PRINT 

7490 Vol=0 

7500 Surf=0 

7510 FOR I=1 TO Nelen 


7520 FOR J=1 TO Nnelen 

7530 Ii=IndicetLiais(1,J)-1 
7540 FCJ)=ACI1) 

7550 NEXT J 

7560 Ti=liais(],i) 

7570 Jj=Liais(l,2) 

7580 Kk=Liais(1,3) 

7590 Delta=R(I,7)%*2 

7600 T(4,4)=CX(Kk)-X(Tj))/Delta 
7610 T(4,2)=(X(11)-X(Kk))/Delta 
7620 T(4,3)={XCTj)-X(11))/Delta 
7630 T(2,4)=-(Y(Jj)-Y(Kk))/Delta 
7640 T(,8)=-(Y(Kk)-Y(11))/Delta 
7650 T(2,3)=-(Y{11)-Y(Jj))/Delta 
7660 MAT Txy=TXF 

7670 Vol=Vol+(F(1)+F(2)+F(3))XDelta/b 
7680 Surf=Surf+Delta/2 

7690 Tmax=SQR(Txy (1) "2+Txy(2)*2) 
7700 Fi=90 

7740 IF Txy(1)=0 THEN 7730 
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77e 
7730 
7740 
7750 
7760 
7770 
7780 
7790 
7800 
7810 


Fi=ATNCTxy(2)/Txy(4)) 
PRINT USING 7740,1,1Txy(4),Txy(2),Tmax,Fi 
IMAGE 3X,DDD,3(4X,MD.DDDDE) , 4X,MDDD. DD 
NEXT 1! 
PRINT LIN(2) 
PRINT ” INERTIE À LA TORSION = ";Volx2 
PRINT 
PRINT ” SURFACE. ............. = ";Surf 
SUREXIT 
SUREND 


6.10 - EXEMPLE : SECTION CIRCULAIRE 


La figure 83 représente la modélisation d’une section circulaire 
pleine en 108 éléments triangulaires. 
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Fig. 83. - Section circulaire pleine. 
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TITRE : Test inertie de torsion d’une section circulaire creuse 


Nombre de noeuds = 67 
Nombre d'éléments= 4108 
Nombre de modules d’élasticité = 2 


COORDONNEES DES NOEUDS 


N° X Y Code appui 
Î b.940 .944 Î 
è 6.467 2.679 { 
3 5.553 4.261 î 
4 A .2bi 5.553 { 
5 2.679 6.467 1 
6 .944 6.940 { 
7 -.914 6.940 Î 
8 2.679 6.467 { 
9 -4.261 5.553 { 
10 -6.553 4.2b{ À 


DEFINITION DES ELEMENTS 


Elément I J K N° Module 
{ 24 { 25 { 
2 25 î 26 { 
3 { ë 86 { 
4 26 è 27 Â 
S è 3 27 { 
b 27 3 28 { 
7 3 4 28 { 
8 28 4 29 { 
9 À 5 29 1 
10 29 S 30 î 
100 bb 59 60 è 
104 66 66 bi a 
102 bi 60 49 è 
103 67 bi 62 è 
104 67 62 63 a 
105 67 63 64 è 
106 67 64 65 a 
107 67 65 66 à 
108 b7 66 bi è 


MODULES D'ELASTICITE 


Module G N° 
Module G N° 2 


.00 
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CONTRAINTES 


N° Elem 


NM EUro 


Tx 


0.0000E+00 
-7.4758E-01 
-1.7788E+00 
-2.4014E+00 
—3.3716E+00 


103 -7.2169E-04 
104 -1.4434E-03 
105 -7.2469E-04 


106 
107 
108 


INERTIE 


SURFACE 


7.2169E-04 
1.4434E-03 
7.2169E-04 


À LA TORSION 


DE CISAILLEMENT 


Ty 


6.7438E+00 
6.0458E+00 
6.b384E+00 
5.5947E+00 
5.8399E+00 


1.2500E-03 
4.0000E-11 
1.2500E-03 


-1.2500E-03 


0.0000E+00 
.2500E-03 


re 


3074.25485598 


152.185598519 


Tnax 


6.7432E+00 
6.0882E+00 
b.8726E+00 
6.0883E+00 
b.7433E+00 


1.4434E-03 
1.4434E-03 
1.4434E-03 
1.4434E-03 
1.4434E-03 
1.4434E-03 


90. 
83. 
75: 
66. 
-60. 
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7 ANNEXES 


1.1 - ALTERNATIVE CONCERNANT LES CONDITIONS AUX LIMITES 


L'option choisie pour prendre en compte les conditions d’appuis 
(ou conditions aux limites) a le mérite de la simplicité : on ne modi- 
fie pas la numérotation des degrés de liberté, qui se déduit de l’ordre 
naturel d'introduction des nœuds de la structure. Cette façon de 
faire a néanmoins trois inconvénients. 


Le premier, mineur, provient du fait que la taille du système 
linéaire n’est pas réduite et les zéros placés sur les lignes et les colon- 
nes correspondant aux degrés de liberté sont effectivement traités 
lors de la procédure de résolution. 


Le deuxième a des répercussions sur l’organisation générale 
du programme. 


En effet, on peut distinguer deux phases bien distinctes : la 
première concerne la structure, la deuxième concerne les charges. 
La jonction de ces deux phases se fait au milieu de la procédure 
’Solve’ : la triangularisation peut être menée sans connaître les char- 
ges appliquées. Or, dans la méthode exposée dans l’ouvrage, on a 
besoin des coefficients de la matrice correspondant aux degrés 
de liberté bloqués pour générer les chargements provenant des 
déplacements d’appui, et après l’appel du sous-programme ‘Modif’, 
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ces coefficients sont détruits. Or il est économique d'effectuer 
la factorisation, de conserver la matrice triangularisée en mémoire, 
puis de résoudre ensuite pour autant de chargements que souhaité. 
Cela nécessite un réarrangement des lignes et des colonnes du systè- 
me linéaire en fonction de la nature des degrés de liberté. 


Le troisième inconvénient concerne le calcul des réactions 
d'appui. Le fait d'interroger chaque nœud de chaque élément pour 
savoir si il contribue à une réaction d’appui est une démarche non 
déterministe, peu satisfaisante par les opérations inutiles qu’elle 
engendre. 


On peut donc commencer par renuméroter les degrés de liberté 
suivant qu’ils correspondent à des appuis ou non. La matrice du 
système linéaire est alors partitionnée en quatre sous-matrices, ainsi 
que l’illustrent les figures 85 et 86. 


4 
RE ne, 
LANGE NP 


123456 7 8 910 


Fig. 85. - Numérotation naturelle. 


3456 7 8 1,2 910 





Fig. 86. - Après réarrangement. 
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Seule la sous-matrice A’ est factorisée, et, si on appelle ‘d' le 
vecteur des déplacements imposés, ‘6’ le vecteur des déplacements 
inconnus, ‘F,' et ‘F;' les forces nodales correspondant aux degrés 
de liberté libres et bloqués et "R' les réactions d’appuis, on a, d’une 
manière formelle : 





R='Bô+ Cd—F, 


Les différentes sous-matrices peuvent être soit rangées dans le 
même tableau soit dans deux tableaux différents, l’un réservé à ‘A 
et l’autre à 'B'et 'C. 


7.2 - INDICATIONS SUR LES ESPACES MÉMOIRE NÉCESSAIRES 


Ainsi qu'il a été dit, les programmes présentés ont été mis au 
point sur micro-ordinateur HP 9845, avec un BASIC interprété. 
Les résultats donnés dans ce paragraphe concernant les encombre- 
ments mémoire n’ont qu’une valeur relative. 


Les différents types de variables utilisés nécessitent les espaces 
mémoire suivants exprimés en octets (1 octet = 8 positions binaires) 





Variable simple Octets 
Réél 10 
Entier 4 
Chaîne de caractères 6 + 1 par caractère. 
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Tableaux Octets 





Rééls 10 + 4 par dimension + 8 par nombre 
Entiers 10 + 4 par dimension + 2 par nombre 
Chaînes de caractères 12 + 4 par dimension + 2 par caractère 


+ nombre de caractères de chaque chaîne 


Le tableau (fig. 87) donne l’encombrement mémoire pour les 
quatres programmes présentés. 


PROGRAMME OSSAT POUTR LSTCT ITRSN 











Nombre total de lignes 1086 1119 1104 


Nombre de lignes commentaires 
Taille mémoire (Kilo-Octets) 123.2 | 1106 


Fig. 87. - Taille mémoire et dimensions. 





Commentaires (Kilo-Octets) 
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7.3 - RÉFÉRENCES ET SOURCES 


— Mathématiques 
Pour la résolution du système linéaire, noyau mathématique des 
problèmes traités, nous nous sommes inspiré de l’ouvrage : 


+ Numerical analysis in the finite element method (Prentice Hall), 
des professeurs KLAUSS-JURGEN BATHE et EDWARD L. WILSON. 


En particulier, la numérotation des coefficients de la matrice 
'A' est celle exposée dans cet ouvrage. 


On pourra trouver des développements analogues dans l’ouvrage 
français . 
e Une présentation de la méthode des éléments finis (Maloine Ed.), 
des professeurs GOURI DHATT et GILBERT TOUZOT. 


Ces deux livres traitent également de la génération des éléments 
isoparamétriques et des problèmes dynamiques. Des sous-program- 
mes en FORTRAN sont donnés. 


— Théorie de l’élasticité 


Les rappels concernant la formulation de l’élément triangulaire 
à champ de déplacement linéaire en élasticité plane sont inspirés du 
livre : 


° La méthode des éléments finis (Ediscience), du professeur D. C. 
ZIENKIEWICZ 


— Résistance des Matériaux 


Pour le problème de la torsion et une éventuelle adaptation 
du programme au calcul des sections réduites des poutres, on se 
reportera à l’ouvrage de référence : 


° Résistance des matériaux (Ed. Dunod), de J. COURBON. 


De nombreuses formules relatives aux matrices de rigidité des 
poutres en flexion se trouvent dans l’ouvrage : 


e Analysis of framed structures (Van Nostrand), des professeurs 
JAMES M. GERE et WILLIAM WEAVER Jr. 
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Imprimerie de la Manutention à Mayenne 
Dépôt légal: Mars 1984 
N° d’Éditeur : 4076 





LES DOSSIERS DE LA CONSTRUCTION 


A l'heure où la micro-informatique pénètre les bureaux 
d'études les plus modestes, et où les Ecoles d'Ingénieurs 
enseignent le calcul des structures, cet ouvrage est consacré 
à la présentation d'un programme utilisant la méthode des 
déplacements. 


Les sujets traités (ossatures planes ou spatiales, poutres 
continues, élasticité plane, torsion) ne sont que les différentes 
facettes d’un même problème, et, moyennant quelques modi- 
fications mineures ce programme pourra être utilisé pour le 
calcul des dalles, des portiques, des coques et des problèmes 
de champ en régime permanent : il suffira de modifier les ma- 
trices de rigidité correspondantes. 


La résolution du système linéaire, dont les coefficients 
sont rangés en utilisant la méthode de la ligne-de-ciel, ainsi 
que la diversité des chargements pris en compte confèrent au 
programme un aspect dépassant le cadre des exercices scolaires 
d'initiation. 

Le BASIC, langage de programmation utilisé, est suffisam- 
ment transparent pour permettre une traduction aisée en 
FORTRAN ou en PASCAL. De nombreux commentaires inclus 
dans le programme contribuent à sa bonne compréhension. 





EYROLLES CAL 
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